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Abstract 


ROELLE,  PAUL  ANDREW.  Oxidized  and  Reduced  Biogenic  Nitrogen  Compound 
Emissions  into  the  Rural  Troposphere;  Characterization  and  Modeling.  (Under  the 
direction  of  Viney  Pal  Aneja.) 

Nitrogen  compound  emissions  are  known  to  have  profound  effects  on  air  quality. 
Consequences  associated  with  increased  emissions  of  oxidized  and  reduced  nitrogen 
species  are  known  to  be  increased  tropospheric  ozone  production,  fine  particulate  aerosol 
production,  nitrate  contamination  of  drinking  water,  eutrophication  and  acidification  of 
soil  and  water  bodies.  It  is  well  recognized  that  soil  emissions  can  contribute  a 
substantial  percent  of  the  total  inventory  for  both  the  oxidized  and  reduced  species,  but 
great  uncertainty  still  exists  in  this  inventory.  Using  a  dynamic  flow-through  chamber 
technique  in  conjunction  with  a  state-of-the-art  mobile  laboratory,  this  research  attempts 
to  characterize  and  model  these  oxidized  and  reduced  biogenic  nitrogen  compound 
emissions  into  the  rural  troposphere. 

North  Carolina  has  relatively  recently  witnessed  the  increased  use  of  both 
municipal  waste  biosolids  and  the  land  application  of  swine  waste  effluentj  two  processes 
which  both  contribute  nitrogen  to  the  ecosystem.  The  first  of  these  processes  involves 
the  land  application  of  municipal  waste  biosolids  as  a  cost  effective  way  to  dispose  of 
these  nutrient  rich  byproducts  of  the  wastewater  treatment  process.  During  the  last  three 
decades  extensive  research  has  been  conducted  on  nitric  oxide  emissions  from 
agricultural  soils  and  consequently  an  extensive  database  has  been  developed  which  is 
used  to  relate  these  emissions  to  various  environmental  parameters.  Biosolid  amended 
soils,  however  remain  a  land-use  type  which  are  comparatively  much  less  studied. 


Therefore,  models  used  to  estimate  nitric  oxide  inventories  often  treat  the  biosolid 
amended  soils  the  same  as  agricultural  soils  amended  with  commercially  derived 
fertilizers. 

A  controlled  experiment  involving  the  application  of  municipal  waste  biosolids  to 
agricultural  soils  was  shown  to  enhance  NO  emissions.  A  more  detailed  analysis 
throughout  several  seasons  found  the  nitric  oxide  emissions  from  biosolid  amended  soils 
to  have  a  strong  temperature  dependence  and  that  their  source  strength  is  much  larger 
relative  to  soils  amended  with  chemically  derived  fertilizers.  Modeling  of  this  source 
strength  using  the  MultiScale  Air  Quality  Simulation  Platform  (MAQSIP)  indicated  that 
when  the  biosolids  are  assumed  to  be  spread  evenly  throughout  the  counties,  no  changes 
in  the  model  output  are  evident  during  daylight  hours,  however  it  is  possible  to  discern 
slight  decreases  of  ozone  during  the  evening.  When  the  same  biosolids  are  concentrated 
in  one  area  of  the  county,  as  opposed  to  being  evenly  distributed,  the  changes  are  more 
pronounced  with  decreases  in  ozone  concentration  reaching  as  high  as  approximately 
12%  and  slight  increases  appearing  (approximately  2%)  during  the  afternoon  hours. 

The  second  process  which  is  contributing  nitrogen  to  the  ecosystem  is  the  land 
application  of  swine  effluent  to  agricultural  soils.  This  ammonia-rich  effluent  has  gained 
wide  spread  attention  in  North  Carolina,  due  to  the  explosive  growth  of  the  swine 
industry  during  the  past  decade.  This  study  revealed  that  while  the  average  source 
strength  of  ammonia  from  soils  is  significantly  smaller  than  that  of  the  lagoons,  the  much 
larger  surface  area  of  the  soils  causes  them  to  also  be  an  important  emissions  source. 
Additionally,  it  was  observed  that  for  time  periods  immediately  following  slurry 


application,  the  NH3  flux  increased  significantly  and  remained  elevated  for  at  least  four 
days  after  application.  Temperature  explained  over  70%  of  the  variability  in  NH3 
emissions  prior  to  being  amended  with  the  lagoon  effluent  and  approximately  40%  of  the 
NH3  emissions  for  time  periods  after  being  amended  with  the  lagoon  effluent.  A 
fundamental  meehanistic  mass  transfer  model  is  presented  and  discussed  in  terms  of  its 
applicability  for  estimating  NH3  flux  and  was  found  to  be  an  effective  predictor  of  the 
NH3  emissions  for  time  periods  immediately  following  slurry  application. 


Chapter  I.  Introduction 


Nitrogen  is  perhaps  the  most  important  nutrient  governing  the  growth  and 
reproduction  of  living  organisms.  Nitrogen  compound  emissions  also  have  a  profound 
effect  on  air  quality.  Two  major  needs  that  drive  the  contemporary  perturbations  of  the 
nitrogen  cycle  are  the  seemingly  insatiable  human  appetite  for  energy,  leading  to  the 
emission  of  nitrogen  oxides  into  the  atmosphere,  and  the  need  for  food  to  sustain  growing 
numbers  of  people  all  over  the  world,  leading  to  the  agricultural  emission  of  ammonia. 
Once  released  into  the  atmosphere  by  either  man-made  (anthropogenic)  or  natural 
sources,  these  nitrogen  compounds  can  undergo  several  different  processes  such  as 
transformation  due  to  atmospheric  reactions  (e.g.  gas-to-particle  conversion),  transport 
associated  with  wind,  and  finally  wet  and  dry  deposition  (Figure  1.1).  All  of  these 
processes  can  perturb  the  environment  with  a  host  of  beneficial  and  detrimental  effects, 
such  as  increased  crop  yields  from  nitrogen  loading  or  decreased  visibility  from  increased 
aerosol  production. 

Scientists  have  focused  recently  on  the  oxidized  species  of  nitrogen 
(N0,<=N0+N02)  and  their  role  as  precursors  to  ozone  (O3)  formation,  and  the  reduced 
species  (NH,.=ammonia+ammonium+amines)  and  their  role  in  nitrogen  enrichment  and 
eutrophication  of  aquatic  ecosystems.  Nitrous  oxide  (N2O),  while  contributing  to  ozone 
destruction  in  the  stratosphere  and  to  climate  change  as  a  greenhouse  gas,  is  relatively 
inert  in  the  troposphere  and  therefore  has  negligible  consequences  in  tropospheric 
photochemistry  (Wameck,  2000).  The  concentrations  of  these  trace  species  in  the 


atmosphere  are  small  (0.01  -  400  ppbV),  however  they  have  been  shown  to  have 
significant  impacts  on  chemical  and  physical  processes  in  the  atmosphere  (Conrad,  2001). 

Although  nitrogen  is  a  critical  nutrient  for  the  survival  of  micro-organisms,  plants, 
humans  and  animals,  it  can  cause  detrimental  effects  when  concentrations  reach  excessive 
levels  (Paerl,  1997;  Erisman  et  ah,  1998).  Figure  1.2  (Gunderson  et  ah,  1992)  illustrates 
this  point  by  showing  how  an  ecosystem  responds  to  increased  N  loadings.  The 
horizontal  line  is  a  crop  which  receives  no  atmospheric  N  deposition,  and  as  indicated  by 
the  vertical  axis,  has  a  stable  index  of  productivity.  However,  as  N  is  initially  added  to 
the  system,  the  index  of  productivity  steadily  increases  to  the  point  of  diminishing 
returns,  where  any  additional  N  loading  actually  reduces  productivity  (Schlesinger, 

1997).  In  addition  to  the  productivity  concerns  of  aquatic  and  terrestrial  ecosystems, 
oxidized  and  reduced  N  compounds  each  play  a  specialized  role  in  degrading  human 
health  and  its  welfare.  Some  of  the  consequences  associated  with  elevated  concentrations 
and  depositions  of  both  oxidized  and  reduced  N  species  are: 

1 .  Respiratory  disease  caused  by  exposures  to  high  concentrations  of: 

-  Tropospheric  ozone 

Other  photochemical  oxidants 

-  Fine  particulate  aerosol  (e.g.  PM  2.5) 

Direct  toxicity  of  NO2  (on  rare  occasions) 

2.  Nitrate  contamination  of  drinking  water 
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3.  Eutrophication,  harmful  algal  blooms  and  decreased  surface  water 
quality 

4.  Climatic  changes  associated  with  increases  in  nitrous  oxide 
(greenhouse  gas) 

5.  Nitrogen  saturation  of  forest  soils 
(Erisman  et  al.,  1998). 

Nitric  Oxide 

Nitric  Oxide  (NO),  which  is  estimated  to  have  a  global  source  strength  of 
approximately  44  Tg  yr  '  (Table  1.1),  plays  an  important  role  in  tropospheric 
photochemistry.  Increasing  NO  emissions,  in  the  presence  of  hydrocarbons  and  sunlight, 
are  thought  to  be  the  cause  of  increased  regional  levels  of  tropospheric  ozone  and  other 
photochemical  oxidants  (Logan,  1983;  Penkett,  1988).  Ozone  photochemistry  in  the 
troposphere  is  regulated  by  oxides  of  nitrogen  (NOx  =  NO  +  NO2).  Currently,  the  only 
known  pathway  for  the  production  of  ozone  is  the  photolysis  of  NOj  (NO2  ->  NO  + 
0(^P)),  which  further  reacts  with  O,  to  produce  ozone  (O3)  by  the  reaction  0(^P)  +  O2  -> 
O3.  In  a  pseudo-photostationary  environment,  the  O3  produced  would  react  with  the  NO 
that  was  generated  via  the  photolysis  of  NOj  in  the  following  reaction:  NO  +  O3  — >  NO2  + 
O2.  Hence,  there  is  no  net  production  of  O3.  However,  in  the  real  atmosphere,  hydroxyl 
radicals  combine  with  volatile  organic  compounds  (VOC’s)  to  produce  new  radicals 
which  preferentially  react  with  NO,  allowing  a  net  O3  accumulation.  The  auto  catalytic 
role  of  NO  in  the  cycle  for  ozone  formation  can  be  seen  in  Figure  1 .3. 
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Yienger  and  Levy  (1995)  developed  an  empirically  based  model  to  estimate  soil 
NOx  emissions  on  a  global  scale.  They  have  reported  that  anthropogenic  land  use  is 
having  a  significant  impact  on  global  soil  NOx  emissions  and  that  soil  emissions  can 
account  for  up  to  75%  of  the  total  NOx  budget  depending  on  location  and  time  of  year.  In 
the  Southeast  U.S.,  which  is  NOx  limited,  an  increase  in  NOx  emissions  is  believed  to 
produce  a  corresponding  increase  in  O3  levels  (Southern  Oxidant  Study,  1993).  O3 
negatively  affects  human  health,  as  well  as  ecological  systems,  such  as  crop  yield. 

Studies  show  that  prolonged  exposure  to  high  ozone  levels  causes  persistent  functional 
cheinges  in  the  gas  exchange  region  of  the  lungs.  Additionally,  ozone  plays  a  critical  role 
in  controlling  the  chemical  lifetimes  and  the  reaction  products  of  many  atmospheric 
species  (National  Research  Council,  1991).  Gaseous  nitric  acid  (HNO3),  the  end  product 
of  NO  reactions  in  the  atmosphere,  combines  with  either  aerosols  or  water  in  the 
atmosphere,  and  is  removed  via  rain,  snow,  or  other  deposition  processes,  as  acidic 
deposition. 

Ammonia 

Ammonia  (NH3)  emissions  are  estimated  to  have  a  global  source  strength  of 
approximately  75  Tg  yf'  (Table  1.2).  These  emissions  have  garnered  increased  interest 
in  the  past  few  years,  due  in  part  to  the  detrimental  effects  of  excess  nitrogen  deposition 
to  nutrient  sensitive  ecosystems  (Aneja  et  al.,  1998a;  Nihlgard,  1985;  van  Breemen, 

1982).  Moreover,  NH3  is  the  most  prevalent  gaseous  base  found  in  the  atmosphere,  and 
is  therefore  fundamental  in  determining  the  overall  acidity  of  precipitation  (Wameck, 


4 


2000),  cloudwater  (Li  and  Aneja,  1992),  and  atmospheric  aerosols  (Lefer  et  ah,  1999). 

The  ecological  impact  of  atmospheric  NHj  deposition  may  he  substantial  as  reduced 
nitrogen  species  are  thought  to  be  the  most  biologically  available  of  nitrogen  species  in 
N-limited  coastal  and  estuarine  ecosystems  (Paerl,  1997).  In  the  atmosphere,  NHj  reacts 
primarily  with  acidic  species  to  form  ammonium  sulfate,  ammonium  nitrate  or 
ammonium  chloride,  or  it  may  be  deposited  to  the  earth’s  surface  by  either  dry  or  wet 
deposition  processes. 

The  spatial  scale  of  a  particular  NHj  source’s  contribution  to  atmospheric  nitrogen 
deposition  is  governed  in  part  by  the  gas  to  particle  conversion  rate  of  NH3  to  NH/. 
Because  of  the  short  lifetime  of  NHj  in  the  atmosphere  (t  =  1-5  days  or  less)  (Warneck, 
2000),  low  source  height,  and  relatively  high  dry  deposition  velocity  (Asman  and  van 
Jaarsveld,  1992),  a  substantial  fraction  (20-40%)  will  likely  deposit  near  its  source. 
However,  ammonium  (NH4'^)  aerosols,  with  atmospheric  lifetimes  on  the  order  of  x  =  1- 
15  days  (Aneja  et  al.,  1998b)  will  tend  to  deposit  at  larger  distanees  downwind  of  sources. 
Ammonia  emissions  from  animal  operations  contribute  substantially  to  atmospheric 
nitrogen  loading  and  may  contribute  the  same  order  of  magnitude  as  emissions  of  NO  in 
some  parts  of  the  world  (Steingrover  and  Boxman,  1996). 

Trends 

Over  the  course  of  the  past  century  fossil  fuel  combustion  has  increased  to  meet 
growing  energy  demands.  The  global  amount  of  fossil  fuel  use  per  person  (Figure  1 .4) 
has  increased  by  more  than  a  factor  of  6  over  the  last  75  years.  At  the  same  time. 
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scientists  have  synthesized  nitrogen-based  fertilizers  to  enhance  crop  development  and  to 
maximize  production  on  limited  land  space.  Before  the  mass  production  of  fertilizers,  it 
can  be  assumed  that  there  was  an  approximate  balance  between  the  relatively  unreactive 
moleeular  nitrogen  (N2-comprises  approximately  79%  of  air)  in  the  atmosphere,  which 
was  naturally  converted  to  forms  used  by  plants  and  animals,  and  the  amount  of  nitrogen 
returned  to  the  atmosphere  via  natural  processes  (Delwiche,  1970).  Currently,  however, 
the  global  production  of  fertilizer  is  approximately  100  million  metric  tons  of  nitrogen 
per  year,  compared  to  approximately  one  million  metric  tons  only  40  years  ago  (The 
Fertilizer  Institute,  2000). 

Coal  combustion  was  once  believed  to  be  the  dominant  source  of  ammonia  in  the 
atmosphere  (Warneck,  2000).  Healey  et  al.,  (1970)  showed  that  the  NH3  contribution 
from  coal  was  small  compared  to  NH3  emitted  from  animal  waste  where  the  mechanism 
responsible  for  the  NH3  release  is  the  bacterial  decomposition  of  urea  to  NH3  and  CO2 
(Warneck,  2000).  The  growth  of  the  intensively  managed  animal  operations  in  some 
areas  of  the  country  has  resulted  in  areas  with  elevated  ambient  NH3  concentrations 
(Nelson,  2000).  In  the  US,  for  example,  the  swine  population  accounts  for  approximately 
10%  of  the  NH3-N  budget  (Figure  1 .6a).  In  NC,  however,  which  has  experienced  an 
explosive  growth  in  the  swine  population  (Figure  1.5)  and  is  currently  the  second  largest 
swine  producing  state  in  the  nation,  the  swine  population  is  estimated  to  account  for 
approximately  46%  of  the  NH3-N  emissions  (Figure  1 .6b).  The  combined  effects  of 
increased  fertilizer  use,  intensively  managed  animal  operations  and  increased  power 
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production  have  reached  a  point  where  the  scientific  community  has  major  concerns 
about  the  fate  of  the  nitrogen  produced. 

Nitrogen  Cycle 

Although  the  atmosphere  is  made  up  of  approximately  79%  N,  plants  and  animals 
can  only  benefit  from  that  portion  which  is  first  fixed,  meaning  the  nitrogen  is  converted 
into  a  form  which  plants  or  animals  can  utilize  (Delwiche,  1970).  Referring  to  Figure 
1 .7,  the  nitrogen  present  in  the  atmosphere  can  be  traced  through  the  various  pathways  in 
the  atmosphere-soil  system.  Under  aerobic  conditions,  certain  bacteria  present  in  the  soil 
can  satisfy  their  energy  needs  through  the  oxidation  of  ammonia  to  nitrate,  or 
nitrification.  During  anaerobic  conditions,  the  nitrogen  fixation  process  is  reversed  and 
N2  is  returned  to  the  atmosphere.  As  shown  in  Figure  1 .7,  NH4^  first  appears  in  the  soil 
through  a  process  called  mineralization,  and  then  depending  on  soil  conditions  can  be 
released  as  NHj  gas  to  the  atmosphere.  NO  can  be  released  to  the  atmosphere  through 
both  nitrification  and  denitrification,  each  of  which  will  be  discussed  in  turn  (Warneck, 
2000). 

Nitrification  is  defined  as  an  aerobic  process  in  which  chemoautotrophic  nitrifying 
bacteria  satisfy  their  energy  needs  via  the  oxidation  of  ammonium  (NH/)  to  nitrite  with 
the  end  product  being  nitrate  fMOj')  (Davidson,  1991).  Among  this  chemoautotrophic 
group  of  bacteria  are  the  nitrosomonas  which  oxidize  NH4^  to  NO2’  and  nitrobacter  which 
oxidize  N02‘  to  NOj'  (Delwiche,  1970).  While  not  apparent  in  Figure  1.7,  several 
researchers  have  used  laboratory  experiments  to  show  that  chemoautotrophic  nitrifiers 
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can  produce  NO  in  the  soil  (Firestone  and  Davidson,  1989).  These  studies  have  consisted 
of  acetylene  or  nitrapyrin  being  used  to  inhibit  NH4^  oxidation  and  ehlorate  to  inhibit 
N02'  oxidation  which  have  shown  that  NO  is  produced  as  a  result  of  the  organisms 
responsible  for  the  oxidation  of  NH/  to  NO2'  (Davidson,  1992;  Tortoso  and  Hutchinson, 
1990;  Holbrook,  1994). 

Referring  to  Figure  1.7,  dentrification  is  defined  as  the  reduction  of  NO3’  to  NO2’ 
and  further  to  N2.  This  process,  which  occurs  when  the  soil  is  almost,  but  not  completely 
anaerobic,  serves  as  the  N2  regeneration  portion  of  the  N  cycle  (Warneck,  2000).  Unlike 
the  relatively  fewer  types  of  bacteria  responsible  for  nitrification,  there  are  several 
different  types  of  bacteria  which  carry  out  the  dentrification  process.  These  denitrifiers, 
pseudomonas  denitrificanas  for  example,  are  aerobic  species  which  have  the  ability  to  use 
the  oxides  of  nitrogen  as  an  O2  source  (Delwiche,  1970).  Anaerobic  conditions  in  the  soil 
are  typically  found  to  be  the  result  of  increased  soil  moisture,  or  waterlogging  of  the  soil 
pores  (Warneck,  2000).  The  water-logging  of  the  soil  pores  acts  to  restrict  the  diffusion 
of  gases  through  the  soil  and  to  increase  the  time  for  NO  to  diffuse  to  the  soil  surface. 

The  slowed  transfer  results  in  further  reduction  and  hence  less  NO  is  emitted  from  the 
soil.  The  effect  of  decreased  NO  emissions  as  %  WFPS  reaches  a  critical  value  can  be 
seen  in  Figure  1.8  which  depicts  Davidson’s  (1991)  model  of  N  production  versus 
different  values  of  %  WFPS.  Although  Figure  1.7  does  show  NO  produced  via 
denitrification,  it  is  not  commonly  found  to  be  a  large  contribution  to  the  total  flux  in 
field  studies.  Most  studies,  in  fact,  find  the  chemoautotrophic  nitrifying  bacteria  to  be  the 
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predominant  source  of  the  NO  emitted  from  the  soils  (Anderson  and  Levine,  1986; 
Tortoso  and  Hutchinson,  1990;  Hutchinson  et  al.,  1993). 

Methods  and  Materials 
Dynamic  Flow-Through  Chamber 

A  dynamic  flow-through  chamber  lined  with  five-mil  thick  fluorinated  ethylene 
propylene  (FEP)  Teflon  was  used  to  measure  NO  and/or  NH3  concentrations  emitted  from 
the  soil.  The  translucent  chamber,  26.5  cm  in  diameter,  and  47.2  cm  high  (volume=26.0 
liters),  fits  inside  a  stainless  steel  metal  ring,  which  is  driven  into  the  ground  to  a  depth  of 
~10  cm  (Figure  1 .9).  Zero  grade  air,  which  is  used  as  a  carrier  gas,  is  passed  through  the 
chamber  at  a  constant  flow  rate  (approximately  5  1pm),  via  a  flow  controller  located  on 
top  of  the  chamber  (Gilmont  Shielded  Industrial  Flowmeter,  Accuracy  ±  5%).  The  air 
inside  the  chamber  is  mixed  by  a  variable-speed,  motor  driven  Teflon  impeller.  The 
sample  exiting  the  chamber  travels  through  Teflon  tubing  (1/4”  outside  diameter,  5/32” 
inside  diameter)  to  the  detection  instruments.  The  entire  measuring  system,  from  the  inlet 
port  on  the  chamber  to  the  point  where  the  stream  is  analyzed  in  the  instrument,  is  coated 
by  Teflon,  stainless  steel  or  gold  to  minimize  further  chemical  reactions  with  the  sample 
stream.  The  sample  lines  do  not  exceed  10  meters.  The  NO  detection  instruments  draw 
approximately  1.5  1pm,  which  resulted  in  a  sample  residence  time  in  the  sample  lines  of 
approximately  5  seconds. 
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Chamber  Effects 

The  techniques  currently  employed  to  measure  fluxes  of  trace  gases  from  soil  and 
water  surfaces  can  be  divided  into  micrometeorological  and  chamber  techniques.  The 
micrometeorological  techniques  consist  of  the  gradient  and  eddy  correlation  methods 
while  the  chamber  techniques  can  be  divided  in  to  the  static  and  dynamic  methods.  Each 
system  has  its  benefits  and  drawbacks.  However  as  part  of  National  Science  Foundation 
(ATM- 94206 10)  “NOx  Emissions  from  Agricultural  Soils:  Intercomparison  of 
Micrometeorological  and  Dynamic  Chamber  Techniques”  project,  these  differences  were 
analyzed  and  determined  not  to  be  statistically  significant  (Li  et  al.,  1999;  Roelle  et  al., 
1999).  Nevertheless  a  major  concern  with  the  chamber  method  is  that  it  is  intrusive  and 
could  possibly  alter  the  microclimate  and  therefore  is  not  representative  of  the  real 
conditions.  To  as  great  an  extent  as  possible,  this  concern  has  been  examined. 

Two  of  the  principle  concerns  with  the  chamber  design  are  increased  temperature 
(“greenhouse  effect”)  and  increased  humidity  in  the  chamber  and  sample  lines. 
Experiments  were  conducted  by  previous  researchers  and  as  part  of  this  study.  Sullivan 
(1996)  reported  that  the  average  temperature  difference  of  the  soil  inside  and  outside  the 
chamber  was  0.23±1.0°C  and  that  air  temperature  difference  inside  and  outside  the 
chamber  was  3.77  ±  2.52  °C.  The  average  of  five  experiments  at  various  research  sites 
were  conducted  to  verify  previous  results  and  found  that  the  average  temperature 
difference  between  the  soil  inside  and  the  soil  outside  the  chamber  was  0.03  ±  0.65  °C 
(Figure  1.10).  The  days  selected  for  these  experiments  were  during  the  summer  to 
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determine  the  temperature  difference  at  maximum  heating  and  all  days  were  either  clear 
or  partly  cloudy  with  daytime  highs  of  32  °C  ±  3  °C. 

Excess  humidity  in  the  chamber  and  sample  lines  is  a  problem  when  using 
ambient  air  as  the  carrier  gas.  This  problem  was  addressed  by  periodically  blowing  out 
any  condensation  in  the  lines  and  chamber  with  zero  grade  air,  and/or  creating  a  water 
trap  to  help  collect  any  moisture  prior  to  the  carrier  gas  entering  the  chamber.  During  all 
measurement  campaigns  presented  here,  the  carrier  gas  used  was  zero  grade  air,  which 
virtually  eliminated  the  issue  of  moisture  in  the  chamber  and  sample  lines.  Occasionally, 
if  the  soil  is  sufficiently  moist,  there  was  still  moisture  that  developed  every  few  hours  in 
the  sample  lines  even  with  the  zero  grade  air.  To  analyze  the  effect  of  moisture,  the 
analyzer  was  spanned  with  known  concentrations  of  NO  and  NH3  with  a  short  (<2  m) 
section  of  Teflon  tubing  and  then  this  piece  of  tubing  was  replaced  with  a  10  m  length  of 
tubing  containing  visible  moisture  droplets.  There  was  no  detectable  loss  of  either  NO  or 
NH3  after  40  minutes  of  sampling.  In  another  study,  the  chamber  lining  was  wetted  so 
that  there  were  visible  droplets  on  the  chamber  walls  and  a  petri  dish  filled  with  water 
was  placed  on  the  bottom  of  a  sealed  chamber.  A  known  concentration  of  NH3  was 
delivered  to  the  chamber  and  losses  in  the  chamber  ranged  from  6-7%. 

Experiments  were  conducted  to  determine  if  the  mixing  speed  of  the  Teflon 
impeller  altered  soil  NO  flux  measurements.  Varying  the  speed  between  20  and  100 
revolutions  per  minute  (rpm)  did  not  produce  any  significant  changes  in  the  calculated 
NO  flux.  The  impeller  was  set  to  50  rpm  for  the  remainder  of  the  experiment.  Outlets  in 
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the  chamber  ensured  that  there  were  no  substantial  pressure  differences  between  the 
outside  atmosphere  and  the  air  within  the  chamber.  Research  conducted  on  similar 
chambers  using  a  tilting  water  manometer  indicate  that  pressure  differences  were  below 
detection  limits  (0.2  mm  H2O)  (Johansson  and  Granat,  1984). 

The  dynamic  chamber  system,  with  the  continuous  stirring  provided  by  the 
impeller,  meets  the  necessary  criteria  for  performance  as  a  Continuously  Stirred  Tank 
Reactor  (CSTR).  For  performance  as  a  CSTR,  the  chamber  needs  to  be  “ideally”  mixed 
(Aneja,  1976).  In  ideal  mixing,  the  composition  of  any  elemental  volume  within  the 
chamber  is  assumed  to  be  the  same  as  that  of  any  other  volume  within  the  chamber. 

Trace  experiments  (Residence  Time  Distribution)  were  used  to  test  the  flow  and  mixing 
characteristics  of  the  system.  The  results  of  these  mixing  studies  indicated  that  the 
dynamic  chamber  behaved  as  a  “perfect”  mixer  with  negligible  stagnancy  or  channeling 
(Aneja  et  al.,  2000).  NHj  fluxes  were  calculated  using  a  mass  balance  approach  as 
described  in  Roelle  et  al.  (1999)  and  Aneja  et  al.  (2000). 

Utilizing  the  power  law  profile,  which  is  frequently  used  in  air  pollution 
applications  (Arya,  1999),  we  are  able  to  estimate  the  wind  velocities  at  a  height  of  0.1  m 
(the  approximate  height  of  the  impeller  above  the  soil  air  interface),  when  the  10  m  wind 
heights  are  known.  The  power  law  profile  is  given  by: 
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where  is  the  wind  velocity  at  a  reference  height  and  m  is  taken  to  be  0. 1  for  this 
smooth  surface  (Arya,  1988).  Mean  wind  velocities  were  between  1  and  4  m  s’'  (at  a 
height  of  10  m)  throughout  this  measurement  period.  Using  the  power  law  profile 
(Equation  1),  and  the  measured  10  m  wind  velocities,  equates  to  wind  speeds  between  0.6 
and  2.4  m  s"'  at  a  height  if  10  cm  (approximate  height  of  impeller).  This  calculated  10  cm 
wind  speed  approximates  the  wind  speeds  measured  inside  of  the  chamber 
(~  1  to  2.5  m  s"’)  using  a  hot  wire  anemometer.  However,  it  should  be  noted  that  equality 
of  the  horizontal  winds  does  not  ensure  equality  of  the  emissions  between  the  inside  and 
outside  of  the  chamber. 

While  the  different  flux  methodologies  have  been  intercompared  and  found  to 
produce  reasonable  agreement  in  emission  estimates,  some  comparison  of  the  factors 
effecting  diffusivity  should  be  examined.  For  example,  in  the  neutral  surface  layer  the 
mixing  length  and  eddy  viscosity  can  be  shown  to  vary  linearly  with  height  and  depend 
on  the  friction  velocity  as  described  in  Arya,  1999.  Evaluation  of  this  parameter  was 
outside  the  scope  of  this  research,  however  future  work  to  verify  equality  between 
conditions  inside  and  outside  the  chamber  should  attempt  to  specify  eddy  length  scales 
and  viscosity. 

Sampling  Scheme 

The  daily  sampling  scheme  consisted  of  measuring  concentrations  of  NO  and/or 
NHj  after  the  sample  exited  the  dynamic  flow-through  chamber.  A  daily  experiment 
consisted  of  placing  the  chamber  on  the  stainless  steel  collar,  which  had  been  inserted 
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into  the  soil  the  previous  evening.  Typically,  all  experiments  were  conducted  a  minimum 
of  50  meters  from  any  field  edges  to  avoid  any  boundary  effects.  The  chamber  was 
placed  on  the  collar  at  approximately  8:00  AM  and  flushed  with  zero  grade  air  for  at  least 
one  hour  before  data  collection  began  at  9:00  AM.  This  technique  ensured  that  the 
concentrations  within  the  chamber  reached  steady  state  prior  to  any  data  acquisition  and 
allowed  for  the  instruments  to  undergo  their  daily  calibrations.  Daily  experiments  ended 
at  approximately  5:00  PM  and  the  stainless  steel  collar  was  relocated  to  a  random 
location  within  a  10m  radius  of  the  mobile  laboratory,  in  preparation  for  the  next  days 
experiment.  This  procedure  allowed  a  minimum  of  16  hours  for  any  effect  on  soil  NO 
flux,  due  to  soil  disturbances  caused  by  the  insertion  of  the  stainless  steel  collar,  to 
dissipate.  Depending  on  research  objectives,  24  hour  diurnal  experiments  were,  on 
occasion,  also  conducted. 

Nitric  Oxide  and  Ammonia  Analyzers 

Nitric  Oxide  (NO)  concentrations  were  measured  using  a  Thermo  Environmental 
Instruments  Incorporated  (TECO)  Model  42S  chemiluminescence,  low  level  NO  analyzer 
(Thermo  Environmental  Instruments,  Inc.,  1992).  The  principle  behind  the  operation  of 
the  NO  instrument  is  the  gas  phase  reaction  between  ozone  (O3)  and  NO,  given  by: 

NO+  Oj^NOj+Oj+hv  (2) 

Light  emissions  from  the  decay  of  NOj  to  lower  energy  states  is  proportional  to  the 
concentration  of  nitric  oxide.  This  decay  is  detected  in  a  photomultiplier  tube  and 
converted  to  a  concentration  measurement  after  calibration  with  known  standards. 
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The  TECO  model  42S,  in  its  standard  eonfiguration,  does  not  measure  NH3. 
Therefore,  the  instrument  had  to  he  modified  so  that  we  could  also  calculate  the  ammonia 
emissions  from  agricultural  soils.  The  modified  instrument  is  still  based  on  the  gas  phase 
reaction  between  O3  and  NO,  and  the  detection  of  the  decay  of  the  excited  NO2  molecule 
as  previously  described  (a  schematic  of  the  modified  42S  is  given  in  Figure  1.11).  In 
flow  path  1,  the  sample  stream  is  mixed  with  high  concentrations  of  O3  to  produce  an  NO 
reading.  In  flow  path  2,  the  sample  stream  bypasses  the  stainless  steel  converter,  and 
passes  through  the  molybdenum  converter  (325  °C)  which  effectively  converts  the  oxides 
of  N  to  NO,  while  allowing  the  NH3  to  pass  through,  producing  an  NOy  reading.  In  flow 
path  3,  the  sample  passes  through  the  stainless  steel  converter  (775  °C),  effectively 
converting  the  oxides  of  N  and  NHj  to  N,  producing  a  value  referred  to  as  Nt.  The 
instrument  is  designed  to  alternate  between  flow  paths,  and  by  subtracting  the  NOy  signal 
from  Ny,  the  concentration  of  NH3  in  the  sample  stream  can  be  determined. 

Calibrations  of  the  instruments  were  conducted  following  protocols  using  a  TECO 
146  dilution/titration  instrument  in  conjunction  with  a  cylinder  of  3.82  ppmV  NO  in  N2 
(National/Specialty  gases)  and  zero  grade  air  (National  Welders).  The  TECO  146  was 
serviced  and  calibrated  to  specifications  by  the  manufacturer  prior  to  the  Summer  1 999 
measurement  campaign.  A  multipoint  calibration  was  conducted  prior  to,  and  at  the 
midpoint  of,  the  measurement  period.  Each  day,  zero  and  span  checks  were  conducted 
according  to  the  operator  manuals. 
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Meteorological  data  were  collected  via  a  Campbell  Scientific  meteorological 
package.  Air  Temperature  and  %  Relative  Humidity  (%RH)  are  measured  inside  of  a 
radiation  shield  at  a  height  of  1 .5  m.  The  accuracy  of  both  the  RH  and  temperature  probe 
is  (+)  3%.  A  Ll-Cor  200SZ  pyranometer  is  used  for  measuring  sun  plus  sky  radiation 
(accuracy  ±3%). 

Data  Acquisition  System  and  Mobile  Laboratory 

All  of  the  gas  detection  and  data  acquisition  instrumentation  is  housed  in  a 
temperature  controlled  mobile  laboratory.  The  mobile  system  consists  of  a  modified  Ford 
Aerostar  van  with  a  13,500  BTU  air  conditioning  unit.  The  temperature  inside  the  van  is 
maintained  within  the  operating  range  of  the  instruments.  The  data  acquisition  system  is 
a  Campbell  Scientific  21X  Micrologger  used  in  conjunction  with  a  Toshiba  laptop 
computer.  Power  for  the  air  conditioning  and  all  of  the  detection  instruments  was 
standard  1 1 0  volt  AC  commercial  power. 
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Flux  Calculation 


The  NO  fluxes  were  calculated  from  a  mass  balance  equation  (Kaplan  et  al.,  1988; 
Kim  et  al.,  1994).  The  mass  balance  equation  is: 

Accumulation  =  Flow  in  -  Flow  out  -  Reactions 


where 


dC  ^[Ca;>] 
V 


p-  V  V 


(3) 


J  =  emission  flux  per  unit  area 

L  =  loss  term  by  chamber  wall  per  unit  area  assumed  first  order  in  [NO] 
q  =  flow  rate  through  the  chamber 
V  =  volume  of  the  chamber 
C  =  NO  concentration  in  the  chamber 

Cajr=  NO  concentration  in  the  ambient  air  immediately  adjacent  to  the  chamber 
(the  inlet  of  the  chamber) 

R  =  chemical  production/destruction  rate  in  the  chamber 
A  =  surface  area  covered  by  the  chamber 
A'  =  surface  area  of  the  chamber  walls 

Assuming  the  chamber  is  well  mixed,  the  concentration  [C]  measured  can  be 
assumed  to  be  the  same  everywhere  within  the  chamber.  Using  zero  grade  air  as  the 
carrier  gas,  the  initial  concentration  or  [C^,>]  is  0.  Additionally,  at  steady  state  conditions, 
the  change  of  concentration,  with  respect  to  time,  will  be  zero.  Equation  (3)  reduces  to: 
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(4) 


JA  LA' 
V  V 


where  [c]  is  the  concentration  measured  at  the  outlet  of  the  chamber. 

In  equation  (4),  the  new  loss  term  (Z')  is  the  sum  of  the  loss  of  NO  through 

reactions  with  the  chamber  walls  and  chemical  reactions  of  NO  with  existing  oxidants  in 

the  carrier  gas  (r),  such  as  ozone  and  peroxy  radicals  (Kim  et  ah,  1994;  Aneja  et  ah, 

1995).  The  total  loss  term(Z')  was  determined  empirically  (five  experiments  were 

conducted  throughout  the  day  and  night)  utilizing  a  method  developed  by  Kaplan  et  al. 

(1988).  Kaplan’s  technique  uses  the  modified  form  of  equation  (3): 

dC  JA  f  LA'  ^ 
dt~  V  I  F  Fy 

and  introduces  an  equilibrium  value  for  NO,  given  by  .  Utilizing  the  equilibrium 


JA 

value,  and  plugging  in  for  —  yields: 


dt 


LA'  q 
«/  V 


=(c„-c 


(6) 


The  solution  to  (6)  is  given  by: 


-In- 


Ceq  —  C 


Ccq-Co 


^LA' 

- +  — 

V  F  F. 


(7) 


Co  is  the  NO  concentration  in  the  chamber  when  NO  reaches  the  first  equilibrium 
state  at  an  initial  flow  rate  and  is  the  NO  concentration  in  the  chamber  after  the  flow 
rate  is  reduced  and  allowed  to  reach  a  second  equilibrium.  From  the  linear  relationship 
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Ce  —C 

between  the  value  of  -  In— - and  time  during  the  experiment,  the  slope  is  found  to 

Ceq  ~  Co 


represent 


^  ^  +— 1 .  The  graphs  of  -  In-^^^^^ — ^plotted  versus  time  for  the  5 


V  V 


C  -C 

^ecj 


different  Kaplan  experiments  is  shown  in  Figure  1.12.  The  total  loss  in  the  chamber  was 
estimated  to  be  0.02  ±  0.007  cm  second  '  from  the  average  of  5  experiments  conducted  at 


different  times  throughout  the  measurement  period.  This  value  of  (Z')  (=.02  cm  second  ') 
agrees  with  that  found  by  Kim  et  al.  (1994),  and  is  directly  used  in  equation  (4)  to 
calculate  the  NO  flux  during  this  measurement  period. 

Soil  Analysis 

Soil  temperature  was  recorded  every  minute,  and  these  values  were  binned  and 
averaged  every  1 5  minutes  using  a  Campbell  Scientific  soil  temperature  probe  (accuracy 
±  3%)  inserted  5  cm  into  the  soil,  adjacent  to  the  chamber.  A  soil  sample  was  taken  from 
the  center  of  the  dynamic  flow-through  chamber  footprint  at  the  end  of  each  measurement 
period  (approximately  1  sample  per  day)  and  analyzed  for  soil  pH,  and  percent  Water 
Filled  Pore  Space  (%WFPS)  by  the  Duke  University  Department  of  Civil  and 
Environmental  Engineering  and/or  the  North  Carolina  State  University  Department  of 
Biological  and  Agricultural  Engineering.  Samples  were  taken  with  a  bucket  auger  which 
removed  a  soil  core  to  a  depth  of  approximately  20  cm. 

Soil  bulk  density,  which  is  the  weight  of  the  soil  solids  per  unit  volume  of  total 
soil,  and  soil  particle  density,  are  used  to  determine  the  %WFPS  of  the  soil  (Troeh  and 
Thompson,  1993).  The  core  method  (345  cm^)  was  used  to  determine  the  soil  bulk 
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density  for  the  research  site  (Blake  and  Harge,  1986).  The  standard  particle  density  is 
2.65  g/cm^  for  most  soils,  however  particle  densities  will  differ  from  this  value  if  the  soils 
have  high  organic  matter  content  or  are  high  in  heavy  minerals  such  as  hydrous  oxides  of 
iron.  The  expression,  %  WFPS,  is  a  measure  of  soil  water  content  and  can  be  expressed 
as  the  percentage  of  pore  spaces  in  the  soil  filled  with  water.  The  %WFPS  is  a 
convenient  expression  to  describe  soil  moisture  because  it  accounts  for  the  differing  bulk 
and  particle  densities  of  soils  and  therefore  allows  for  the  comparison  of  soil  moisture 
from  different  soil  types. 

The  total  extractable  nitrogen  was  calculated  by  summing  the  extractable  fractions 
of  ammonium  (NH4'^)  and  nitrate  (N03')  determined  from  the  soil  samples.  Extractable 
NH4''  and  NOj'  were  determined  using  a  1  M  KCL  soil  extract  (expressed  on  a  weight 
basis)  (Keeney  and  Nelson,  1982)  and  standard  autoanalyzer  techniques  (Lachat 
Instruments,  1990).  The  Total  Kjeldahl  Nitrogen  (TKN)  was  determined  using  a 
digestion  procedure,  which  converts  all  the  N  in  the  sample  to  NH4'^  whose  concentration 
is  then  determined  using  colorimetry.  The  total  soil  water  content  was  calculated  as: 

[initial  weight  -  oven  dry  weight  ( 1 05  °C)]  /  oven  dry  weight  (105  °C)  (8) 
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Objectives 


The  objectives  of  this  study  can  be  stated  as  follows: 

1 .  Design  a  field  study  to  measure  nitric  oxide  emissions  from  biosolid  amended 
soils. 

2.  Evaluate  nitric  oxide  emissions  from  soils  amended  with  biosolids  and 
develop  an  algorithm  for  predicting  its  source  strength  in  North  Carolina. 

3.  Include  biosolid  amended  soils  as  a  new  land-use  class  in  biogenic  NO 
emission  inventories.  Compare  and  contrast  inventories  with  and  without  this 
new  land-use  class. 

4.  Use  an  air  quality  model  to  analyze  effects  of  including  this  new  land-use 
class. 

5.  Design  a  field  study  to  measure  ammonia  emissions  from  agricultural  soils, 
before  and  after  they  have  been  sprayed  with  swine  waste  effluent. 

6.  Develop  an  observationally  based  model  relating  ammonia  emissions  to  soil 
parameters  as  well  as  present  a  fundamental  mechanistic  mass  transport 
model.  Assess  both  models  with  field  data. 
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Source 

NO,' 

Fossil  Fuel  Combustion 

22 

Biomass  Burning 

7.9 

Lightening  Discharges 

5.0 

Release  from  Soils 

7.0 

NH3  Oxidation 

0.9 

Stato  sphere 

0.64 

High-flying  aircraft 

0.85 

Total  Sources* 

44 

Sink 

NO,* 

Wet  deposition,  land 

19 

Wet  deposition,  oceans 

8 

Wet  deposition,  combined 

27 

Dry  deposition 

16 

Total  Sinks 

43 

(1  Tg=10'2g). 
Source: 

^  Lee  et  al.,  1997 
^  Logan,  1983 


Table  1.1.  Global  tropospheric  NOx  budget.  Source:  modified  from  Warneck,  2000. 
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Source 

NHj^ 

Coal  Combustion 

2 

Automobiles 

0.2 

Biomass  Burning 

5 

Domestic  animals 

32 

Wild  animals 

- 

Human  excrements 

4 

Soil/plant  emissions 

10 

Fertilizer  losses 

9 

Oceans 

13 

Total  Sources* 

75 

Sink 

NH,' 

Wet  deposition 

46 

Dry  deposition 

10 

Reaction  with  OH 

1 

Total  Sinks 

57 

*  It  is  accepted  that  the  apparent  difference  between  total  NH3  sources  and  sinks 
represents  uncertainties  in  identified  budget  terms,  not  atmospheric  accumulation. 

[-]  Indicates  insignificant  or  unavailable  terms. 

(lTg=10’'  g). 

Source: 

^  Schlesinger  and  Hartley,  1992 


Table  1.2.  Global  tropospheric  NH3  budget.  Source:  modified  from  Warneck,  2000. 
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Figure  1.1.  Atmospheric  emissions,  transport,  transformation  and  deposition  of  trace  gases. 
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Figure  1.2.  Hypothetical  growth  curve  for  an  ecosystem,  given  different  lengths  of  exposure  to  nitrogen.  Source:  Gundersen,  1992. 


Figure  1.3.  Nitric  Oxide  cycle  for  ozone  formation.  Source:  Fourth  Report  of  the 
Photochemical  Oxidants  Review  Group,  1 997. 
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Figure  1.5.  Swine  population  in  North  Carolina.  Source:  North  Caro 


Hogs  and  Pigs 
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Figure  1.6a.  Percent  of  Ammonia-Nitrogen  from  various  sources  in  the  US  for  1994.  Source:  Battye  et  al.,  1994. 
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Figure  1.6b.  Percent  of  Ammonia-Nitrogen  from  various  sources  in  North  Carolina  for  1996.  Source:  Aneja  et  al.,  1998a. 
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Figure  1.7.  Biological  nitrogen  cycle  in  the  atmosphere-soil  system.  Source:  Wameck,  2000 


In-Situ  Analysis  via 
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Figure  1.9.  Schematic  of  Dynamic  Flow-Through  Chamber. 


Figure  1.10.  Temperature  Difference  between  the  inside  and  outside  of  the  dynamic  flow-through  chamber.  This  data 
represents  the  average  of  5  experiments  measured  at  a  depth  of  approximately  5  cm. 
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Figure  1.11.  Schematic  of  modified  Thermo  Environmental  Instruments  Incorporated  (TECO)  Model  42S  chemiluminescence 
analyzer.  Instrument  is  now  capable  of  also  measuring  ammonia. 


♦  Run  #1 
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Figure  1.12.  Summary  of  5  Kaplan  experiments  conducted  during  Summer  1999  study.  Average  loss  term  (L')=0.02cm  s' 
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Abstract 

Emissions  of  nitric  oxide  (NO)  were  determined  during  late  spring  and  summer 
1995  and  the  spring  of  1996  from  four  agricultural  soils  on  which  four  different  crops 
were  grown.  These  agricultural  soils  were  located  at  four  different  sites  throughout  North 
Carolina.  Emission  rates  were  calculated  using  a  dynamic  flow-through  chamber  system 
coupled  to  a  mobile  laboratory  for  in-situ  analysis.  Average  NO  fluxes  during  late  spring 
1995  were:  50.9  ±  47.7  ng  N  m'^  s  ‘  from  soil  planted  with  com  in  the  lower  coastal  plain. 
Average  NO  fluxes  during  summer  1995  were:  6.4  ±  4.6  and  20.2  ±  19.0  ng  N  m  ^  s  * 
respectively  from  soils  planted  with  com  and  soybean  in  the  coastal  region;  4.2  ±  1 .7  ng 
N  m'^  s'*  from  soils  planted  with  tobacco  in  the  piedmont  region;  and  8.5  +  4.9  ng  N  m'^ 
s'*  from  soils  planted  with  corn  in  the  upper  piedmont  region.  Average  NO  fluxes  for 
spring  1996  were:  66.7  ±  60.7  ng  N  m'^  s'*  from  soils  planted  with  wheat  in  the  lower 
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coastal  plain;  9.5  ±  2.9  ng  N  s'*  from  soils  planted  with  wheat  in  the  coastal  plain;  2.7 
±  3.4  ng  N  m'^  s'*  from  soils  planted  with  wheat  in  the  piedmont  region;  and  56.1  ±53.7 
ng  N  m'^  s'*  from  soils  planted  with  corn  in  the  upper  piedmont  region.  An  apparent 
increase  in  NO  flux  with  soil  temperature  was  present  at  all  of  the  locations.  The 
composite  data  from  all  the  research  sites  revealed  a  general  positive  trend  of  increasing 
NO  flux  with  soil  water  eontent.  In  general,  increases  in  total  extractable  nitrogen  (TEN) 
appeared  to  be  related  to  increased  NO  emissions  within  each  site,  however  a  consistent 
trend  was  not  evident  across  all  sites. 
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Introduction 


Ozone  is  produced  in  the  troposphere  through  photochemical  processes  involving 
oxides  of  nitrogen  (NOx  “  NO  +  NO2)  and  volatile  organic  compounds  (VOCs). 
Tropospheric  ozone  is  an  important  photochemical  air  pollutant  which  increases 
respiratory  illness,  damages  crops,  and  causes  other  environmental  problems.  Currently, 
the  only  known  pathway  for  the  production  of  ozone  is  the  photolysis  of  NO2  (NO2 
NO  +  0(^P)),  which  further  reacts  with  O2  to  produce  ozone  (O3)  by  the  reaction  0(  P)  + 
O2  ->•  O3.  In  a  pseudo-photostationary  environment,  the  O3  produced  would  react  with 
the  NO  that  was  generated  via  the  photolysis  of  NO2  in  the  following  reaction:  NO  +  O3 
NO2  +  O2.  Hence,  there  is  no  net  production  of  O3.  However,  in  the  real  atmosphere, 
hydroxyl  radicals  combine  with  VOC’s  to  produce  new  radicals  which  preferentially 
react  with  NO,  allowing  a  net  O3  accumulation. 

Recent  estimates  report  that  on  a  global  scale  biogenic  emissions  of  NO  are 
comparable  to  anthropogenic  sources  (~  20  Tg  NO-N  year  ’)  (Li  et  al.,  2000).  Regions 
such  as  the  southeast  U.S.,  which  comprise  approximately  40%  of  the  US  non-attainment 
areas,  are  classified  as  NOx  limited  and  increased  emissions  of  NO  into  the  troposphere 
are  likely  to  produce  increased  O3  concentrations  (Fehsenfeld  et  al.,  1993).  Over  40%  of 
the  ozone  non-attainment  areas  of  the  U.S.  are  in  the  Southeast.  Therefore,  in  order  to 
formulate  more  successful  ozone  control  strategies  in  the  southeast  U.S.,  the  source 
strength  of  ozone  precursors  (i.e.,  NO)  and  the  parameters  which  control  their  release 
must  be  better  understood.  An  estimate  of  the  global  source  of  biogenic  NO  emissions 
conducted  by  Davidson  and  Kingerlee  (1997)  reports  a  slight  improvement  over 
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estimates  made  by  Galbally  and  Roy  (1978).  In  their  paper,  Davidson  and  Kingerlee 
report  that  further  improvements  to  their  estimate  will  most  likely  come  from  an 
increased  biogenic  NO  dataset  and  by  gaining  a  greater  understanding  of  the  processes 
which  govern  the  release  of  NO  from  the  soil.  To  these  ends,  the  objectives  of  this  study 
were  to  increase  the  NO  dataset  and  to  attempt  to  verify  previously  reported  relationships 
between  NO  and  environmental  variables  (i.e.,  soil  temperature,  soil  water  content, 
applied  N  fertilizer). 

Methods  and  Materials 
Field  Sites  and  Sampling  Scheme 

NO  concentration  measurements  were  made  from  four  different  agricultural  soils 
(on  which  corn,  soybean,  tobacco  and  wheat  were  cultivated)  during  late  spring  and 
summer  1995  and  spring  1996  (refer  to  Tables  2.1  &  2.2  for  research  sites,  dates,  soil 
parameters  and  NO  flux  values).  All  of  the  sites,  except  for  Plymouth  were  operated  by 
the  North  Carolina  State  University  in  conjunction  with  the  North  Carolina  Agricultural 
Research  Service.  All  of  the  sites  were  managed  using  practices  typical  for  their 
respective  crops  and  physiographic  locations. 

The  design  of  this  experiment  was  to  measure  soil  emissions  of  NO  at  four 
different  locations  in  North  Carolina  during  two  different  times  of  the  year.  Each  site 
was  sampled  twice,  once  in  the  late  spring/summer  and  again  the  following  year  in  the 
spring.  The  measurement  locations  at  each  site  were  essentially  the  same  as  the  previous 
season  (within  a  distance  of  5  m).  The  daily  sampling  scheme  consisted  of  measuring 
ambient  NO  concentrations  at  ground  level  before  entering  the  dynamic  flow-through 
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chamber  and  immediately  after  exiting  the  chamber.  A  Toshiba  laptop  computer  using 
Labview  software  (National  Instruments)  was  utilized  as  the  data  acquisition  system. 

The  system  produced  60  second  rolling  average  concentrations,  which  were  then  binned 
and  averaged  every  1 5  minutes.  A  daily  experiment  consisted  of  placing  the  chamber  on 
the  stainless  steel  collar  approximately  one  hour  prior  to  data  acquisition.  At  the 
conclusion  of  each  daily  experiment,  the  collar  was  relocated  to  a  random  position  within 
a  10m  radius  of  the  mobile  laboratory,  allowing  a  minimum  of  12  hours  for  any  effect  on 
the  soil,  due  to  insertion  of  the  collar,  to  dissipate.  The  randomization  of  the  collar  and 
chamber  system  placement  resulted  in  it  being  located  both  in  and  between  the  rows  of 
the  row  crops  (corn,  soybean,  tobacco),  however,  plants  were  never  inside  the  enclosure. 
Instrumentation  and  Flux  Calculation 

The  chamber  design,  associated  mass  balance  equation  and  calibration  procedures 
are  described  in  full  in  Chapter  I. 

Soil  Analysis 

Soil  samples  and  soil  analyses  are  described  in  full  in  Chapter  I. 

Results  and  Discussion 

Soil  Temperature 

Figure  2.1  is  an  hourly  averaged  diurnal  profile  of  NO  flux  versus  time  of  day  for 
each  site  and  crop  type.  This  graph  has  been  split  to  reflect  2  different  ranges  in  NO  flux 
so  that  diurnal  variations  can  be  better  discerned.  In  general,  both  graphs  reveal  that 
maximum  NO  emissions  occur  in  the  afternoon  when  soil  temperatures  are  typically  at  a 
maximum  and  reduce  to  lower  values  in  the  late  afternoon  hours  when  soil  temperatures 
are  typically  lower.  The  bottom  graph  also  reveals  that  several  sites  have  a  local 
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maximum  of  NO  emissions  in  the  morning.  This  morning  peak  has  been  observed  by 
other  researchers  (Johansson  and  Granat,  1984;  Roelle,  1996),  however  a  solid 
understanding  as  to  the  cause  of  this  phenomenon  has  yet  to  be  determined.  It  has  been 
hypothesized  that  this  morning  peak  is  a  result  of  plant  roots  which  may  exude  organic 
substances  in  the  morning  hours  which  denitrifying  bacteria  can  utilize  to  produce  NO 
(Johansson  and  Granat,  1984).  Aneja  et  al.  (1995)  reported  a  negative  correlation 
between  NO  flux  and  soil  temperature,  which  was  explained  by  a  combination  of 
moisture  and  heat  stress  on  the  soil  microorganisms  in  the  top  20  cm  of  soil.  This 
explanation  is  not  sufficient  in  this  study,  because  although  there  was  a  brief  period  when 
the  correlation  was  negative,  the  greatest  values  of  NO  flux  were  in  the  afternoon,  when 
soil  temperatures  were  at  a  maximum. 

It  is  generally  found  that  in  the  range  of  temperatures  between  288-308  °K, 
biochemical  reactions  will  rise  exponentially  with  temperature  (Warneck,  2000).  Several 
researchers  have  investigated  this  relationship,  and  in  fact  have  observed  there  to  be  an 
exponential  dependence  of  NO  flux  on  soil  temperature  (Williams  and  Fehsenfeld,  1991; 
Sullivan  et  al.,  1996;  Thornton  et  al.,  1997;  Roelle  et  al.,  1999).  The  research  conducted 
during  this  experiment  revealed  similar  results,  however,  the  exponential  dependence  of 
NO  flux  on  soil  temperature  was  not  consistent  at  all  sites. 

Plymouth,  NC,  during  the  late  spring  1995  experiment  displayed  some 
dependence  of  NO  flux  on  soil  temperature  [log  (NO  flux)  =  0.083(Soil  T)  -  0.41;  R^ 
=0.27]  for  the  entire  16  day  period.  However,  this  research  site  was  characterized  by  an 
additional  application  of  N  fertilizer  midway  through  the  experiment  allowing  the  data  to 
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be  segregated  into  periods  ‘before’  and  ‘after’  fertilization.  For  the  period  before  the 
additional  N  fertilization,  there  was  virtually  no  relationship  between  temperature  and 
flux  [Log  (NO  flux)  =  0.023(Soil  T)  +  0.89;  =0.10].  However,  for  the  period  after  the 

additional  N  fertilization  (Excluding  one  day  due  to  the  remnants  of  Hurricane  Allison 
flooding  the  field),  the  relationship  increased  [Log  (NO  flux)  =  0.23(Soil  T)  -  3.823; 
=0.58].  Additionally  there  was  some  dependence  of  NO  flux  on  soil  temperature 
observed  at  a  soybean  crop  located  in  Kinston,  NC  during  the  summer  1995  research 
period  when  one  day  was  excluded  due  to  heavy  rains  the  previous  day  [Log  (NO  flux)  = 
0.45(Soil  T)  - 1 1.33;  R^  =.35].  There  was  also  some  dependence  observed  at  the  tobacco 
crop  in  Oxford,  NC  [Log  (NO  flux)  =  0.09(Soil  T)  -  2.00;  R^  =.56].  At  the  remainder  of 
the  research  sites,  this  same  exponential  dependence  on  soil  temperature  was  difficult  to 
detect. 

Figure  2.2  is  a  composite  of  the  data  from  each  site  with  hourly  averaged  soil 
temperature  plotted  versus  hourly  averaged  NO  flux.  This  graph  and  regression  analysis 
reveals  a  significant  improvement  in  R^  values  for  data  which  was  collected  at  Plymouth 
and  Reidsville,  which  were  being  cultivated  with  corn  during  approximately  86%  of  the 
measurement  period.  These  results  are  more  in  line  with  similar  analyses  conducted  by 
Thornton  et  al.,  1997,  who  reported  R^  values  of  0.87  for  a  dataset  comprised  of 
approximately  9900  observations,  of  which  94%  were  from  soils  where  corn  was  being 
cultivated.  This  deviation  from  the  often-cited  relationship  between  NO  flux  and  soil 
temperature  at  Kinston  and  Oxford  suggests  that  factors  other  than  soil  temperature  are 
acting  to  control  the  flux  of  NO.  Researchers  in  California  conducting  extensive  research 
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on  a  very  robust  NO  data  set  report  that  depending  on  location  and  time,  other  variables 
(i.e.  crop  type,  soil  type,  water  filled  pore  space,  N  content)  can  have  stronger  influences 
than  soil  temperature  on  NO  flux  (Matson  et  ah,  1997).  Additionally,  Sullivan  et  al. 
(1996)  reported  that  differences  in  a  crops  growth  stage  and  the  age  and  amount  of  root 
biomass  further  act  to  influence  the  NO  flux/temperature  dependence. 

Total  Extractable  Nitrogen 

Another  variable,  which  has  also  been  found  to  control  the  NO  flux,  is  total 
extractable  nitrogen  (TEN)  (NH4^  +  NO3').  Figure  2.3  is  a  graph  of  the  total  average  NO 
flux  (stippled  bar)  (0900-1700)  and  TEN  on  the  secondary  axis  plotted  for  each  research 
site.  Three  sampling  days  were  removed  from  this  analysis,  which  occurred  during  or 
immediately  after  rain  events  at  the  site.  This  graph  does  reveal  a  trend  of  NO  emissions 
responding  to  increasing  and  decreasing  amounts  of  TEN,  however  a  consistent 
relationship  was  unable  to  be  detected.  Researchers  conducting  similar  studies  have 
reported  on  the  effects  of  organic  and  inorganic  nitrogen  content  in  soils  on  NO 
emissions  and  in  general  have  found  that  soils  with  higher  N  content  produce  higher  NO 
emissions  (Davidson,  1991  a,b;  Cardenas  et  ah,  1993;  Potter  et  al.,  1996;  Sullivan  et  ah, 
1996;  Roelle  et  ah,  1999).  If  the  research  conducted  in  spring  1995  at  Plymouth,  NC  is 
taken  as  an  example,  it  is  evident  that  TEN  alone  is  not  controlling  the  release  of  NO 
from  the  soil.  At  this  site,  a  small  increase  in  TEN  was  associated  with  a  much  larger 
increase  in  NO  emissions.  A  plausible  explanation  as  to  the  lack  of  a  consistent 
relationship  in  our  data  between  TEN  and  NO  emissions  could  be  the  result  of  our  soil 
sampling  method.  One  method  of  N  fertilizer  application  is  to  use  applicator  nozzles 
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which  drag  across  the  soil  surface  (2  cm  wide  strip  of  concentrated  liquid  fertilizer).  This 
applieation  method  was  used  on  the  eom  erop  at  the  Plymouth  site.  Upon  drying,  there  is 
no  way  to  discern  where  exactly  this  fertilizer  band  is  in  the  interrow  (1  m  width).  It  is 
possible  that  the  fertilized  strip  was  in  the  chamber  footprint  (diameter  27  cm)  but  our 
soil  core  sample  (diameter  8.3  cm)  missed  the  band.  Further,  it  is  possible  that  a  soil 
sample,  whieh  occurred  immediately  following  fertilization,  colleeted  the  fertilized  strip 
which  had  yet  to  penetrate  the  soil  surfaee. 

Statistieal  analysis  relating  TEN  with  NO  in  Figure  2.3  reveals  a  noisy  data  set 
with  an  R-squared  value  of  only  0.15.  However,  this  relationship  does  improve  and 
becomes  highly  significant  (p<0.0001)  when  those  flux  values  greater  than  100  ng  N  m‘^ 
s"'  are  removed  from  the  regression.  The  justifieation  for  removing  3  (from  Spring  95 
data  set)  of  the  5  data  points  greater  than  100  ng  N  m’^  s’*  is  that  the  data  were  eollected 
immediately  (within  2  weeks)  after  fertilization.  In  contrast,  at  all  the  remaining 
measurement  sites,  there  were  no  data  collected  within  at  least  one  month  of  fertilization. 
These  results  suggest  that  TEN  can  be  used  effectively  in  NO  flux  models,  although 
further  research  will  be  needed  to  determine  under  what  physiologieal  and  environmental 
conditions  this  relationship  can  be  best  applied.  Sullivan  et  al.  (1996),  during  a  field 
study  with  an  equivalent  experimental  design,  experienced  similar  results  and  reported 
that  the  differences  in  a  plants  physiological  growth  stages,  during  which  nitrogen 
demands  of  the  plants  will  vary,  confound  the  relationship  between  NO  flux  and  TEN. 
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Soil  Moisture 

The  role  of  soil  moisture  on  NO  flux  has  been  examined  by  several  researchers 
(Cardenas  et  al.,  1993;  Valente  and  Thornton,  1993;  Ormeci  et  al.,  1999).  Typically, 
each  soil  type  will  have  a  range  of  soil  moisture  which  optimizes  NO  flux.  For  example, 
Cardenas  et  al.  (1993)  reported  that  the  optimum  range  of  soil  water  content  for  sandy 
loam  soils  was  between  9-18%.  Moisture  values  above  a  soil’s  optimum  zone  will 
generally  decrease  NO  flux  due  to  pore  spaces  filling  with  water  and  inhibiting  gas 
transport.  Moisture  values  below  a  soils  optimum  zone  will  generally  lead  to  decreased 
NO  emissions  as  a  result  of  moisture  stress  to  the  soil  microbes.  Figure  2.4  is  a  graph  of 
the  daily  averaged  (0900-1700)  NO  flux  values  versus  %  soil  moisture  for  the  four 
research  sites  where  the  experiments  were  conducted.  These  results  partially  support 
Cardenas’  et  al.  (1993)  research  in  that  NO  flux  tends  to  respond  positively  to  increased 
moisture,  and  then  begins  to  decrease  again  at  moisture  levels  greater  than  20%. 

However,  it  is  also  evident  that  with  R-squared  values  of  0.2,  this  variable  alone  cannot 
adequately  predict  NO  emissions  from  biogenic  processes.  Using  the  same  rationale  as  in 
Figure  3,  the  data  points  greater  than  100  ng  N  m‘^  s"'  were  removed  from  the  data  set 
and  then  reanalyzed.  However  unlike  Figure  2.3,  this  approach  only  resulted  in  a  slightly 
better  model  prediction  (R  =.3). 

NO  Response  to  Soil  Parameters 

Chameides  et  al.  (1988  and  1994)  have  reported  on  the  importance  of 
including  both  anthropogenic  and  biogenic  NO  emissions  when  attempting  to  predict  the 
concentrations  of  tropospheric  ozone.  One  of  the  methods  used  in  air  quality  models  for 
estimating  soil  NO  emissions  is  the  Biogenic  Emissions  Inventory  System  (BEIS2) 
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model.  This  model  utilizes  soil  temperature  data  and  an  emission  factor  which  is  based 

on  crop  type,  the  type  and  amount  of  fertilizer  applied,  and  other  chemicals  applied  to  the 

soil  (Birth  et  al.,  1995).  The  NO  flux  is  then  calculated  by; 

NO  Flux  (ng N  s’')  =  A  *  EXP(0.71  *Ts) 

where  A  =  experimentally-derived  coefficient  which  is 

dependent  on  the  land  use  category 

Ts  =  soil  temperature 

(Williams  et  al.,  1992). 

This  temperature  dependence  model  has  been  questioned  by  several  researchers  who 
have  found  that  temperature  alone  does  not  adequately  explain  the  flux  of  NO  in  their 
measurements  (Matson  et  al.,  1997;  Roelle  et  al.,  1999;  Li  et  al.,  1999).  In  fact,  Matson 
et  al.  (1997)  reported  that  for  most  crops  in  their  research,  percent  water-filled  pore  space 
(WFPS)  was  equally  as  important  as  temperature  in  predicting  NO  emissions. 

Table  2.3  lists  Williams’  values  for  A  as  well  as  those  calculated  for  this  study. 
The  range  of  soil  temperatures  (Ts)  with  the  corresponding  fluxes  produced  using  the 
BEIS2  model  versus  the  actual  data  measured  during  this  study  is  also  on  this  table. 

There  was  reasonable  agreement  between  A  factors  for  the  com  crop,  where  Williams’  A 
value  is  9  ng  N  m'^  s''  and  the  A  value  from  this  study  is  7  ng  N  m'^  s’*.  However,  the 
remaining  crops,  (soybean,  tobacco,  and  wheat)  varied  from  Williams’  model  by  factors 
of  15,  6.7,  and  3.7  respectively.  Research  conducted  in  the  San  Joaquin  Valley  in 
California  produced  even  larger  deviations  from  Williams’  model,  where  these 
researchers  calculated  A  values  for  corn  between  0.0001  and  0.0005  ng  N  m'^  s’'  (Matson 
et  al.,  1997).  As  the  table  indicates,  using  the  BEIS2  model  to  estimate  NO  emissions 
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can  produce  values  different  by  up  to  an  order  of  magnitude  than  the  NO  emissions  that 
were  calculated  during  this  study.  < 

Conclusions  and  Recommendations 

NO  emissions  and  soil  properties  were  studied  from  several  croplands  in  order  to 
gain  a  better  understanding  of  the  chemieal  and  physical  properties  of  soil  which 
influence  NO  flux,  and  to  provide  much  needed  data  to  the  biogenie  NO  emission  dataset. 
Utilizing  a  dynamic  flow-through  chamber,  the  flux  of  NO  from  soil  was  determined  for 
four  different  regions  and  crop  types  during  both  spring  and  summer.  NO  flux  was 
generally  found  to  follow  a  diurnal  profile  with  maximum  emissions  eoinciding  with 
maximum  soil  temperatures.  Additionally,  several  sites  displayed  a  morning  peak  in  NO 
emissions,  which  is  occasionally  observed  by  some  researchers  but  has  yet  to  be 
explained.  The  exponential  dependence  of  NO  flux  on  soil  temperature  existed  at  all 
sites,  but  to  different  levels  of  significance.  It  was  also  observed  that  NO  flux  did 
respond  to  varying  amounts  of  both  total  extractable  nitrogen  (TEN)  in  the  soil  and  soil 
moisture  eontent. 

Although  relationships  between  soil  parameters  (soil  moisture,  TEN,  soil 
temperature)  were  evident,  no  one  variable  or  combination  of  variables  has  yet  been 
found  to  adequately  model  the  flux  of  NO  from  agricultural  soils.  Parkin  (1993) 
addressed  the  difficulty  in  modeling  a  system  with  such  high  variability,  yet  stressed  the 
importance  of  continued  research,  as  improved  estimates  will  only  be  achieved  as  we 
gain  a  greater  understanding  of  the  processes  causing  the  variability.  The  data  from  this 
observationally  based  study  raise  concerns  about  the  current  practice  of  basing  emission 
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estimates  solely  on  temperature  and  land  cover  type.  Errors  in  soil  emission  estimates 
caused  by  ignoring  the  influence  of  parameters  such  as  soil  moisture  and  TEN,  which 
was  pointed  out  in  the  comparison  of  BEIS2  model  estimates  to  actual  calculated  fluxes, 
may  also  hinder  the  ability  of  ozone  models  to  simulate  VOC/NOx  emission  control 
scenarios.  With  the  addition  of  this  data  set  reported  here,  and  with  data  such  as  reported 
by  Thornton  et  al.  (1997),  the  scientific  community  should  combine  other  observational 
data  to  build  a  model  which  includes  a  broader  set  of  environmental  parameters  (i.e., 
%WFPS,  TEN,  soil  type),  which  will  likely  lead  to  better  estimates  of  NO  emissions 
from  agricultural  soils. 
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Location 

Crop 

Soil  Temp 
CC) 

TEN 

mg  N 

%  Soil 
Moisture 

NO 

(ng  N  m'^  s*b 

(Kg  dry  Soil) 

Plymouth 

Corn 

Avg 

24.1 

51.0 

21.1 

50.9 

N=16 

Std  Dev 

3.2 

26.2 

2.8 

47.7 

n=759 

Min 

16.4 

24.0 

17.7 

4.2 

Max 

32.7 

116.0 

27.7 

264.7 

Kinston 

Com 

Avg 

25.6 

9.8 

12.4 

6.4 

N=5 

Std  Dev 

2.6 

1.7 

1.4 

4.6 

n=203 

Min 

21.1 

8.0 

10.7 

2.1 

Max 

32.6 

12.1 

14.0 

37.2 

Kinston 

Soybean 

Avg 

25.8 

14.2 

12.8 

20.2 

N=4 

Std  Dev 

3.2 

3.4 

0.8 

19.0 

n=276 

Min 

21.5 

11.0 

11.6 

1.7 

Max 

31.9 

19.0 

ISA 

96.8 

Oxford 

Tobacco 

Avg 

27.4 

8.0 

5.6 

4.2 

N=6  " 

Std  Dev 

2.2 

2.5 

2.3 

1.7 

n=285 

Min 

23.5 

6.0 

2.7 

1.0 

Max 

32.5 

13.0 

8.1 

13.0 

Reidsville 

Com 

Avg 

23.0 

12.8 

11.3 

8.5 

N-4 

Std  Dev 

2.5 

'  11.9 

2.4 

4.9 

n=289 

Min 

19.7 

4.0 

10.0 

1.4 

Max 

29.0 

32.0 

15.6 

20.5 

Plymouth 

Wheat 

Avg 

14.3 

9.4 

22.7 

66.7 

N=4 

Std  Dev 

3.2 

4.4 

1.2 

60.7 

n=229 

Min 

5.5 

5.3 

21.3 

2.7 

Max 

19.2 

13.9 

23.8 

175.6 

Kinston 

Wheat 

Avg 

17.1 

3.2 

9.4 

9.5 

N=4 

Std  Dev 

2.9 

0.8 

2.5 

2.9 

n=166 

Min 

9.5 

2.2 

6.3 

4.3 

Max 

21.5 

4.0 

11.7 

19.9 

Oxford 

Wheat 

Avg 

15.4 

1.8 

6.9 

2.7 

N=4 

Std  Dev 

2.4 

0.5 

1.3 

3.4 

n=187 

Min 

10.7  ! 

1.3 

5.8 

0.0 

Max 

19.3  i 

2.4 

8.7 

25.5 

Reidsville 

Corn 

Avg 

19.1 

19.7 

21.7 

56.1 

N=4 

Std  Dev 

4.9 

13.0 

2.4 

53.7 

n=161 

Min 

11.0 

8.6 

19.8 

4.5 

Max 

28.4 

36.5 

24.9 

191.9 

Table  2.2.  Soil  parameters  and  NO  flux  for  various  agricultural  soils  on  which  different 
crops  were  grown.  N=number  of  sampling  days;  n=number  of  NO  concentration 
measurements.  TEN=total  extractable  nitrogen 
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Corn 

Soybean 

Tobacco 

Wheat 

Emission  Factor  (A)  [this  study] 
(ng  N  m'^  s'*) 

7.0 

3.0 

0.6 

11.0 

Williams’  Emission  Factor  (A) 
(ng  N  m'^  s'*) 

9.0 

0.2 

4,0 

3.0 

Tsoii  Range  (  C) 

11.0-32.7 

21.5-31.9 

23.5-32.5 

5.5-21.5 

NO  Flux  [using  BEIS2] 

(ng  N  m'^  s'*) 

20.0-92.0 

1.0-2.1 

21.0-40.0 

4.0-14.0 

NO  Flux  [this  study] 

(ng  N  m'^  s'*) 

Range: 

2.0-265.0 

Average: 

30.0 

Range: 

2.0-97.0 

Average: 

30.0 

Range: 

1.0-13.0 

Average: 

4.0 

Range: 

0.0-176.0 

Average: 

26.0 

Table  2.3.  NO  emissions  for  this  study  and  the  results  using  the  BEIS2  model  for  the 
corresponding  range  of  soil  temperatures  observed  during  the  study. 
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Figure  2.1.  Hourly  averaged  NO  flux  (0900-1700)  for  each  measurement  site  and  crop  type.  Note  the  change  in  scale  on  the 
vertical  axes-graphs  have  been  split  so  that  diurnal  trends  would  be  easier  to  discern.  The  1  &  2  attached  to  each  site  name 
indicates  the  season  the  site  was  sampled,  l=late  spring/summer  95;  2=spring  96. 


♦  Reidsville 
■  Plymouth 
▲  Kinston 
50.0  -  »  Oxford 


Figure  2.2.  Hourly  averaged  NO  flux  plotted  versus  h( 
lata  from  the  four  sites  where  the  data  was  collected. 
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Figure  2.3.  NO  flux  between  0900-1700  (stippled  bars)  and  total  extractable  N  (on  secondary  axis)  plotted  for  each  site.  The  1  & 
2  attached  to  each  site  name  indicates  the  season  the  site  was  sampled,  l=late  spring/summer  95;  2=spring  96.  Error  bars  represent 
one  standard  deviation.  One  sampling  day  at  Kinston  1 -Soybean,  Plymouth^- Wheat,  and  Reidsville2-Com,  were  not  included  in 
this  dataset  due  to  rainshowers  occurring  at  the  site. 


■  Plymouth  | 


LO 

CO 


s  ^.UI  N  Bu)  xnid  ON 


64 


%  Soil  Moisture 

Figure  2.4.  Average  NO  flux  (0900-1700)  vs  %  Soil  Moisture.  Percent  soil  moisture  was  determined  from  a  soil  sample 
taken  from  the  center  of  the  chamber  footprint  at  the  end  of  each  measurement  period.  Data  have  been  segregated  into  sites. 


Chapter  III.  Nitric  Oxide  Emissions  from  Soils  Amended  with 
Municipal- Waste  Biosolids 
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Department  of  Marine,  Earth,  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Published  in:  Atmospheric  Environment,  in  press,  2001. 

Abstract 

Spreading  of  nitrogen-rich  municipal  waste  biosolids  (NOj'-N  <  256  mg-N/kg  dry 
weight,  NHj-N  -23,080  mg-N/kg  dry  weight.  Total  Kjeldahl  N  ~  41,700  mg-N/kg  dry 
weight)  to  human  food  and  non-food  chain  land  is  a  practice  followed  throughout  the 
US.  This  practice  may  lead  to  the  recovery  and  utilization  of  the  nitrogen  by  vegetation, 
but  it  may  also  lead  to  emissions  of  biogenic  nitric  oxide  (NO),  which  may  enhance 
ozone  pollution  in  the  lower  levels  of  the  troposphere.  Recent  global  estimates  of 
biogenic  NO  emissions  from  soils  are  cited  in  the  literature,  which  are  based  on  field 
measurements  of  NO  emissions  from  various  agricultural  and  non-agricultural  fields. 
However,  measurements  of  biogenic  emissions  of  NO  from  soils  amended  with  biosolids 
are  lacking.  Utilizing  a  state-of-the-art  mobile  laboratory  and  a  dynamic  flow-through 
chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were  measured  during  the 
Summer/Fall  of  1999  and  Winter/Spring  of  2000  from  an  agricultural  soil  which  is 
routinely  amended  with  municipal  waste  biosolids.  The  average  NO  flux  for  the 
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Summer/Fall  time  period  (9  June  1999-  3  September  1999)  was  57.8  ±  34.6  ng  NO-N  m'^ 
s''.  Biosolids  were  applied  during  September  1999  and  the  field  site  was  sampled  again 
during  Winter/Spring  2000  (28  February  -  9  March  2000),  during  which  the  average  flux 
was  3.6  ±  3.9  ng  NO-N  m'^  s'.  Summer  accoimted  for  60%  of  the  yearly  emission  while 
fall,  winter  and  spring  accounted  for  20%,  4%  and  16%  respectively.  Field  experiments 
were  conducted  which  indicated  that  the  application  of  biosolids  increases  the  emissions 
of  NO  and  current  techniques  to  estimate  biogenic  NO  emissions  would  have 
underestimated  the  NO  flux  from  this  field  by  a  factor  of  26.  Soil  temperature  and  % 
Water  Filled  Pore  Space  (%WFPS)  were  observed  to  be  significant  variables  for 
predicting  NO  emissions,  however  %WFPS  was  found  to  be  most  significant  during  high 
soil  temperature  conditions.  In  the  range  of  pH  values  found  at  this  site  (5.8  ±  0.3),  pH 
was  not  observed  to  be  a  significant  parameter  in  predicting  NO  emissions. 
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Introduction 


Although  there  have  been  many  experiments  conducted  that  have  measured  NO 
emissions  from  various  soil  types  (Johansson  and  Granat,  1984;  Kaplan  et  al.,  1988; 
Aneja  et  al.,  1995;  Williams  and  Fehsenfeld,  1991;  Kim  et  al.,  1994,  Roelle  et  al.,  2001), 
relatively  few  have  included  intensively  managed  agricultural  soils,  or  continued 
measurements  for  substantial  periods  of  time  (Anderson  and  Levine,  1987;  Sullivan  et  al., 
1996).  In  the  Southeast  U.S.,  which  is  NOx  limited,  an  increase  in  NOx  emissions  is 
believed  to  produce  a  corresponding  increase  in  O3  levels  (Fehsenfeld  et  al.,  1993). 

Ozone  adversely  affects  human  health  and  has  also  been  shown  to  reduce  crop  yield. 
Studies  show  that  prolonged  exposure  to  high  ozone  levels  causes  persistent  functional 
changes  in  the  gas  exchange  region  of  the  lungs.  Additionally,  ozone  plays  a  critical  role 
in  controlling  the  chemical  lifetimes  and  the  reaction  products  of  many  atmospheric 
species  (National  Research  Council,  1991).  Gaseous  nitric  acid  (FINO3),  the  end  product 
of  NO  reactions  in  the  atmosphere,  combines  with  either  aerosols  or  water  in  the 
atmosphere,  and  is  removed  via  rain,  snow,  or  other  deposition  processes,  as  acidic 
deposition. 

Davidson  and  Kingerlee  (1997)  reviewed  over  60  papers  in  order  to  develop  a 
global  inventory  of  NO  emissions  from  soils.  Their  revised  estimate  included  emissions 
from  grasslands,  desert,  marshlands,  tundra,  forests,  fertilized  and  unfertilized  croplands. 
However  none  of  the  papers  addressed  NO  emissions  from  land  treated  with  biosolids 
from  wastewater  treatment  facilities.  A  sample  of  the  112  different  NO  flux  values 
reported  by  Davidson  and  Kingerlee  (1997)  is  presented  in  Table  3.1.  Although  eight  of 
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the  sixty  papers  had  larger  emissions  than  this  field  study,  most  of  the  sites  reported 
lower  NO  emissions.  In  fact,  the  overall  average  of  1 12  reported  NO  flux  averages  was 
17.4  ng  N-NO  yf',  while  the  average  from  this  field  study  was  57.8  ng  N-NO  yf'.  In  the 
U.S.  alone,  more  than  6  million  dry  metric  tons  of  biosolids  are  generated  annually,  and 
approximately  30%  is  spread  or  injected  into  the  land  (Peirce  and  Aneja,  2000).  Given 
that  there  is  a  lack  of  field  data  from  soils  amended  with  biosolids  may  indicate  a  possible 
error  in  the  global  estimates  of  NO. 

Methods  and  Materials 
Biosolids 

The  biosolids  applied  to  the  field  in  this  study  were  from  a  wastewater  treatment 
facility  located  in  North  Carolina.  The  population  which  this  facility  serves  is 
approximately  sixty  thousand  people.  The  flow  of  wastewater  through  the  plant  is 
approximately  30.3  x  10®  liters  day  '.  The  wastewater  treatment  facility  consists  of 
preliminary  treatment,  primary  treatment,  secondary  treatment,  disinfection  and  sludge 
treatment  and  disposal  (Figure  3.1).  Preliminary  treatment  is  the  part  of  the  process 
where  the  wastewater  is  first  screened  for  large  objects  via  mechanical  filters  and  then 
enters  grit  chambers  where  sand  and  inorganic  solids  settle  out.  During  primary 
treatment,  the  heavier  solids,  which  have  settled  to  the  bottom,  and  the  lighter  solids 
which  have  risen  to  the  top,  are  removed  from  the  tanks  by  skimming/scraping  processes. 
The  wastewater  emitted  from  the  primary  treatment  then  enters  the  secondary  treatment. 
During  this  step,  the  liquid  is  aerated  and  microbial  biomass  is  introduced  to  induce 
nitrification.  At  the  conclusion  of  the  secondary  treatment,  the  liquid  is  treated  with 
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chlorine  to  remove  any  infectious  organisms  and  then  the  water  is  released  to  a  local 
stream.  The  heavier  and  lighter  solids  combined  together,  called  biosolids  or  sludge,  are 
then  sent  to  the  treatment  part  of  the  facility  where  it  is  prepared  for  land  application.  At 
this  stage  the  biosolids  are  then  transported,  via  tanker  trucks,  to  facilities  where  it  is 
stored  until  it  can  be  applied  to  various  fields.  The  biosolids  that  were  used  for  simulated 
fertilization  experiments  were  collected  from  the  purging  faucet  of  a  digested  sludge 
storage  tank  outlet.  The  biosolids  were  discarded  within  two  weeks  of  being  collected  to 
avoid  possible  changes  in  their  biological  and  chemical  composition  (Droste  1997). 

Table  3.2  is  a  description  of  the  chemical  composition  of  the  biosolids.  The  biosolids 
were  also  analyzed  for  metals  with  the  results  on  a  dry  weight  basis  as  follows:  arsenic 
(2.8  mg/kg),  cadmium  (1.8  mg/kg),  chromium  (29.5  mg/kg),  copper  (342  mg/kg),  lead 
(48.5  mg/kg),  mercury  (2.9  mg/kg),  molybdenum  (5.7  mg/kg),  nickel  (8.4  mg/kg), 
selenium  (2.0  mg/kg),  and  zinc  (509  mg/kg).  No  biosolids  were  applied  to  the  crop 
during  the  in-situ  measurement  period,  therefore  fertilization  was  simulated  by  applying 
approximately  400  ml  of  biosolids  (an  amount  equivalent  to  what  the  dynamic  chambers 
footprint  area  would  typically  receive)  onto  the  soil  surface  30  minutes  prior  to  any 
measurements.  Field  measurements  were  made  during  summer  1999  (June- August), 
winter  2000  (February-March)  and  spring  2000  (May- June).  The  field  studies  on  the 
biosolid  amended  soils  were  conducted  during  the  spring  and  summer  measurement 
campaigns  to  represent  the  actual  times  when  the  crop  would  receive  biosolid 
applications. 
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Instrumentation  and  Flux  Calculation 


The  chamber  design,  associated  mass  balance  equation  and  calibration  procedures 
are  described  in  full  in  Chapter  I. 

Physiographic  Location 

NO  concentration  measurements  were  made  on  a  property  which  is  owned  by  a 
municipal  water  and  sewer  authority  located  in  central  North  Carolina.  Through  a 
cooperative  agreement,  a  local  farmer  cultivates  the  field  and  the  water/sewer  authority 
applies  nutrients  in  the  form  of  biosolids  and  limes  as  necessary.  The  crop  grown  during 
our  measurement  campaign  was  a  small  grain  including  Kenland  Red  Clover  and 
Hallmark  Orchard  Grass,  which  was  planted  in  fall  1 998  and  amended  with  biosolids  in 
January  1999  with  approximately  80  kg  N/ha.  (Note:  this  year’s  fertilizer  data  was  not 
available  and  these  values  are  from  fertilizer  records  of  the  preceding  year).  The 
measurements  at  this  field  site  occurred  9  June  —  3  September  1999  and  28  February-9 
March  2000).  The  field  site  was  mowed  on  July  27,  limed  at  the  end  of  August  and  then 
received  another  application  of  biosolids  in  mid  September  (~80  kg  N/ha). 

Sampling  Scheme 

For  a  complete  description  of  the  sampling  method  used  see  Chapter  I  (Sampling 
Scheme). 

Soil  Analysis 

The  mineral  soils  at  this  research  site  are  classified  in  the  soil  survey  as  Lignum 
silty  loam,  which  is  described  as  a  moderately  well  drained,  slowly  permeable  soil  on  the 
uplands  with  slopes  of  0  to  3%  (Dunn,  1977).  However,  soil  surveys  generally  cover 
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broad  areas,  and  a  soil  texture  determination  of  our  specific  soil  cores  within  the  larger 
research  site  identified  our  samples  soils  as  a  sandy  loam  (Tabachow,  2000).  Soil 
temperature  was  recorded  every  minute,  and  these  values  were  binned  and  averaged  every 
15  minutes  using  a  Campbell  Scientific  soil  temperature  probe  (accuracy  ±  3%)  inserted  5 
cm  into  the  soil,  adjacent  to  the  chamber. 

Results  and  Discussion 

Effects  of  Soil  Nitrogen  Content  and  pH 

The  data  set  plotted  in  Figures  3.2  (a-c)  represent  daily  average  NO  flux  values 
(0900-1700)  versus  the  daily  sampled  nitrogen  content  as  Total  Kjeldahl  Nitrogen,  NH3 
and  NOj'  respectively.  As  indicated  by  the  data  in  Figure  3.2a,  TKN  is  not  found  to  be  a 
significant  variable  in  predicting  NO  emissions.  However,  Figures  3.2b  and  3.2c  reveal 
that  the  addition  of  ammonia  or  nitrate  to  the  soil  produces  a  corresponding  increase  in 
NO  emissions.  These  results  are  to  be  expected  based  on  the  predominant  NO  production 
pathways,  namely  nitrification  and  denitrification  processes,  and  the  fact  that  the  organic- 
N  must  first  be  converted  to  inorganic  forms  prior  to  nitrification  or  denitrification 
processes  to  occur  (Warneck,  2000;  Troeh  and  Thompson,  1993).  The  fact  that  NO 
increases  in  response  to  both  NH3  and  N03'  indicates  that  both  nitrifiers  and  denitirifiers 
are  present  in  the  soil.  However,  the  steeper  slope  in  Figure  3.2c  may  indicate  that 
denitrification  is  a  slightly  more  significant  process  in  this  biosolid  amended  soil  system. 

This  field  site  is  typically  amended  with  biosolids  twice  per  year  and  the  last 
application  prior  to  the  summer  1999  measurement  campaign  was  in  January  1999.  In 
order  to  examine  how  biosolid  applications  immediately  affect  NO  emissions,  three 
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experiments  were  conducted  in  which  one  served  as  a  control  and  two  experiments  had 
biosolids  applied  to  the  soil.  Previous  experiments,  both  by  the  North  Carolina  State 
University  Air  Quality  Research  Group  and  other  researchers  have  shown  that  adjacent 
plots  with  seemingly  homogeneous  soils  can  differ  in  NO  emissions  by  more  that  an 
order  of  magnitude  (Williams  and  Fehsenfeld,  1991;  Valente  et  ah,  1995).  In  order  to 
make  the  experimental  plot  as  homogeneous  as  possible,  a  depression  was  dug  (13”  x 
38”)  large  enough  to  accommodate  two  chambers  side-by-side,  and  deep  enough  to 
simulate  the  depth  to  which  the  soil  is  tilled.  Each  time  an  experiment  was  conducted  on 
this  plot,  the  soil  was  removed,  thoroughly  mixed  and  then  returned  where  it  was  left 
undisturbed  for  48  hours  prior  to  any  measurements.  On  both  biosolid-amended 
experiments,  the  biosolids  were  collected  directly  from  the  purging  spigot  of  the  biosolid 
storage  tank  outlet  which  fills  the  tankers  used  to  transport  the  biosolids  to  the  field  sites. 
The  biosolids  were  applied  to  the  experimental  plots  by  pouring  the  biosolids  (~80  kg 
N/ha  -  an  amount  equivalent  to  what  this  size  plot  would  typically  receive)  onto  the  soil 
surface  30  minutes  prior  to  any  measurements. 

Figure  3.3a  represents  the  first  experiment  in  which  no  biosolids  were  applied  to 
either  plot.  The  purpose  of  this  experiment  was  to  see  if  mixing  the  soil  had,  in  fact, 
produced  a  fairly  homogeneous  environment.  Although  the  graph  does  show  a  slight 
divergence  between  12:00  -  2:00  pm,  in  general  the  NO  emissions  measured  in  the  two 
plots  shows  similar  trends.  For  example,  relative  peaks  appeared  at  both  plots  at 
approximately  the  same  time  (12:00,  2:30,  4:15).  Further  similarities  were  that  emissions 
from  both  plots  decreased  during  the  first  2  hours  of  measurements  which  is  surprising 
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because  soil  temperatures  were  actually  increasing  during  this  time  period.  Given  the  fact 
that  adjacent  soils  can  differ  by  orders  of  magnitude,  and  that  the  largest  difference 
between  these  two  plots  never  exceeded  45%,  we  believed  that  we  had  created  an 
environment  similar  enough  to  examine  differences  between  two  side-by-side  chambers. 

Figure  3.3b  represents  the  first  of  2  experiments  conducted  on  a  plot  amended 
with  biosolids.  The  experiment  started  similarly  to  the  unamended  experiment  (Figure 
3.3a)  in  that  both  plots  began  with  fluxes  of  NO  that  were  within  45%  difference  of  each 
other.  However,  unlike  Figure  3.3a,  the  emission  from  the  amended  plot  immediately 
began  to  increase  while  that  in  the  unamended  plot  maintained  a  fairly  steady  value. 

From  the  results  of  this  one  experiment,  it  appears  as  if  the  plot  amended  with  biosolids 
did  affect  the  NO  emissions.  Figure  3.3c  is  a  graph  of  the  second  experiment  in  which 
one  of  the  plots  was  amended  with  biosolids.  This  experiment  began  with  NO  emissions 
at  almost  the  same  exact  NO  flux  value.  Similarly  to  Figure  3.3a,  both  plots  follow  very 
similar  diurnal  profiles,  with  relative  minimum  and  maximum  emissions  occurring  at 
approximately  the  same  time  throughout  the  profile.  The  amended  plot,  however, 
appears  to  respond  to  the  biosolids  application  with  a  greater  amplitude  in  emissions  than 
the  non-amended  plot,  especially  in  the  late  afternoon  when  the  biosolids  have  had  a 
chance  to  permeate  into  the  top  layer  of  the  soil  and  come  into  contact  with  more  of  the 
bacteria  which  are  presumed  to  be  responsible  for  the  production  of  NO. 

Relationships  between  soil  pH  and  NO  flux  in  field  studies  have  been  investigated 
in  the  past  and  have  proven  difficult  to  discern  due  to  the  fact  that  pH  values  remain  in  a 
fairly  tight  range  in  intensively  managed  agricultural  fields  (Roelle  et  ah,  1999).  Average 
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pH  values  for  this  measurement  period  were  5.8  ±  0.30.  Statistical  analysis  of  the  data 
measured  during  this  campaign  revealed  no  apparent  relationship  and  that  pH  was  not  a 
predictor  of  NO  flux  at  this  field  site.  These  results  are  corroborated  by  the  findings  of 
both  the  Ormeci  et  al.  (1999)  and  Matson  et  al.  (1997)  for  similar  ranges  of  pH.  Ormeci 
et  al.  (1999)  reported  that  NO  emissions  reach  maximum  values  for  near  neutral  to  lower 
soil  pH  values.  Under  more  acidic  conditions,  the  chemical  decomposition  of  NOj  and 
HONO,  otherwise  known  as  chemodenitrification,  has  been  suggested  to  be  responsible 
for  the  increased  availability  of  NO  in  the  soil  (Galbally,  1989).  Controlled  laboratory 
experiments,  in  which  pH  was  varied  while  maintaining  constant  soil  temperature  and 
soil  moisture  resulted  in  highest  NO  emissions  when  soil  pH  was  at  its  lowest  value  (4.3) 
(Ormeci  et  al.,  1999).  In  an  extensive  field  study  conducted  in  the  San  Juaquin  Valley, 
CA,  Matson  et  al.  (1997)  found  the  highest  NO  emissions  occurring  for  pH<  4.0.  Both 
studies  found  NO  emissions  to  be  unaffected  by  changes  in  pH  in  the  range  between  5 
and  8. 

Effects  of  Soil  Temperature  and  Soil  Water  Content 

An  important  parameter  controlling  NO  emissions,  which  has  been  studied  by 
many  investigators,  is  soil  temperature  (Williams  et  al.,  1992;  Kim  et  al.,  1994;  Sullivan 
et  al.,  1996;  Roelle  et  al.,  1999).  Biochemical  rates  in  the  soil  have  been  shown  to  rise 
exponentially  with  soil  temperature  in  the  range  between  15-35  °C  (Warneck,  2000).  The 
relationship  most  cited  in  the  literature,  and  which  is  also  used  in  the  Environmental 
Protection  Agencies’  (EPA)  Biogenic  Emissions  Inventory  System  2  (BEIS2)  is  the 
model  proposed  by  Williams  et  al.  (1992)  which  is  given  by: 
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Flux  (ng  N  m'^  s"')  =  A*Exp(0.071*T5) 


(1) 


where  A  =  experimentally  derived  coefficient  associated  with  land  use 

categories 

Ts  =  soil  temperature  ('’C) 

In  order  to  examine  similar  exponential  relationships  with  this  dataset,  the 
protocol  used  by  Thornton  et  al.  (1997)  was  adopted.  Soil  temperatures  were  segregated 
into  1.5  °C  spans  and  the  corresponding  fluxes  were  averaged  to  produce  one  mean  NO 
flux  for  each  temperature  span.  The  results  of  this  procedure  are  plotted  in  Figure  3.4, 
with  the  corresponding  regression  equation  and  values.  Because  of  the  strong 
dependence  of  NO  flux  on  soil  temperature,  the  highest  emissions  are  found  during  the 
summer.  Utilizing  the  regression  equation  in  Figure  3.4,  and  the  air-to-soil  temperature 
conversion  factors  described  in  Williams  et  al.  (1992),  we  estimate  that  Summer  (June- 
August)  accounts  for  60%  of  the  yearly  emissions,  while  Fall  (September-November), 
Winter  (December-February),  and  Spring  (March-May)  account  for  20%,  4%,  and  16% 
respectively.  The  use  of  the  BEIS2  model  would  yield  a  different  distribution  of  NO 
emissions  of  44%,  24%,  10%  and  22%  for  Summer,  Fall,  Winter  and  Spring  respectively. 
The  BEIS2  could  underestimate  the  yearly  NO  emissions  from  biosolid-amended  soils 
(Figure  4.4)  by  a  factor  of  26. 

The  relationship  between  emissions  of  NO  and  soil  temperature  discussed 
previously,  is  further  confounded  by  the  dependence  of  NO  emissions  on  soil  moisture 
content.  Researchers  have  found  that  nitrification  is  optimized  for  moisture  contents 
between  30  and  65%  WFPS  (Davidson  and  Swank,  1986;  Parton  et  al.,  1988).  Recent 
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laboratory  experiments  conducted  by  Ormeci  et  al.,  (1999)  found  that  the  range  for 
maximum  NO  emissions  occurred  between  20-45%. 

In  Figure  3.5,  the  data  are  binned  into  values  of  %  WFPS  of  25,  35  and  45  and 
soil  temperatures  of  1 5,  25,  35  "C.  NO  flux  follows  a  general  trend  of  increasing 
emissions  as  %  WFPS  increases  within  a  given  soil  temperature,  except  at  25  °C. 
Likewise,  within  a  given  %  WFPS,  NO  emissions  increase  as  soil  temperatures  increase 
except  at  35%  WFPS.  Soil  temperature  and  %WFPS  were  both  found  to  be  statistically 
significant  parameters  for  predicting  NO  emissions,  however  %WFPS  was  only 
significant  during  high  (T=35  "C;  R^=0.46)  and  low  (T=15  °C;  R^=.21)  soil  temperatures. 
The  highest  average  NO  flux  occurred  at  45%  WFPS  and  35  "C  soil  temperature  (1 19.3 
ng  N  m'^  s’’),  which  is  expected  when  considering  the  temperature  dependence  alone. 
However,  relatively  equivalent  NO  flux  values  occurring  at  35%  WFPS  and  25  “C  (1 14.3 
ng  N  m'^  s"')  seems  to  indicate  that  these  conditions  may  also  be  optimum  for  maximizing 
NO  production  from  this  biosolid  amended  field  site. 

Conclusions  and  Recommendations 

NO  emissions  from  fields  amended  with  biosolids  were  studied  during 
Summer/Fall  1 999  and  Winter/Spring  2000,  in  order  to  examine  what  environmental 
parameters  might  control  these  emissions  and  to  determine  the  impact  of  biosolids 
applications  on  soil  NO  emissions.  Soil  temperature  and  %WFPS  are  both  found  to  be 
significant  parameters  for  predicting  NO  emissions,  however  %WFPS  is  only  significant 
during  high  (T=35  "C;  R^=0.46)  and  low  (T=15  “C;  R^=.21)  soil  temperatures.  When  NO 
was  modeled  versus  Total  Kjeldahl  Nitrogen  and  the  inorganic  components,  NO 
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emissions  were  found  to  be  dependent  only  on  the  inorganic  N  species.  Further,  at  this 
field  site,  denitrification  appeared  to  be  the  dominant  process  for  NO  production.  The 
average  NO  emissions  from  this  small  grain  field  were  up  to  a  factor  of  26  higher  than 
what  the  EPA’s  currently  used  biogenic  emissions  model  would  predict  for  similar  crop 
types.  This  leads  to  the  conclusion  that  NO  emissions  are  currently  being  underestimated 
for  these  biosolid-amended  fields.  Further,  the  use  of  the  BEIS2  model  estimated  NO 
emissions  in  photochemical  models  may  also  be  erroneous. 

The  practice  of  land  applying  biosolids  not  only  serves  as  an  economical  means  of 
disposal  but  also  saves  farmers  the  expense  of  purchasing  chemically-derived  fertilizers. 
In  North  Carolina,  biosolid  amended  soils  are  applied  in  localized  areas  and  represent  a 
small  fraction  of  total  crop  land  soils  (<1%),  and  therefore  will  likely  represent  only  a 
small  fraction  of  the  total  biogenic  NO  budget.  However,  these  biosolid  amended  soils 
may  act  as  significant  sources  of  localized  O3  production,  especially  during  the  hot  and 
stagnant  periods  of  the  summer  when  biogenic  NO  emissions  and  photochemical  activity 
are  at  a  maximum.  Therefore  future  work  should  consist  of  a  modeling  study  to  examine 
the  localized  effects  biosolid  amended  soils  have  on  O3  production. 
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Cultivated  Land  -  Temperate 

Location 

Description 

Mean  Flux 
(ng  N  m'^  s'^) 

Source 

TN,  USA 

Corn  -  unfertilized 

3.06 

Thornton  et  al.,  1996 

TN,  USA 

Com  -  fertilized,  1 1 1  kg 

N/ha 

44.4 

Valente  and  Thornton, 
1993 

NC,  USA 

Cotton  -  fertilized 

84  kg  N/ha 

3.9 

Anejaetal.,  1995 

NC,  USA 

Corn  -  fertilized,  173  kg 

N/ha 

8.1 

Anejaetal.,  1995 

Italy 

Agricultural  area,  heavily 
fertilized  with  urea 

91.7 

Kessel  et  al.,  1992 

VA,  USA 

Com  and  Barley  -  fertilized, 
196.4  kg  N/ha 

6.9 

Anderson  and  Levine, 
1987 

NC,  USA 

Sludge  Amended  soil 

57.8 

This  Study 

Grassland  /  Woodland  -  Temperate 

CO,  USA 

Wheat  grass 

8.9 

Delany  et  al,  1986 

CO,  USA 

Grassland 

3.1 

Williams  et  al.,  1987 

Table  3.1.  NO  emission  data  from  different  land  classes.  Source:  Davidson  and 
Kingerlee,  1997. 
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Table  3.2.  Characteristics  of  grab  samples  of  the  soil  at  the  research  site  and  of  the  biosolids  which  are  applied  to  the  field. 


Sedimentation  Tanks/Primary  Clarifiers 


Schematic  of  Waste-Water  Treatment  Facility. 
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Figure  3.2a.  NO  Flux  (ng  N  m’^  s'*)  plotted  versus  Total  Kjeldahl  Nitrogen.  TKN=organic  N+NH3-N+NH, 


♦  Plot  1-Unamended 
■  Plot  2-Unamended 


O  Flux  vs  Time  of  Day  on  an  unamended  experimental  plot  (July  21,  1999). 


Unamended 


s  J..IU  N  6u)  xny  ON 


88 


Time  of  Day 

Figure  3.3b.  Plot  of  NO  Flux  vs  Time  of  Day  on  an  amended  (~  75  ml  of  biosolids)  experimental  plot  (July  23, 1999). 
One  plot  was  amended  with  biosolids  approximately  30  minutes  prior  to  sampling.  The  other  plot  did  not  receive  any 
additional  biosolids. 
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Time  of  Day 

Figure  3.3c.  Plot  of  NO  Flux  vs  Time  of  Day  on  an  amended  (~  75  ml  of  biosolids)  experimental  plot  (July  29,  1999).  One 
plot  was  amended  with  biosolids  approximately  30  minutes  prior  to  sampling.  The  other  plot  did  not  receive  any  additional 
biosolids. 
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Soil  Temperature  (°C) 

Figure  3.4.  NO  flux  plotted  versus  soil  temperature.  NO  flux  has  been  binned  into  1 .5  °C  temperature  spans  and  the 
average  flux  for  each  temperature  span  is  plotted. 


temperature  has  been  binned  into  values  of  15, 25  and  35  °C. 


Chapter  IV.  Modeling  of  Nitric  Oxide  Emissions  from  Biosolid 
Amended  Soils 

Abstract 

Utilizing  a  state-of-the-art  mobile  laboratory  in  conjunction  with  a  dynamic  flow¬ 
through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were  measured 
during  the  summer,  winter  and  spring  of  1 999/2000.  The  field  site  where  these 
measurements  were  conducted  was  an  agricultural  soil  amended  with  biosolids  from  a 
municipal  wastewater  treatment  facility.  The  average  NO  flux  from  this  biosolid 
amended  soil  was  found  to  be  exponentially  dependent  on  soil  temperature  [NO  Flux  = 
1.07  exp(0.14*Tsoii);  R^=0.81].  Comparing  this  relationship  to  results  of  the  widely 
applied  Biogenic  Emissions  Inventory  System  (BEIS)  model  revealed  that  for  this  field 
site,  if  the  BEIS  model  was  used,  the  NO  emissions  would  have  been  underestimated  by 
a  factor  of  26.  Using  this  newly  developed  NO  flux  algorithm,  combined  with  North 
Carolina  Division  of  Water  Quality  statistics  on  how  many  biosolid  amended  acres  are 
permitted  per  county,  county-based  NO  inventories  from  these  biosolid  amended  soils 
were  calculated.  Results  from  this  study  indicate  that  when  the  majority  of  North 
Carolina  is  taken  into  consideration,  NO  emissions  from  biosolid  amended  soils 
represent  approximately  1%  of  the  total  biogenic  NO  emissions.  However,  county-level 
biogenic  NO  emissions  can  increase  by  as  much  as  1 8%  when  biosolid  amended  soils 
are  included  and  depending  on  location  and  time  of  day  the  biogenie  emissions  can  be 
larger  than  the  anthropogenic  NO  emissions.  The  Multiseale  Air  Quality  Simulation 
Platform  (MAQSIP)  was  then  used  to  model  differences  in  ozone  (O3)  and  odd-reactive 
nitrogen  (N)  products  (NOy)  which  might  be  caused  by  differences  in  the  estimated  NO 
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emissions.  Results  showed  that  during  the  daytime,  when  mixing  heights  are  typically  at 
their  greatest,  any  changes  in  O3  or  NOy  concentrations  predicted  by  the  model  were 
insignificant  and  unrecognizable.  However,  in  some  locations  during  late  evening/early 
morning  hours  and  depending  on  how  the  acreages  of  biosolid-amended  soil  are 
dispersed  throughout  the  model,  the  changes  in  concentration  of  these  species  can  be  as 
large  as  10%. 
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Introduction 


Although  it  is  well  established  that  the  primary  source  of  nitric  oxide  (NO) 
emissions  into  the  troposphere  is  anthropogenic  activity,  soil  emissions  can  also  make  a 
significant  contribution  to  the  NO  budget  (Wameck,  2000).  Depending  on  which  NO 
budget  is  considered,  soil  emissions  and  biomass  burning  rank  either  second  or  third 
behind  fossil  fuel  combustion  as  the  largest  source  of  NO  emissions  (Table  1.1).  In  fact, 
in  some  regions  the  soil  contribution  has  been  reported  to  be  approximately  equal  to  the 
emissions  from  anthropogenic  sources  (Yienger  and  Levy,  1995).  Because  of  its  role  in 
tropospheric  ozone  (O3)  formation,  accurate  inventories  of  NO  are  required  to  confidently 
run  air  quality  models  and  to  design  and  implement  O3  control  strategies. 

Williams  et  al.  (1992)  reported  on  the  strong  dependence  of  biogenic  NO 
emissions  on  soil  temperature  and  land  use  type.  The  algorithm  which  Williams  et  al. 
(1992)  developed  is  one  of  the  principal  methods  which  researchers  currently  use  to 
derive  global  inventories  of  soil  NO  emissions  (Yienger  and  Levy,  1995;  Thornton  et  al., 
1997).  In  fact,  many  current  air  quality  models  derive  the  biogenic  NO  input  data  from 
the  Biogenic  Emissions  Inventory  System  (BEIS2)  model  which  is  based  on  a 
temperature  and  land-use  algorithm  proposed  by  Williams  et  al.  (1992).  While  many 
studies  have  been  conducted  on  fallow  soils,  forested  soils,  grassland,  golf  course  soils, 
agricultural  soils,  etc.,  there  are  comparatively  fewer  studies  on  emissions  from  biosolid 
amended  soils  and  therefore  no  land-use  type  for  this  category  currently  exists  in  the 
BEIS2  model. 
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Biosolids  are  often  applied  to  the  soil  for  their  beneficial  nutrient  content  and  as  a 
cost  effective  way  to  dispose  of  this  byproduct  of  the  wastewater  treatment  process  (See 
Chapter  III  for  detailed  discussion  of  biosolids).  Currently  in  the  United  States, 
approximately  6,000,000  metric  tons  are  generated  annually  (Peirce  and  Aneja,  2000). 

In  this  paper,  a  flux  algorithm  which  has  been  previously  developed  (Chapter  III)  to 
estimate  NO  emissions  from  biosolid  amended  soils,  will  be  applied  to  those  acreages  in 
North  Carolina  which  receive  biosolid  applications.  This  updated  NO  inventory  will 
then  be  used  in  the  Multi  Scale  Air  Quality  Simulation  Platform  (MAQSIP)  to  model 
ozone  (O3)  production  on  a  4  x  4  km  grid  scale  resolution  (modified  case).  These  results 
will  then  be  compared  to  estimated  O3  concentrations  using  the  existing  NO  inventory 
(base  case).  The  oxidized  nitrogen  species  (NOy  = 

NO+NO2+HNO3+HONO+NO3+N2O5+HNO4+PAN+RONO2+ROONO2)  and  ratios  of 
these  species  will  also  be  analyzed  to  help  determine  the  end  products  of  the  increased 
NO  emissions.  This  comparison  will  help  to  determine  whether  a  significant  error  may 
exist  in  air  quality  models  by  neglecting  biosolid-amended  soils  as  a  land-use  class  when 
developing  the  NO  emission  inventories. 

Methods  and  Materials 
Instrumentation  and  Flux  Calculation 

The  chamber  design,  associated  mass  balance  equation  and  calibration  procedures 
are  described  in  full  in  Chapter  I. 
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Model 


MAQSIP  is  a  publicly  available,  fully  modularized,  three  dimensional  modeling 
framework  which  has  been  developed  through  a  cooperative  agreement  between  the 
North  Carolina  Supercomputing  Center  and  the  US  Environmental  Protection  Agency’s 
(EPA)  National  Exposure  Research  Laboratory  (NERL).  The  current  MAQSIP  model 
was  developed  with  various  options  for  describing  regional  and  urban-scale  air  quality 
through  physical  and  chemical  processes.  MAQSIP  was,  in  part,  developed  as  a  result  of 
the  dissatisfaction  with  the  black-box  type  models,  which  only  allowed  the  user  to  modify 
the  inputs  or  adjust  a  few  of  the  model  parameterizations  to  optimize  the  performance  of 
the  model  (MAQSIP  Users  Manual,  1998).  The  open  approach  of  MAQSIP  allows 
scientists  to  develop  subroutines  as  new  science  becomes  available,  apply  them  in  the 
model  and  then  evaluate  their  effectiveness.  In  “open”  type  models,  the  rate  constants 
and  chemical  mechanisms  are  evaluated  with  environmental  chamber  data.  Different,  or 
improved  chemical  mechanisms  can  be  substituted  into  the  model  after  evaluation. 
However  the  mechanisms  themselves  are  not  adjusted  once  they  are  placed  within  the 
modeling  framework. 

The  model  is  designed  so  that  fine  grids  (2-6  km)  for  urban  scale  simulations  can 
be  nested  within  coarser  grids  (18-80  km)  for  regional  scale  simulations,  and  consists  of  a 
vertical  domain  of  up  to  30  layers.  The  spatial  resolution  used  in  this  study  was  a  4  x  4 
km  grid  cell,  38  m  high,  where  one  concentration  per  species  was  calculated  for  each  grid 
cell.  The  episode  being  examined  is  a  13-day  period  from  June  19  -  July  1, 1996,  and  the 
model  domain  (480  km  x  276  km)  covers  the  majority  of  North  Carolina  (see  Figure  4.1). 
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Mixing  ratios  of  species  are  determined  in  the  model  by  solving  the  mass 
conservation  equation  for  a  given  volume.  The  equation  can  be  represented  as: 

^  =  -V-(VC)  +  V-(KVC)  +  R  +  S  (1) 

dt 

where 

C :  mixing  ratio  of  species  of  interest 

V  :  velocity  vector  (3 -Dimensional  at  each  point  in  model  domain) 

K :  eddy  diffusivity  used  to  parameterize  fluxes  of  tracer  species  at  subgrid-scale 

S :  losses  due  to  sources  and  sinks  (to  include  dry  deposition) 

R  :  changes  in  concentration  due  to  chemical  reactions 

(MCNC,  2001;  Arya,  1999) 

For  this  simulation  MAQSIP  was  configured  to  predict  lower  tropospheric  ozone 
and  NOy  species  using:  Carbon  Bond  IV  chemical  mechanism,  K-Theory  for  turbulent 
redistributions  of  pollutants  in  the  vertical  and  a  dry  deposition  scheme  (McHenry  et  al., 
1999). 

Carbon  Bond  IV  Mechanism  (CBM  IV) 

Computational  power  currently  limits  the  ability  to  explicitly  represent  the 
chemistry  of  the  hundreds  of  organic  species  in  the  troposphere  (Gery  et  al.,  1989). 
Therefore,  to  be  able  to  represent  these  reactions,  the  common  chemistry  of  the  reactive 
hydrocarbons  are  grouped  together.  Two  different  grouping  approaches  have  been 
attempted,  namely  structural  and  molecular  grouping,  where  the  former  involves 
grouping  according  to  the  type  of  bond  (ie,  single  carbon  bond,  double  carbon  bond,  etc.), 
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and  the  latter  groups  common  reactions  of  entire  molecules  (Gery  et  al.,  1989).  Grouping 
by  bond  type  has  been  found  to  be  a  favorable  approach  as  there  are  many  fewer 
categories  to  consider  as  opposed  to  the  large  number  of  organic  species  in  the 
atmosphere.  The  carbon  bond  approach  categorizes  the  species  as  follows: 

1)  inorganic  species 

2)  organic  species  that,  because  of  their  unique  chemistry,  are  treated  explicitly 

3)  organic  species  that  are  represented  by  carbon  bond  surrogates 

4)  organic  species  that  are  represented  by  molecular  surrogates 
(Gery  etal.,  1989). 

This  approach,  which  contains  over  200  chemical  reactions,  was  found  to  be  too 
computationally  extensive  and  was  paired  down  to  the  approximately  80  reactions  now 
present  in  CB  IV.  Both  the  larger  mechanism  and  CB  IV  were  tested  in  smog  chamber 
studies  and  both  performed  well  without  any  consistent  bias  between  the  two  models. 

The  current  Carbon  Bond  IV  principle  mechanism  is  listed  in  Table  4.1. 

K  Theory  Scheme  for  Turbulent  Vertical  Re-Distribution  of  Pollutants: 

Accurate  representation  of  the  atmospheric  boundary  layer  is  important  to  air 
quality  modelers,  as  conditions  in  this  layer  dictate  the  vertical  extent  to  which  the 
pollutants  can  be  mixed.  The  model  uses  the  general  form  of  the  following  dynamic  and 
thermodynamic  equations: 
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dq,  dq^ 

=  -u  -  V 

dx  dy 


u :  eastward  component  of  the  wind 
V :  northward  component  of  the  wind 
Mg :  eastward  component  of  the  geostrophic  wind 
Vg :  northward  component  of  the  geostrophic  wind 


(4) 

(5) 


(  ) :  turbulent  fluxes 

z :  altitude 

/ :  Coriolis  parameter 

q,, :  mixing  ratio  of  water  vapor 

T :  temperature 

Explicitly  solving  for  the  turbulent  fluxes  is  both  computationally  expensive  and 
time  consuming,  making  it  impractical  for  real-time  air  quality  forecasting.  Therefore  to 
solve  for  these  turbulent  fluxes,  they  are  assumed  to  be  related  to  the  mean  gradients,  as 
initially  proposed  by  J.  Boussinesq  in  1877  (Arya,  1999),  and  solved  using  gradient 
theory. 

Dry  Deposition  Scheme 

Dry  deposition  is  defined  as  airborne  gaseous  and  particulate  matter  being 
transferred  to  the  earth’s  surface  (Arya,  1999).  Early  air  quality  models  used  constant 
deposition  velocities  (Vj)  in  either  time  or  space,  which  various  researchers  have  shown 
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to  be  a  poor  assumption  (Hicks  et  al.,  1982;  Wesely  et  al.,  1985).  These  researchers  have 
shown  that  deposition  of  various  trace  gases  and  particulates  have  considerable 
variability  depending  on  meteorology,  land-use  type,  insolation  and  season,  particulate 
size,  etc.  Deposition  velocity  is  used  to  calculate  the  flux  of  the  pollutant  in  question  to  a 
surface  through  the  following  expression: 


The  constant  (V^)  approach  led  to  inaccurate  deposition  estimates  of  trace  gases 

and  particulates  when  applied  to  large  regions  with  different  topography  and 
meteorological  conditions.  Walcek  et  al.,  (1986)  have  described  a  method  for  computing 
maps  of  deposition  velocities  which  account  for  meteorology  and  land  use  type.  In  their 
model,  they  use  fine  resolution  meteorological  model  output  and  land  use  data  from  the 
United  States  Environmental  Protection  Agency  to  adjust  the  deposition  velocities  that 


are  then  used  in  the  MAQSIP  model. 


Emissions  Data 

Emission  inputs  into  MAQSIP  are  provided  by  the  Sparse  Matrix  Operator  Kernel 
Emissions  (SMOKE)  modeling  system.  The  submodels  for  SMOKE  are  biogenic, 
mobile,  area  and  point  emissions  (Figure  4.2),  which  eventually  are  merged  to  produce 
emissions  that  are  ready  for  input  into  the  models  timeline  and  grid  system.  The  biogenic 
emissions  processing  submodel  is  essentially  the  Biogenics  Emissions  Inventory  System 
2  (BEIS2)  model  which  has  been  slightly  modified  to  be  compatible  with  SMOKE.  The 
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biogenic  sources  include  Volatile  Organic  Compounds  (VOCs)  from  plants  and  trees  and 
NOx  emissions  from  soil  processes. 

Prior  to  any  emission  inventories  being  produced,  accurate  representation  of  the 
meteorological  conditions  must  be  produced  as  many  of  the  emission  processes  are 
dependent  on  the  temperature  fields.  Meteorological  files  for  MAQSIP  are  provided  by 
the  Pennsylvania  State  University  (PSU)  and  National  Center  for  Atmospheric  Research 
(NCAR)  mesoscale  model  5  (PSU/NCAR  MM5).  A  schematic  of  the  inputs  and  outputs 
of  MAQSIP  can  be  seen  in  Figure  4.3. 

Currently,  biogenic  nitric  oxide  emissions  are  estimated  in  BEIS2  through  a 
temperature  and  land-use  relationship  described  in  Williams  et  al.  (1992),  where 

NO  Flux  (ng  N  m'^  s’*)  =  A*Exp(0.07 1  *Ts)  (7) 

The  coefficient.  A,  is  an  emission  factor  that  is  experimentally  derived  and  based 
on  land  use.  Ts  (”C)  is  soil  temperature  which  is  computed  via  air  temperature  (Ta  in  "C) 
through  the  following  algorithm,  given  by  Williams  et  al.  (1992): 


Grasslands 

Ts=0.66Ta+8.8 

(8) 

Forests 

Ts=0.84Ta+3.6 

(9) 

Wetlands 

Ts=0.92Ta+4.4 

(10) 

Agriculture 

Ts=0.72Ta+5.8 

(11) 

For  each  of  the  4  x  4  km  grid  cells  (=1600  ha)  in  the  model  domain,  statistics  are 
available  which  describe  the  land-use  in  that  grid  cell.  An  example  of  these  statistics  are 
shown  in  Table  4.2.  Each  4  letter  code  corresponds  to  a  crop,  tree  or  land  use  type  and 
has  an  associated  acreage  and  an  emissions  factor  that  is  used  for  calculating  the  NO  flux. 
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In  general,  the  emission  factors  will  be  highest  for  those  crops  which  receive  the  highest 
applications  of  nitrogen  fertilizer  (ie,  corn,  cotton,  wheat),  as  increased  N  fertilization  has 
been  shown  to  increase  NO  emissions  (Sullivan  et  ah,  1996;  Roelle  et  ah,  1999).  While 
this  land  use  approach  is  a  significant  improvement  over  the  earlier  methods  which 
assumed  constant  NO  emissions  across  all  crops  or  forest  types,  it  still  fails  to  capture  the 
emissions  from  all  land-use  types.  For  example,  if  a  rye  crop  were  fertilized  with 
biosolids,  it  would  be  treated  in  the  model  as  having  an  NO  factor  of  12.8  pg  m  ^  h  ^  An 
earlier  study  (Chapter  III)  showed  that  a  rye  crop  amended  with  biosolids  generates  an 
NO  emission  factor  closer  to  the  value  for  com  (550.5  pg  m'^  h'*).  Using  daily  averaged 
temperatures.  Figure  4.4  shows  how  the  BEIS2  model  would  have  underestimated  the 
NO  emissions  during  all  seasons.  Therefore  this  study  proposes  to  modify  the  current 
NO  inventory  (base  case)  so  that  these  biosolid  amended  acreages  are  represented  in  the 
biogenic  NO  inventory  (modified  case). 

Source  Apportionment 

The  process  of  applying  biosolids  to  soils  in  North  Carolina  requires  a  permit 
from  the  North  Carolina  Division  of  Water  Quality.  The  Division  of  Water  Quality 
maintains  a  database  of  all  existing  permits  on  a  county-basis  and  although  restrictions  do 
exist  that  dictate  where  and  when  the  biosolids  can  be  applied,  there  is  no  restriction  upon 
which  crop  they  are  applied  (Barnett,  2000).  For  the  purposes  of  this  study,  the  acreages 
of  biosolid  amended  soils  within  each  county  are  assumed  to  be  equal  to  the  existing 
permits  issued.  Figure  4.5  represents  the  area  of  North  Carolina  within  the  model  domain 
and  the  respective  aereages  of  biosolid  amended  soils  within  each  county. 
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It  is  not  possible  to  ascertain  where  within  each  county  the  biosolids  are  applied, 
therefore  the  procedure  used  to  apportion  the  bio  solids  within  a  county  was  to  distribute 
them  so  that  each  4  x  4  km  grid  cell  within  the  county  had  the  same  percentage  of 
biosolid  amended  soil  in  its  respective  cell.  County  boundaries  do  not  follow  the  same 
boundaries  as  the  grid  cells,  and  therefore  there  were  often  cells  that  contained  acreages 
in  two  counties.  In  these  instances,  the  county  surrogate  data,  which  identifies  the 
fraction  of  a  grid  cell  within  a  county,  were  used  to  apportion  the  biosolids  accordingly. 

It  is  estimated  that  V2  of  the  total  cost  of  wastewater  treatment  is  tied  to  disposal 
(Viessman  and  Hammer,  1993).  Consequently,  transporting  distances  will,  in  part, 
dictate  where  these  biosolids  are  applied  and  evenly  distributing  the  biosolids  throughout 
the  county  may  not  be  the  most  likely  scenario.  Therefore,  in  another  scenario  we  have 
also  analyzed  the  impact  of  concentrating  the  biosolids  in  only  one  area  of  the  county. 
These  two  scenarios  should  be  viewed  as  bounds  which  can  be  used  to  assess  the 
significance  of  biosolid  application  in  terms  of  changes  in  ozone  or  NOy  concentrations. 

The  new  biosolid  amended  acreage  for  each  grid  cell  could  not  simply  be  added 
back  into  Table  4.2,  because  this  would  cause  the  total  acreages  within  each  grid  cell  to 
be  too  large  (ie,  total  acreages  in  each  cell  would  now  be  1600  ha  +  biosolid  amended 
acres).  Because  the  crops  or  land-use  types  receiving  the  biosolids  could  not  be 
identified,  the  biosolid  category  was  substituted  for  the  Mscp  (Miscellaneous  Crop) 
category  in  all  instances.  For  example  (in  Table  4.2),  if  grid  cell  001  001  was  determined 
to  contain  50  ha  of  biosolid  amended  soils,  then  category  Mscp  would  be  reduced  to 
484.8  ha  and  a  new  category  of  biosolid  amended  soils  would  be  added  with  an  acreage 
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of  50  and  the  emission  factor  for  the  biosolid  amended  soil  (calculated  previously)  would 
be  applied  to  this  acreage.  It  was  now  possible  to  account  for  the  biosolid  amended  soil 
and  determine  a  new  NO  inventory  for  each  grid  cell  in  the  model  domain. 

The  increase  in  biogenic  emissions  of  NO  resulting  from  the  biosolids 
applications  was  determined  for  each  of  the  counties  in  the  model  domain  to  produce  a 
modified  NO  inventory.  The  average  increase  of  biogenic  NO  in  each  of  the  North 
Carolina  counties  in  the  model  domain  can  be  seen  in  Figure  4.6.  For  example,  in  Gaston 
County  NC,  approximately  18%  of  the  biogenic  NO  is  due  to  these  biosolid  amended 
soils,  and  in  Brunswick  County,  Forsyth  County,  Granville  County  and  Wake  County,  the 
contribution  from  biosolid  amended  soils  is  10%,  14%,  9%  and  7%  respectively. 
Similarly,  applying  all  the  biosolids  in  concentrated  areas  resulted  in  1-2  grid  cells  being 
modified  in  each  of  the  5  counties  with  the  highest  biosolid  land  application  rate.  The 
cells  which  were  selected  and  the  factors  by  which  the  biogenic  NO  emissions  increased 
in  each  of  these  cells  can  be  seen  in  Figure  4.7. 

Results  and  Discussion 

Figure  4.8  shows  the  relative  contributions  to  the  NO  inventory  from  biogenic  and 
anthropogenic  sources  for  the  entire  model  domain  and  for  the  5  counties  with  the  largest 
application  rates  of  biosolids.  Over  the  entire  model  domain,  the  biogenic  emissions 
represent  approximately  5%  of  the  anthropogenic  emissions.  Surprisingly,  for  the  5 
counties  with  the  greatest  applications  of  biosolids  (with  the  exception  of  Granville 
County)  the  biogenic  emissions  make  up  an  even  smaller  fraction  than  the  overall  model 
average  of  5%  (Brunswick  County~l%;  Gaston  County<l%;  Granville 
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County~l  l%;Wake  County~2.5%;  Forsyth  County~l%).  While  this  may  at  first  seem 
contradictory,  this  relationship  can  be  explained  by  a  combination  of  population  centers 
and  power  generation  locations  (Figure  4.9). 

Biosolids  are  produced  as  a  byproduct  of  the  municipal  waste-water  treatment 
process  (Chapter  III),  therefore  those  areas  with  the  densest  populations  will  produce  the 
greatest  quantity  of  biosolids.  In  these  more  densely  populated  counties  (or  adjacent 
counties),  if  land  application  is  the  chosen  method  of  handling  (as  opposed  to  the  costlier 
methods  of  land-filling  or  incineration),  there  will  be  greater  land  application  rates  (ie, 
Gaston  County  and  Wake  County,  Forsythe  County).  These  densely  populated  urban 
areas  have  greater  vehicular  activity,  and  in  some  instances  have  large  power  production 
facilities  and  therefore  have  NO  budgets  dominated  by  anthropogenic  sources.  This 
effect  can  clearly  be  seen  in  Figure  4.10,  which  shows  the  strong  mobile  NOx  source 
strength  near  Gaston  County,  Wake  and  Forsyth  Counties,  and  to  a  lesser  extent, 
Brunswick  and  Granville  Counties.  Although  Brunswick  County  is  less  affected  by 
mobile  sources,  the  biogenic  NO  contribution  is  still  overshadowed  by  the  utility  sources 
(Figure  4.9). 

When  the  data  from  Figure  4.8  is  hourly  averaged  over  the  entire  model  domain 
(Figure  4.1 1),  several  interesting  observations  become  apparent.  First,  the  peaks  in  both 
curves  are  in  phase,  occurring  in  the  late  afternoon,  which  is  to  be  expected  based  on  the 
biogenic  temperature  dependence  and  vehicle  driving  patterns  which  reach  a  maximum 
during  the  late  afternoon.  The  two  relatively  lower  peaks  in  anthropogenic  emissions 
which  occur  at  approximately  75,100,  250  and  270  hours  are  during  weekends  when 
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vehicular  traffic  is  generally  at  a  minimum  (Figure  4. 12).  What  is  most  revealing  in 
Figure  4.1 1  is  the  large  standard  deviation  in  anthropogenic  emissions  indicating  that 
relatively  few  of  the  grid  cells  in  the  model  domain  are  responsible  for  the  majority  of  the 
emissions.  Upon  inspection  of  the  individual  grid  cells  it  was  found,  as  expected,  that 
those  grid  cells  responsible  for  the  majority  of  the  emissions  were  located  in  and  around 
where  the  power  plants  were  located  (Figure  4.9).  Furthermore,  the  anthropogenic 
emissions  vary  by  as  much  as  a  factor  of  four  throughout  the  diurnal  profile,  whereas  the 
biogenic  emissions  vary  by  only  30%.  The  combination  of  these  two  effects  (location 
within  the  model  domain  and  time)  indicates  that  on  a  localized  scale,  the  biogenics  could 
represent  a  significant  fraction  of  the  total  NO  emissions  and  may  therefore  be  important 
in  terms  of  developing  accurate  air  quality  models. 

An  example  of  the  temporal  and  spatial  relationship  can  be  seen  in  Figure  4.13  (a- 
b)  which  shows  the  variation  in  the  ratio  of  biogenic  NO  emissions  to  the  anthropogenic 
emissions.  In  Figure  4.13a  (weekday,  mid-aftemoon),  the  biogenics  in  Granville  County 
are  at  most  40%  of  the  anthropogenic  source  strength,  and  only  in  a  few  remote  locations. 
The  white  line  extending  from  the  southwest  comer  of  the  county  to  the  northeast  reflects 
Interstate  85  and  the  dominance  of  mobile  sources.  Figure  4.13b  represents  a  weekend 
early-morning  episode  which  shows  how  the  biogenics  now  begin  to  dominate  the 
anthropogenics,  in  some  instances  by  an  order  of  magnitude,  throughout  the  more  rural 
portions  of  the  county.  Similar  analyses  were  conducted  for  other  counties  with  the  same 
results  and  therefore  in  the  analysis  of  the  MAQSIP  model  mns,  those  areas  and  time 
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periods  in  which  biogenics  make  the  most  significant  contribution  to  the  total  NO  budget 
will  be  most  closely  scrutinized. 

In  Figure  4.14  (a-c),  the  change  in  O3  concentration  (A[03])  (Modified  Case  - 
Base  Case)  in  pptV  is  plotted  beginning  the  evening  (8:00  pm)  of  Wednesday,  June  26 
and  ending  after  sunrise  (8:00  am)  on  the  following  day.  In  Figure  4.14a,  the  areas  of 
greatest  A[03]  are  observed  to  coincide  with  the  areas  of  greatest  biosolid  application 
(Figure  4.5).  Furthermore  these  areas  are  all  negative  in  value,  as  shown  by  the  scale  to 
the  side  of  the  figure,  indicating  ozone  depletion.  During  the  early  morning  (Figure 
4.14b),  the  same  pattern  of  O3  depletion  occurs  but  in  greater  magnitude.  As  daylight 
approaches  (Figure  4.14c),  the  same  regions  of  maximum  biosolid  application  continue  to 
have  the  largest  -A[03],  however  the  magnitudes  begin  to  get  smaller. 

Similarly  to  the  weekday  episode,  a  weekend  morning  (Figure  4.15a)  shows  O3 
being  consumed  in  the  largest  quantities  where  the  biosolid  application  is  the  greatest. 
Interestingly,  while  the  large  urban  centers  continue  to  indicate  ozone  depletion,  some  of 
the  more  rural  areas  which  also  had  high  biosolid  application  rates  are  beginning  to  show 
modest  increases  in  O3  production.  This  relationship  could  confirm  that  some  areas  in 
North  Carolina  may  be  NOx  limited,  meaning  an  increase  in  the  NOx  emissions  results  in 
increased  [O3].  The  early  morning  plot  (Figure  4.15b)  shows  approximately  the  same 
trends  as  the  weekday  plot,  where  regions  of  greatest  biosolid  application  contribute  the 
most  to  the  O3  depletion.  Throughout  the  daylight  hours  (~  8:00  am  -  8:00  pm),  no 
consistent  changes  in  [O3]  (other  than  patterns  similar  to  the  one  shown  in  Figure  4.15a) 
were  observed  on  any  of  the  modeled  days.  This  same  diurnal  pattern  of  decreasing  [O3] 
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throughout  the  evening,  reaching  a  maximum  -A[03]  before  sunrise  and  then  gradually 
returning  to  unrecognizable  changes  throughout  the  daylight  hours  was  evident  during  all 
of  the  modeled  days  and  can  be  explained  by  the  combined  effects  of  meteorology  and 
chemistry. 

Figure  4.16  shows  the  diurnal  variation  of  the  mixed  layer,  depicting  how  it  grows 
to  its  maximum  in  the  afternoon  coinciding  with  the  maximum  daytime  heating.  It  is 
therefore  expected  that  the  greatest  impact  of  any  increased  NO  emissions  would  occur 
during  the  morning/evening  hours  when  they  are  confined  to  the  smallest  mixing  volume. 
Any  increased  emissions  during  the  daytime  would  be  rapidly  diluted  as  the  species  mix 
through  a  significantly  greater  volume  of  the  troposphere.  The  influence  of  chemistry  is 
based  on  the  following  equation:  NO  +  O3— >  NO2+O2,  which  acts  to  remove  O3 
throughout  the  night.  During  the  evening  when  the  boundary  layer  is  confined,  and  the 
biogenic  NO  source  strength  is  increased  as  a  result  of  the  biosolids,  ozone  consumption 
is  increased  and  most  evident  in  areas  of  greatest  biosolid  application,  as  shown  in 
Figures  4.14  and  4.15.  However,  as  was  shown  in  Figure  4.15a,  there  may  be  a  few 
hours  prior  to  the  boundary  layer  reaching  its  maximum  that  the  increased  NO  results  in 
slight  increases  of  [O3]. 

It  should  be  pointed  out  that  all  of  the  A[03]  data  presented  in  Figures  4.14  and 
4.15  are  in  units  of  pptV  whereas  ambient  concentrations  of  O3  during  this  same  time 
period  are  typically  on  the  order  of  1-100  ppbV.  At  no  time  during  the  daylight  hours 
does  the  A[03]  approach  the  same  order  of  magnitude  as  the  ambient  concentrations  and 
is  typically  several  orders  of  magnitude  smaller  indicating  that  during  periods  of  greatest 
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mixing  the  increase  in  NO  has  negligible  consequences  on  overall  ozone  formation.  The 
simulation  with  the  overall  greatest  change  in  [O3]  is  plotted  in  Figure  4.17  along  with  a 
time  series  of  the  data  throughout  a  diurnal  cycle.  This  relationship  highlights  that 
daytime  ozone  concentrations  are  unaffected  by  the  increased  biogenic  NO  emissions  and 
that  any  change  at  night  consists  of  slight  ozone  depletion  (<1%). 

Applying  all  biosolids  in  a  concentrated  region  of  the  county,  as  described  in 
Figure  4.7,  produces  more  pronounced  results.  Figure  4.18  shows  the  percentage  change 
in  [O3]  and  the  corresponding  time  series  of  this  change  throughout  a  diurnal  cycle  at  a 
Gaston  County  location.  Like  the  other  plots,  the  biggest  change  in  [O3]  occurs  at  night 
and  consists  of  ozone  depletion.  Whereas  in  the  earlier  method  of  evenly  distributing  the 
biosolids  throughout  the  county  resulted  in  changes  of  less  than  1%,  the  highly 
concentrated  biosolid  areas  now  see  ozone  being  depleted  by  as  much  as  11%.  Further, 
for  the  first  time,  consistent  trends  of  increased  ozone  production  (>1%)  during  afternoon 
hours  were  evident,  albeit  at  values  less  than  3%. 

In  addition  to  O3  production,  as  NO  is  oxidized  it  also  gives  rise  to  new 
compounds  which  are  often  grouped  together  in  the  family  called  NOy.  This  family  of 
odd-reactive  nitrogen  species  consists  of  NOx  (NO+NO2),  nitric  acid  (HNO3),  nitrous 
acid  (MONO),  the  nitrate  radical  (NO3),  dinitrogen  pentoxide  (N2O5),  peroxynitric  acid 
(HNO4),  peroxyacetyl  nitrate  (PAN)  (RC(0)00N02),  alkyl  nitrates  (RONO2)  and 
peroxyalkyl  nitrates  (ROONO2).  Figure  4.19  shows  the  percentage  increase  in  NOy 
concentrations  between  the  modified  case  and  the  base  case  where  the  bioslids  have  been 
evenly  distributed  throughout  all  grid  cells  in  the  county.  Also  included  in  Figure  4.19  is 
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the  time  series  of  the  data  with  the  greatest  percentage  increase  throughout  a  diurnal  cycle 
which  was  observed  to  be  in  Brunswick  County.  As  shown  in  the  graph,  any  change  in 
[NOy]  during  daylight  hours  was  imperceptible  within  several  orders  of  magnitude,  and 
the  maximum  change  throughout  the  evening  was  an  increase  in  [NOy]  of  ~  2%. 

NOy,  which  is  a  quasi-conserved  quantity,  is  of  significant  atmospheric  interest 
because  through  ratio  analysis  it  is  possible  to  examine  the  fate  of  the  increased  nitrogen 
and  also  to  determine  relative  aging  of  an  airmass  (Roberts,  1995).  For  example.  Figure 
4.20  is  a  plot  of  the  ratio  of  NOx/NOy  which  shows  that  near  the  large  urban  centers,  the 
majority  of  NOy  exists  as  NOx.  However,  in  more  rural  areas,  the  NOx  emitted  into  the 
atmosphere  upstream  begins  to  be  converted  into  organic  and  inorganic  nitrates  and 
therefore  makes  up  a  smaller  fraction  of  the  NOy.  During  the  nighttime,  the  increase  in 
[NOy]  observed  in  Figure  4.19  can  almost  solely  be  attributed  to  the  increased  NO,  as 
investigation  of  the  changes  in  the  ratios  of  all  the  individual  species  to  NOy  resulted  in 
changes  of  less  than  .01%. 

Conclusions  and  Recommendations 

Comparing  the  NO  flux  algorithm  developed  in  an  earlier  study  (Chapter  III)  to 
results  using  the  BEIS2  algorithm  (Base  case)  revealed  that,  for  one  particular  field  site 
the  BEIS2  model  would  have  underestimated  NO  emissions,  on  a  yearly  average,  by  a 
factor  of  26.  Applying  this  new  flux  algorithm  to  the  biosolid  acreage  data,  a  modified 
NO  inventory  was  developed  (Modified  case).  It  should  be  noted  that  the  data  used  to 
produce  this  revised  inventory  was  from  only  one  field  site  and  has  been  extended  to  all 
soils  throughout  the  model  domain  receiving  the  biosolids.  The  exponential  temperature 


no 


relationship  reported  here,  however  has  been  consistently  reported  throughout  the 
literature  and  therefore  provides  some  basis  for  this  study  to  be  extended  to  these  various 
soil  and  crop  types  throughout  North  Carolina.  It  should  not  be  assumed,  however  that 
this  temperature  dependence  can  be  extended  to  all  temperatures,  as  temperatures  outside 
the  range  of  1 5-35  °C  are  often  foimd  to  alter  the  often-cited  exponential  relationship. 

The  results  of  this  research  indicate  that  on  a  broad  scale  (entire  model  domain  ~ 
132,000  km^),  the  contribution  of  NO  from  biosolid  amended  soils  to  the  total  biogenic 
emissions  inventory  in  North  Carolina  is  approximately  1%.  It  can  be  argued  that  when 
the  entire  model  area  is  taken  into  consideration,  biogenic  NO  emissions  are  less  than  5% 
of  the  anthropogenic  emissions,  and  therefore  any  modest  increase  in  the  biogenic  source 
strength  will  likely  have  negligible  consequences  on  tropospheric  air  quality.  However, 
the  majority  of  the  anthropogenic  emissions  are  concentrated  in  or  around  areas  which 
contain  large  power  plants  or  large  urban  centers.  Consequently,  in  these  industrial  and 
urban  areas,  biogenics  are  a  very  small  fraction  of  the  total  NO  inventory.  In  the  more 
remote/rural  areas  however,  biogenics  can  be  as  much  as  an  order  of  magnitude  larger 
than  the  anthropogenic  NO  emissions.  Therefore  an  underestimation  of  the  biogenic  NO 
in  these  remote  regions  would  result  in  a  significant  bias  in  the  emissions  inventory  for 
these  areas. 

Results  from  the  two  model  scenarios  (modified  case  and  base  case)  revealed  that 
any  increased  NO  from  biosolid  amended  soils  produced  A[03]  several  orders  of 
magnitude  smaller  than  background  concentrations  during  daytime  hours.  During  the  late 
evening/early  morning  hours  when  the  mixing  volumes  are  at  their  smallest,  A[03]  were 
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found  to  be  reduced  at  most,  by  approximately  1%.  During  one  particular  episode,  both 
ozone  production  and  ozone  depletion  were  observed  during  the  same  time  period, 
possibly  indicating  that  some  of  the  more  remote  areas  of  North  Carolina  are  NOx 
limited.  NOy  species  showed  slightly  larger  changes,  although  most  of  the  increase  could 
be  attributed  to  the  increased  NO  rather  than  to  any  of  the  oxidized  products. 

Depending  on  how  uniformly  the  biosolids  are  distributed  throughout  the  model 
domain  can  cause  larger  impacts  in  the  model  results.  By  applying  the  biosolids  in 
concentrated  areas  of  the  county  resulted  in  ozone  depletion  of  as  much  as  1 1%  and  for 
the  first  time,  consistent  trends  of  increased  ozone  production  during  afternoon  hours  was 
evident.  While  both  approaches  (concentrating  biosolids  in  one  area  versus  evenly 
distributing  them  throughout  the  county  grid  cells)  are  inherently  flawed,  they  do  acts  as 
bounds,  within  which  the  true  situation  likely  exists.  More  importantly,  this  study 
addresses  the  importance  of  biosolid  amended  soils  as  a  land-use  class  and  that  more 
detailed  monitoring  is  required  so  that  they  can  be  accounted  for  in  future  emission 
inventories. 
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//  Updated  Corbon  Bond  Four  Prinoiple  Mechanism  Uersion  99 

// 

//  I  *=  INORGANIC  CHEMISTRV  =«==s==s====»==»=«=a=«ass*«ss«c»a*===B*»==s 
//  a)  N02  photolysis - — - 

NAMES 

PhotoRatelDs  +«  {  N02_to_03P  }; 

Rllal]=  N02  -hv~>  NO  +  03P 

R[Ia2]=  03P  +  02  +  M - >  03  +  M 

R[Ia3]«  03  +  NO  - >  N02  +  02 

R[Ia4]»  03P  +  N02  >  NO  +  02 

R{Ia5l«  03P  +  N02  >  N03  +  02 


Rlla6]«  03P  +  NO  - >  N02 


R{Ia7l=  NO  +  NO  +  02  - >  2.0*NO2 

//  b)  N03  CHEMISTRV  - 

NAMES 

PhotoRatelDs  +=  {  N03_toJH0,  N03„toji02  }; 

Rtibll  «  03  +  N02 - >  N03  +  02  8  1.2E-13*EXP(-2450.0/TK) J 

Rilb2]  -  N03  -hv->  NO  +02  8  J  [N03„to_N0l ; 

Rllb2b]=  NOS  -hv->  N02  +  03P  8  j  [N03_to Ji02l ; 

Rllb3]  =  NOS  +  NO  >  2.0*NO2  8  1.S0E-11*EXP(170.0AK)  J 

Rllb4]  «  NOS  +  N02 - >  NO  +  N02  +  02  8  4.50E-14*EXP("1260.0AK); 

Rtlb3f]*=  N03  +  N02  HI— >  N205  8  TROE(2.20E-38*T_300^~3.9, 

1.50E-12*T  300^-0.7, 
biM],  0.6)  ; 

Rtlb5r]«  N205  >  NOS  +  N02  8  k[Ib5f]/C2.7E~27*EXPC  11000. 0AK))  | 

//  IbSf  NnSfi97,  T2; 

//  IbSr  Ke*2.7E“27*EXP(11000A) 

Rllb7la  N205  +  H20  >  2.0+HNO3  8  l.SE-Sl;  //  homogeneous  rate  only 

//  c)  OZONE  photolysis  - 

NRMES 

PhotoRatelDs  +-  {  03.to_03P,  03.toJJlD  }; 

R[Ioll  =  03  'hu“>  03P.  +  02  8  j  l03_to_03Pl  j 

R  [Ic2i  *  03  H^v->  OID  +  02  8  j  [03_to_01Dl  \ 

RlIoS]  -  OID  +  M  >  03P  +  M  8  1.92E-11*EXP(  12B.0AK)  J 

//  awe  of  N2  and  02  rates 

R[Ic4]  -  OID  +  H20 - >  2.0*OH  8  2.20E-10; 

R(Ic5]  -  03  +  OH  >  H02  +  02  8  l.eBE-12*EXP(-940.0AK)  \ 

RlIo6]  -  03  +  H02 - >  OH  +  2.0  *  02  8  l.i0E~14*EXP(-5B0.0AK>  J 

//  d)  HOMO  CHEMISTRV  - 

NAMES 

PhotoRatelDs  +«  {  H0N0_to_0H  )5 

Rlldl]  *  NO  +  N02  +  H20 - >  2.0+HONO  8  4.4E-40  f 

Rlld2]  “  HONO  +  HONO  >  NO  +  N02  +H20  8  i.0E"’20  ; 

Rlld3l  -  OH  +  NO  HI— >  HONO.  8  TROE(7.00E-31+T_300^-2.6, 

3.60E-lin_300^-0.1, 
blM],  0.6)  } 

Rlld4l  *  HONO  -hw->  OH  +  NO  8  j  iHONO.toJDH)  ; 

RtldSl  *  OH  +  HONO  >  N02  +  H20  8  1.80E-11*EXPC-390.0AK) J 

//  e)  N0/N02  with  H02 - ^ - 

Rllel]  »  H02  +  NO - >  OH  +  N02  8  3.50E“12*EXP(250.0AK)  J 

R[l62fl  *  H02  +  N02  -K— >  PNR  8  TROE(l.B0E-31*Tj300^-3.2, 

4.70E-12*T_300^-1.4, 
b[Ml,  0.6) J 

R[Ie2r]  =  PNR  HI— >  H02  +  N02  8  krie2fl/(2.1E-27+EXPCl0900.0AK)); 

//  Ie2f  NRSR97,  T2;  Ie2r  Ke=2.1E“27+EXP(10900A) 

R[Ie3]  *  OH  +  PNR - >  N02  +  H20  +  02  8  1.30E-12*EXP(3B0.0AK)  J 


Table  4.1.  Carbon  Bond  IV  Principle  Mechanism. 


a  j [N02_to_03P] ; 
a  6.0E»34*T_300^-2.3; 

8  2.0E-12*EXP(“1400.0AK)i 

//  Ia2,  Ia3  NRSR97,  Ti 
8  6.5E-12*EXPC120.0/TK); 

//  Ia4  NflSR97,  Tl 
8  TROE(9.00E-32*T_300^-2.0, 
2.2E-il, 
b[M]»  0.6); 

8  TROE(9.00E-32*T_300^-1.5, 
3.0E-11, 
biMl,  0.6); 

//  laS  and  loB  Nfl5fl97,  T2 
8  3.30£~39*EXP(530.0AK); 

//  Ia7  IUPflC97 
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//  f)  HD2  TERMINRTION  REACTIONS 


mnES 

PhotoRotelDs  +=  {  H202  lo  OH  }; 


R[Ifi1  «  OH  +  ND2 

R[If2]  «  OH  +  HN03 

R[If3]  -  H02  +  H02 

R[If6]  =  H202 

R[If7l  »  OH  +  H202 

R[If8j  =  OH  +  H02 


-M“->  HN03 
- >  N03  +H20 

- >  H202  +  02 

-hv->  2.0*OH 

- >  H02  +  H20 

- >  02  +  H20 


fi  TROE(2.6E-30*T  300^-2.9, 

7.5E-11*T  300^-0.6,  //lUPflC  B:07 
b[l1I,  0.41)j 

&  7.20E-15*EXP(705.0AK)  + 

U1HW(1.90E-33*EXPC  725. 0AK), 

4 . 10E‘-16’*EXP  ( 1440 . 0 AK) , 
bin]) I  //  NflSfl97 
a  (  2.2E-13*EXP(  600.0/TK)  + 

1.9E-33*EXP(  9B0.0AK)  ♦  b[M]  )  ♦ 

(  1.0  + 

1.4E-21*EXP(2200.0AK)  ♦  blH20l  ); 
//  lUPflC  97 
a  j [H202  lo_OHl ; 

a  2.g0E-i2*EXp(-i6B.0AK)  ; 
a  4.80E-11*EXP(  250. 0AK)  ; 


//  IfS—IfS  NRSR97,  Tl,  not®  813 

//  g)  OH  basic  PROPRGRTION  REfiCTIONS  - 

R[Igll  =  OH  +  CO  - >  H02  +  C02  B  1.50E-13  ♦  (1.0+0.B*Patm)  j 

R[Ig2]  =  OH  +  CH4  - >  X02  +  HCHO  +  H02  0  2.45E-12*EXP(-1775.0AK)  i 

it  NRSfl  97 


it  ORGANIC  CHEMISTRV  «=«===*=== 
it  Cl  =«  formal dhyde  chemistry 


NRHES 

PhotoRaUlDs  +- 


HCH0_toJi02,  HCH0_to_H2  }; 


R[Ci  1] 
RIC1_2] 
RlCl_3] 
RIC1„41 
RIC1_5] 


HCHO 

HCHO 

noD 

-hv-> 

-hv-> 

2.0*HO2  +  CO 

CO  +  H2 
nu  X  uno  X  rn 

HCHO 

Uor 

nu 

wri  “  fiiJX  ^  wu 

X  rn 

HCHO 

HCHO 

+ 

OH 

N03 

HH03  +  H02  +  CO 

a  j lHCHO_to  H021  ; 
a  j  IHCHO^to  J!2]  ; 
a  3.4E“11*EXP(~1600.0AK)  \ 
a  B.6E-12*EXPC20.0AK)  i  /AUPflC,97 
a  2.0E-12*EXP(-2430.0AK)  ! 


//  Cl_5  is  from  IUPflC97  and  is  equal  to  5.BE-16  e29aK 


//  C2  higher  aldehyde  chemistry 


NRMES 

PhotoRotelDs  +=  {  CH3CH0_to_HC0  }j 


REC2^1  «= 
R[C2_2]  « 
R[C2„3J  « 
R[C2^4I  « 


CCHO 

X02  + 

2.0+HO2  +  CO  + 
OH 

CCHO 

rrun 

* 

03P 

nu 

— — >  C203  + 

H20 

HN03 

LLnU 

CCHO 

+ 

Un 

N03 

- >  C203  + 

HCHO  a  j [CH3CH0_to_HC0l  } 
a  i.eE~ii*EXPC“ii00.0AK)  \ 
a  5.6E“12*EXP(  270. 0AK)  J 

a  1.4E-12*EXP(-1060.0AJ<)  J 
//all  NRSR97,  IUPflC97 


//  PAN  PAN  chemistry - - - - 

/♦ 

♦  fill  PRN  Rates  from  NRSfl  1997. 

♦  The  reoommended  rates  for  PRNJL  are  based  on  neu  data  in  NflSfl97 

♦  that  tuas  not  used  in  IUPflC97.  These  are  consistent  with  PRNJZ 

♦  The  reoommended  rates  for  PRN_2f  (see  Note  Table2|  051)  are  the 

♦  same  as  IUPflC97,  but  with  F_c=0.B  instead  of  0.3;  here  we  use 

♦  the  original  ref's  values  which  are  those  of  IUPRC97. 

♦  The  equilibrium  constant  used  in  PflN_2r  is  also  bosed  on  new  data  In  NflSfl97 

♦  that  is  not  cited  in  IUPflC97. 

RfpRN  11  «  C203  +  NO - >  N02  +  X02  +  HCHO  +  H02  B  5.3E-12^EXP(360.8/TK)j  //NflSfl97 

R  PAN  2f]»  C203  +  N02 - >  PAN  8  TROE(2.7E~28*TJ300'‘-7.1, 

1.2E-ll*T„300^-0.9, 

biM],  0.3);  /AUPflCO7«NRSfl97a0.6 

R[PflN_?rl*  PAN - >  C203  +  N02  (1  k [PflN_2f) /(9.0E-29*EXP( 14000. 0AK) ) J  //NASAQ 


NAMES 

IgnSpoIDs  +**  {  CC03H, 


CC02H  };  //  carbonyl  peroxide  and  aoetio  acid 


RIPRN  3]  =  C203  +  H02  - >  0.75*(CCO3H  +  02)  +  0.25*(CCO2H  +  03) 

R[PflN“4l  -  C203  +  C203  - >  2.0*(XO2  +  HCHO  +  H02) 

//  PAN  3  and  PRN  4  from  IUPflC97  ( including  produots  In  PAN_3) 


a  4.3E-13*EXPC1040.0AK); 
a  2.8E-’12*EXPC530.0AK)} 


//  Organic  formulation  from  Gary,  et  al.,  1989  unless  otherwise  specified 
NAMES 


Table  4.1.  Continued. 
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IgnSpoIDs  +»  {  KETENE  }j 

//  PRRRFFIN  CHEMI5TRV  - 

SCHLRRS 

const  sjaarjL  *  -0.11,  Bj3ar_2  «  -2.1j 

R[PflRJ.l  *  OH  +  PAR - >  0.e7*XO2  +  0.i3*XO2N  +  0.11*HO2  +  0.11*CCHO  + 

0.76*ROR  +  ■_papJ.*PRR  Q  8.14E-13; 

R[PflR  2]  *  ROR - >  1.10*CCHO  ♦  0.96*XO2  +  0.94*HO2  +  s^r^*PflR  + 

0.04*XO2N  +  0.02*RaR  8  1.0E+15*EXP(-8000.0/TK); 

R[PflR  3]  «  ROR - >  H02  +  KETENE  8  1.6E+3J 

R|PfiRl43  ■  ROR  +  N02  - >  fiLKN03  8  1.5E-11| 


//  ETHENE  CHEMISTRV  - 

//  Rfii^ised  Ethane  kinetics  (except  ETH_1)  based  on  data  in 
//  fitkinson,  J.Phys.Chem. Ref .Data,  wol.20,no.2,  1997 


RIETHJ.} 

R[ETH_2] 

R[ETH_3] 


03P  +  ETH  - >  0.49*HCHO  ♦  0.60+XO2  +  0.95*CO  + 

1.55*H02  +  0.35*OH  8  1.04E-ll*EXP(-792.0/TK); 

OH  +  ETH  - >  X02  +  1.56+HCHO  +  H02  + 

0.22*CCHO  8  TROEC7.0E-29*TJ300-'-3.1, 

9.0E-12,  btMl,  0.7)5 

03  +  ETH  - >  1.03*HCHO  +  0.325*CO  +  0.08*HO2  + 

0.02*H2O2  +  0.08*OH  8  O.lE-15»EXPC-25a0.0AO J 


//  OLEFIN  CHEMISTRV  - 

//  Olefin  OzonoXysis  and  reaction  with  0  lumping  reevaluated  based  on  data  in 
//  Atkinson,  J.Phys.Chem.Ref .Dato,  vol.26,no.2,  1997 


SCRLRRS 

const  sjole  ■  -1.0  5 


RIOLEJ] 


R[0LEJ2] 


fl[0LE.3j 


R[0LE^4) 


03P  +  OLE 


OH  +  OLE 


03  +  OLE 


N03  +  OLE 


•>  0.49*  CCHO  +  0.29*HO2 

0.20*  HCHO  +  0.007*XO2N 


■>  HCHO  +  CCHO  +  X02  +  H02 


- >  0.52*CCHO 

0.42*HO2 

s_ole*PflR 

0.0595*CH4 

— ->  0.9i*(XO2 
8.09*XO2N 


+0.86*  HCHO 
+  0.45*XO2 
+  0.0e*H2O2 


+  N02)  +  HCHO  + 
+  s_ole*PflR 


+  0.19*XO2  +  0.20+CO  + 

+  0.Ol*PflR  +  0.10*OH 
8  4.0E-12  5 
+  m  ole*PRR 

a  TROEC  8 . 0E-27*T_300^-3 . 5 . 
3.0E-11,  b[Ml,  0.5) 
//  lUPflC,  97 
+  0.3947*CO  + 

+  0.31*OH  + 

+  0.1948*C02  '*■ 


I 


a  5.5E-15*EXPC-18B0.0AK)  ; 
CCHO  + 

8  4.eE-13*EXPC-115S.0/TK)  5 


//  ISOPRENE  CHEMISTRV— CONDENSED - 

//  SRI  version  of  the  one-product  condensed  chemistry  from  Cor ter  (1995) . 


R[IS0_1]  «  03P  +  ISOP 

0.25*HO2  +  0.25*XO2  +  0.75*ISPD  + 

0.50*HCHO  ♦  0.25*PflR  +  0.25*C2O3 

8  3.60E-11  5 

RllSO^l  »  OH  +  ISOP 

> 

0.991*XO2  +  0.629*HCHO  +  0.912*HO2  + 

0.088+XO2N  +  0.912*ISPD 

8  2.54E“11*EXP(  407. 6 AK)  5 

RlISD_3l  «  03  +  ISOP 

0.60*HCHO  +  0.15*CCHO  +  0,35*PRR  + 
e.065*CO  *  0.e66*HO2  *  0.266«OH  + 

0.20*C2O3  +  0.20*XO2  +  0.65*ISPD 

a  7.86E-15*EXP(-1O12.0/TK)  | 

RtIS0_4]  »  N03  +  ISOP 

> 

X02  +  0.65*ISPD  +  0.80*ISPNO3  +  0.00*HO2  + 

0.20*NO2  +  0.80*CCHO  +  2.40*PflR 

a  3.03E-12*EXP(  -448. 0AK)  5 

R[IS0_5I  -  N02  +  ISOP 

0.B0*CCHO  4  2.46*PRR  4  0.80*ISPNO3  4 

X02  4  0.80*HO2  4 

0.20*ISPD  4  0.20*NO  a  1.50E-195 

RIISPD_U  »  ISPD  +  OH 

- > 

1.565*PflR  4  0.167*HCHO  4  0.713*XO2 

4  0.503*HO2  4  0.334*CO 

4  0,158*MOLV  4  0.273*CCHO 

4  0.49a*C2O3  a  3.36E-li5 

R[ISPD_21  «  ISPD  +  03 

-> 

0.114*C2O3  4  0.1504HCHO  4  0.850*MGLV 

4  0.i54*HO2  4  0.268*OH 

4  0.064*XO2  4  0.020*CCHO 

4  0.360*PflR  4  0.225*CO  8  7.IIE-I85 

RlISPD_3l  «  ISPD  +  N03 

- > 

0.357*CCHO  4  0.282*HCHO  4  1.2B2*Pflfl 

4  0.925*HO2  4  0.643*CO 

Table  4.1.  Continued. 
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+  0.850*ISPND3  +  0.075*0203 
+  B.075*XO2  +  0.075+HNO3  8  1.0E-15j 


NAMES 

PhotoRatelDs  +*  {  flCn0_to_R02  )j 


B[ISPDJH  ■  ISPD 

//  TOLUENE  CHEMISTRV  ~ 
RtTOL^l  -  OH  +  TOL 


R[TOL^]  =  T02  +  NO 
R[T0L_3]  «  T02 


-hw->  0.333*00  +  0.067*CCHO  +  0.90B*HCHO 

+  0.832*PRR  +  1.033*HO2 

+  0.700*XO2  +0.967*0203  B  j [flCR0_loJ^02l j 


->  0.0B*(XO2  +  H02)  +  0.36* (ORES  +  H02)  +  0.56*TO2 

0  1.8E“12*EXPC355.0/TK)  j 


->  0.90*(NO2  +  H02  +  OPEN)  +  0.1*flRONO3 
->  ORES  +  H02 


Q  8.1E~i2  ; 
3  4.2  { 


R[CRll  =  OH  +  ORES  • 
RICR21  «  N03  +  ORES  ■ 
R[CR3J  «  CRO  +  N02 


->  0.40*CflO  +  0.60*(XO2  +  H02)  +  0.30*OPEN  3  4.10E-11  ; 


->  CRO  +  HN03 
->  RRONDS 


3  2.20E-li  J 
3  1.4BE-11  j 


SCflLRRS 
k0PEN_R  i 

R[RF_ll  ' 
R[RF_2]  ^ 


6.0J  //  Soal9d  back  from  Gery's  EPR  rpt«9.04xjHOHOr  by  Baas/Moor«0 . 70 

3  jlHCH0„toJ<02l*k0PENJ^  j 
3  3.0E-11  ; 


OPEN  -hw->  C203  +  H02  +  CO 

OPEN  +  OH  - >  X02  +  2.8*00  +  2.0*HO2  +  0203  +  HCHO 


R[RF  3]  »  OPEN  +  03 - >  0.03*CCHO  +  0.62*0203  +  0.70*HCHO  + 

0.03+XO2  +  0,69*00  +  0.08*OH  ♦ 

0.76*HO2  +  0.20+MDLV 


3  5.40E“17*EXP(~500.0/T1C)  j 


//  XVLENE  CHEMISTRV 
RtXVL^l]  «»  OH  +  XVL 


->  0.10*(/*  XL02«  */  H02  +  X02  +  PRR)  + 

0.20*(CRES+HO2+PflR)  + 

0.30*TO2  + 

0.40*(/*  XINT=  */  MOLV  +  MGLV  +  PAR  +  PRR  +  H02) 

3  1.70E-11*EXP(11B.0AK)  J 


SCflLRRS 
kMGLVJ^  = 

R[XVL_2] 

R[XVL_3l 


6.01  //Scaled  back  from  Gary's  EPR  rpt  9,64xJ  HCHO  r  by  Bass/Moor “0.78 


MOLV  -hv->  C203  +  H02  +  CO 
OH  +  MGLV  - >  X02  +  C203 


3  jlHCH0_toJ102l*kMGLV_R  j 
8  1.70E-11  i 


//  NOTE  WELL:::  MGLV  here  is  ♦not*  msant  to  ba  real  MGLV,  but  a  ring 
//  ring  fragmentation  product  of  XVL  luhioh  has  been  tuned  based 

//  upon  original  Bass  HCHOr  rates  and  assumed  yields.  Most 

//  importantly  the  sigma-phi (MGLV)  of  lUPRC  can  not  be  used 

//  to  compute  the  j (MOLV) . 


r[mnjlT 

l1eN02 - > 

H02  +  HCHO  +  NO 

8  J [HONO] 

R[flF^l 

“  OH  + 

MeOH - > 

HCHO  +  H02  +  H20 

3  6 . 7E-12*EXP (-600 . 0 AK) } 

RlflF_2l 

*  OH  + 

EtOH - > 

0.945*CCHO  +  H02  + 
0.055*XO2  +  0.11*HCHO 

8  7.0E-12*EXP(-235.0AK)i 

•end 

//  Operator  Chemistry  - 

RIX0J.1  =  X02  +  NO 

RlX0_2j  *  X02N  +  NO 

RtXO  31  -  X02  +  X02 

fl[X0  41  =  X02  +  H02 

RIX0“51  =  X02N  +  H02 

R[X0l6l  -  XD2N  +  X02N 

R[X0  71  “  X02  +  X02N 


- >  N02 

- >  XRN03 

> 


8  3.0E-12*EXP(280.0AK)|  //NflSft97,Ti 
3  kIXOJ.1  I 

3  2.5E-13*EXP(190.0AK)i  //  1.70E-14*EXP(1300 
8  3.8E-13*EXP(800.0AK)}  //  7,68E-14*EXP(1300 
8  krX0_4l  ; 

8  kiXOJBl  J 
8  ktXO„3l*2.0| 


//  Declare  species  NR,  which  is  used  to  prouide 
//  mass  balance  with  ambient  inputs  of  propane,  eto. 

RINR J.1  =  - >  NR  3  ( 

//  end  of  principle  mechanism  file. 


Table  4.1.  Continued. 
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Genus 


Grid  Cell  Area  in  Cell 

I  i 

NO  Emission  Factor 
(|xg  m-2  h-') 

41 7.385,  27  M - Total  hectares  in  the  following  27  categories 


•Acer’ 

12.84 

4.5 

Acer  (maple) 

’Alnu^ 

0.00 

4.5 

Alnus  (European  alder) 

'Betu’ 

1.30 

4.5 

Betula  (birch) 

•Carp’ 

0.61 

4.5 

Carpinus  (hombean) 

•Cary’ 

26.67 

4.5 

Carya  (hickory) 

'Cere’ 

0.39 

4.5 

Cercis  (redbud) 

'Coru’ 

11.72 

4.5 

Cornus  (dogwood) 

'Dios’ 

1.22 

4.5 

Diospyros  (persimmon) 

'Fagu’ 

2.33 

4.5 

Fagus  (american  beech) 

'Frax’ 

1.29 

4.5 

Fraxinus  (ash) 

'Ilex’ 

1.91 

4.5 

Llex  (holly) 

•Jugl’ 

1.51 

4.5 

Juglans  (black  walnut) 

'Juni’ 

2.53 

4.5 

Juniperus  (east.  Red  cedar) 

'Liqu’ 

8.54 

4.5 

Liquidambar  (sweetgum) 

’Liri’ 

11.30 

4.5 

Liriodendron  (yellow  poplar) 

’Magn’ 

0.22 

4.5 

Magnolia 

•Nyss’ 

4.87 

4.5 

Nyssa  (blackgum) 

'Ofor’ 

83.90 

’Oxyd’ 

23.29 

4.5 

Oxydendrum 

’Pinu’ 

95.59 

4.5 

Pinus  (pine) 

'Prun’ 

0.89 

4.5 

Prunus  (cherry) 

'Quer’ 

120.28 

4.5 

Quercus  (oak) 

'Robi’ 

0.87 

4.5 

Robinia  (black  locust) 

•Sali’ 

0.19 

4.5 

Salix  (willow) 

'Sass’ 

0.16 

4.5 

Sassafras 

'Tsug’ 

2.42 

4.5 

Tsuga  (Eastern  hemlock) 

'Ulmu’ 

0.41 

4.5 

Ulmus  (American  elm) 

0.000,  0 

766.468, 

7  ◄ - 

-  Total  hectares  in  the  following  7  crop  types 

’Corn’ 

6.90 

577.6 

Corn 

'Cotf 

0.71 

256.7 

Cotton 

'Hay  ’ 

202.39 

12 

Hay 

’Msep’ 

539.48 

12.8 

Misc  Crops 

'Rye  ’ 

0.60 

12.8 

Rye 

'Soyb’ 

5.25 

12.8 

Soybean 

'Whea’ 

11.11 

192.5 

Wheat 

416.147, 

4  ◄ - 

—  Total  hectares  in  the  following  4  grass  types 

'Barr’ 

92.66 

0 

Barren 

'Gras’ 

1.85 

57.8 

Grass 

'Othe’ 

321.34 

57.8 

Other  (unknown,  asuume  grass) 

’Scru’ 

0.28 

57.8 

Scrub 

Table  4.2.  Example  of  land-use  statistics  that  are  available  for  each  4  x  4  km  grid  cell 
in  the  model  domain.  This  data  is  then  used  in  calculating  the  biogenic  NO  budget. 
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North  Carolina  1996  MAQSIP  Domain  (9/99) 
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Figure  4.1.  Model  domain  used  for  this  study  consists  of  a  480  km  x  276  km  area  covering  most  of  North  Carolina.  Grid  scale 
resolution  is  4  x  4  km. 
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Figure  4.2.  Schematic  of  the  Sparse  Matrix  Operator  Kernel  Emissions  (SMOKE)  modeling  system.  Source:  MCNC, 
North  Carolina  Supercomputing  Center,  Environmental  Programs. 
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Figure  4.3.  Schematic  of  the  Multiscale  Air  Quality  Simulation  Platform  (MAQSIP).  Source;  MCNC,  North  Carolina 
Supercomputing  Center,  Environmental  Programs. 


□  BEIS  Estimates 
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Figure  4.4.  Estimated  NO  flux  at  field  site  sampled  throughout  1999-2000.  Clear  bar  represents  estimates  using  existing 
BEIS2  model.  Stippled  bar  represents  estimated  NO  flux  using  algorithm  developed  during  this  study. 
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Biosolids  applied  on  between  2000-3000  hectares 
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Figure  4.6.  NO 
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Figure  4.7.  Cells  in  the  model  domain  receiving  all  of  the  biosolids  which  are  applied  within  the  respective  coimties.  The 
scale  indicates  how  much  the  biogenic  NO  emissions  increase  over  the  base  case  scenario. 
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Figure  4.9.  North  Carolina  NO^  utility  sources  for  1995. 
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Figure  4.10. 
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Figure  4.11.  Hourly  averaged  NO  emissions  (moles/hour)  of  all  the  4x4  km  grid  cells  in  the  model  domain. 
Anthropogenic  NO  is  the  solid  line  (primary  axis)  and  the  biogenics  are  represented  by  the  dashed  line  (secondary  axis). 
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Figure  4.12.  Weekly  profile  of  rural  and  urban  mobile  sources.  Source:  MCNC,  2001. 
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Figure  4.13b.  Ratio  of  biogenic  NO  to  anthropogenic  NO  in  Granville  County,  NC  (weekend,  early  morning). 


re  4.14.  Change  in  O3  concentrations  throughout  Wednesday,  June  26  and  Thursday,  June  27, 1996  (Modified  case 
case)  plotted  at  8:00  pm  (a),  6:00  am  (b),  and  8:00  am  (c). 
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Figure  4.15.  Change  in  O3  concentrations  throughout  Saturday,  June  29  and  Sunday,  June  30, 1996  (Modified  case  -  Base 
case)  plotted  at  10:00  am  (a)  and  4:00  am  (b). 


Figure  4.17.  Percent  change  of  ozone  in  the  modeling  domain.  Time  series  plot  represents  the  average  percent  change 
throughout  the  24  hour  period  (8:00  am  Saturday  -  8:00  am  Sunday).  Negative  values  indicate  that  ozone  levels  are  lower 
in  the  modified  case  than  in  the  base  case. 
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Figure  4.18.  Percent  change  of  ozone  in  the  modeling  domain  with  biosolids  concentrated  within  the  counties.  Time  series 
plot  represents  the  average  percent  change  throughout  the  24  hour  period  (8:00  am  Saturday  —  8:00  am  Sunday).  Negative 
values  indicate  that  ozone  levels  are  lower  in  the  modified  case  than  in  the  base  case. 
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Figure  4.20.  Ratio  of  NO^  to  NO^  in  MAQSIP  modeling  domain. 


Chapter  V.  Characterization  of  Ammonia  Emissions  from  Soils 
IN  THE  Upper  Coastal  Plain,  North  Carolina 

by 

Paul  A.  Roelle  and  Viney  P.  Aneja 
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Abstract 

A  dynamic  flow-through  chamber  system  was  used  to  measure  fluxes  of 
ammonia-nitrogen  (NH3-N,  where  NH3-N=(14/17)*NH3)  from  soil  surfaces.  The 
research  site  was  located  in  eastern  North  Carolina  (35.9°N  Latitude;  77.7°W  Longitude) 
and  measurements  were  conducted  during  spring  and  winter  2000,  in  order  to  assess  the 
NH3  source  strength  of  intensively  managed  agricultural  soils  and  the  physiochemical 
properties  which  control  these  emissions.  Soil  temperature  (Tsoii),  soil  pH,  soil  moisture. 
Total  Kjeldahl  Nitrogen  (TKN=organic  N+NH3-N-^NH4'^-N)  were  monitored  throughout 
both  research  periods.  Soil  temperature  was  fovmd  to  explain  the  largest  variability  in 
soil  NH3  emissions  [Log,oNH3-N  Flux  (ng  N  m‘^  s-')=  0.054*Tsoii  (°C)  +  0.66;  R^=0.71), 
suggesting  that  an  approach  similar  in  design  to  the  Biogenic  Emissions  Inventory 
System  (BEIS2)  land  use  and  temperature  model  for  NO  emissions,  might  be  effective 
for  modeling  biogenic  NH3  emissions.  Soil  nitrogen  was  also  significant  in  predicting 
NH3  flux  [NH3  Flux  =55.5*(NH3-N)-160,  rM.86;  NH3  Flux=0.6*(TKN)-410,  R^=0.27], 
but  only  after  the  two  days  with  the  heaviest  rainfall  were  removed  from  the  regression. 
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emphasizing  the  role  of  soil  moisture  in  controlling  the  transfer  of  gases  across  the 
soil/air  interface.  Soil  pH  remained  relatively  constant  throughout  both  research  periods 
and  therefore  did  not  serve  as  a  useful  predictor  of  NH3  flux.  A  rain  event,  followed  by  a 
drying  period  produced  a  characteristic  pulse  in  ammonia  emissions.  This  pulsing 
phenomena  has  been  observed  for  other  trace  gases  by  various  researchers.  This  research 
location  was  the  site  of  a  commercial  hog  operation,  which  allowed  for  the  comparison  of 
soil  and  lagoon  emissions  (lagoon  emissions  were  based  on  an  algorithm  developed  by 
Aneja  et  al.,  2000).  An  analysis  of  the  source  strengths  confirmed  that  lagoon  emissions 
are  a  larger  flux  source  (average  lagoon  flux  ~  1 8,137  ng  N  m'^  s'';  average  soil  flux  ~  54 
ng  N  m'^  s''),  however  soil  surfaces  make  up  a  larger  fraction  of  a  commercial  hog 
operation  than  the  lagoon  surfaces,  and  as  a  result  they  can  not  be  neglected  when 
developing  and  apportioning  NH3  emissions.  A  yearly  average  of  ammonia  emissions  at 
this  site  revealed  that  soil  emissions  represent  approximately  28%  of  the  lagoon 
emissions. 
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Introduction 


Ammonia  (NH3)  is  an  important  atmospheric  trace  species,  both  in  terms  of  its 
effect  on  tropospheric  chemistry  and  due  to  its  impact  on  ecosystems.  Ammonia,  which 
is  the  most  abundant  alkaline  specie  in  the  atmosphere,  is  critical  to  neutralizing  acids 
formed  through  the  oxidation  of  sulfur  dioxide  (SO2)  and  nitrogen  oxides  (NOx)  (Asman 
et  al.,  1 998).  When  NH3  is  deposited  onto  the  soil,  it  is  both  taken  up  by  plants  and 
converted  by  bacteria  into  nitrate  (NO3')  (nitrification)  (Lekkerkerk  et  al.,  1995).  The 
nitrification  process  forms  hydrogen  ions  leading  to  acidification  of  the  soil,  leaching  of 
N03‘  to  groundwater,  and  possible  deficiencies  of  other  plant  nutrients  such  as  potassium 
(K"^)  and  magnesium  (Mg^'^)  (Asman  et  al.,  1998).  In  addition  to  acidification,  excess 
nitrogen  loading  can  lead  to  over  enrichment  of  both  land  and  water  ecosystems.  Further, 
excess  N  deposition  can  cause  the  above  ground  portion  of  the  plant  to  grow  rapidly, 
leaving  the  root  system  relatively  smaller  and  weaker  and  more  susceptible  to  disease  and 
harsh  weather  conditions  (Lekkerker  et  al.,  1995). 

Unlike  oxidized  nitrogen  and  sulfur  compounds  (NOx  and  SOx  respectively), 
which  are  predominately  emitted  from  industrial  processes,  NH3  is  primarily  emitted  by 
agricultural  sources  and  therefore  requires  different  control  strategies  (Sutton  et  al., 

1993).  A  review  of  the  current  literature  revealed  that  soil  and  plant  emissions  account 
for  anywhere  between  11%  and  28%  of  the  global  NH3  budget  (Table  1.2).  In  North 
Carolina,  which  is  currently  the  second  largest  pork  producing  state  in  the  US,  the 
percentage  of  ammonia  emissions  from  swine  is  approximately  46%  compared  to  the  US 
average  of  10%  (Battye  et  al.,  1994)  (Figure  1.6  a,b). 
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Accurate  inventories  of  NH3  are  needed  to  model  both  its  transport  and  deposition 
(Misselbrook  et  al.,  2000).  Beyond  quantifying  this  overall  source  strength  and  the 
effects  to  ecosystems,  there  needs  to  be  an  accurate  budget  and  reliable  source 
apportionment  of  the  various  NH3  emission  pathways  from  these  intensively  managed 
animal  operations  (Figure  5.1).  Current  inventories  that  are  then  used  in  air  quality 
models  are  determined  using  emission  factors  (ie  tons  NH3/animal/year  or  %  NH3-N 
emitted/tons  of  N-fertilized  applied),  most  of  which  are  based  on  European  studies. 
Therefore,  the  objectives  of  this  study  are  to  quantify  and  relate  the  soil  NH3  flux  at  an 
eastern  North  Carolina  site  to  different  environmental  variables  in  the  soil,  and  determine 
the  source  strength  of  these  soil  emissions.  Moreover,  the  soil  ammonia  emissions  will 
be  compared  and  contrasted  to  ammonia  emission  measurements  made  over  animal 
(swine)  waste  treatment  and  storage  lagoons  (Aneja  et  al.,  2000). 

Methods  and  Materials 
Sampling  Site  and  Sampling  Scheme 

The  NH3  flux  measurements  were  made  at  the  Upper  Coastal  Plain  Research 
Station  (Figure  5.2),  located  in  Edgecombe  County,  North  Carolina  (See  Table  5.1  for 
site/soil  characteristics).  This  facility  is  operated  with  typical  agronomic  and  husbandry 
practices  for  the  respective  crops  and  animals  and  contains  approximately  178  hectares, 
101  of  which  are  cropland  soils.  The  facility  also  maintains  a  farrow-to-finish  hog 
operation  with  approximately  1250  hogs  on  site.  The  waste  from  the  animals  (urine  and 
feces)  is  flushed  from  the  hog  production  houses  into  two  uncovered  anaerobic  waste 
treatment  and  storage  lagoons  (a  primary  and  secondary  lagoon,  total  acreage  ~  1 
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hectare).  The  effluent  from  these  lagoons  is  periodically  sprayed  to  the  crops  as  a 
nutrient  source.  A  corn  crop  was  planted  on  April  5*'^,  2000  and  spring  measurements 
began  on  April  26*  and  ended  on  May  14*,  2000.  Due  to  several  rain  events  and 
equipment  malfunctions,  only  5  days  were  suitable  for  conducting  the  measurements,  all 
of  which  are  presented  in  Table  5.1. 

Although  the  field  sampled  is  typically  sprayed  with  lagoon  effluent,  it  was  not 
sprayed  during  the  2000  measurement  campaign  and  instead  received  approximately  146 
kg  N/ha  on  May  24*,  2000.  The  corn  crop  was  harvested  on  August  21**  and  the  stalks 
were  shredded  and  left  on  the  soil  surface.  No  cover  crop  was  planted  and  the  winter 
measurements  were  conducted  from  Dec  13*’’- 19*'*,  2000.  Rain  events  prevented 
measurements  on  two  of  the  days,  although  measurements  were  made  on  the  remaining  5, 
all  of  which  are  presented  in  Table  5.1. 

NH3  concentration  measurements  were  made  on  1 0  random  sampling  plots 
located  within  a  10  m  radius  of  a  mobile  laboratory  (Modified  Ford  Aerostar  Van, 
temperature  controlled  to  within  the  operating  range  of  the  instruments).  The  sampling 
scheme  consisted  of  measuring  concentrations  of  NH3  after  the  sample  exited  the 
dynamic  flow-through  chamber  system.  The  system  recorded  60-second,  rolling 
averages  of  NH3  concentrations.  These  values  were  then  binned  and  averaged  every  15 
minutes.  The  15-minute  binned  averages  were  used  in  all  flux  calculations.  N  in  Table 
5.1  refers  to  the  number  of  these  binned  averages  (to  also  include  soil  temperature  as 
discussed  in  Soil  Analysis). 
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A  daily  experiment  consisted  of  placing  the  chamber  on  the  stainless  steel  collar, 
which  had  been  inserted  into  the  soil  the  previous  evening.  The  collars  were  all  located 
on  bare  soil  with  no  plants  being  enclosed  within  the  collar  or  chamber  system.  The 
chamber  was  placed  on  the  collar  at  approximately  8:00  AM  and  flushed  with  zero  grade 
air  for  at  least  one  hour  before  data  collection  began  at  9:00  AM.  This  sampling  scheme 
ensured  that  the  concentrations  within  the  chamber  reached  steady  state  prior  to  any  data 
acquisition  and  allowed  for  the  instruments  to  undergo  their  daily  calibrations.  Daily 
experiments  ended  at  approximately  5:00  PM  and  the  stainless  steel  collar  was  relocated 
to  a  random  location  within  a  1 0  m  radius  of  the  mobile  laboratory,  in  preparation  for  the 
next  days  experiment.  This  procedure  allowed  a  minimum  of  16  hours  for  any  effect  on 
soil  NH3  flux,  due  to  soil  disturbances  caused  by  the  insertion  of  the  stainless  steel  collar, 
to  dissipate. 

Instrumentation  and  Flux  Calculation 

The  chamber  design,  associated  mass  balance  equation  and  calibration  procedures 
are  described  in  full  in  Chapter  I. 

Soil  Analysis 

A  soil  sample  was  taken  from  the  center  of  the  chamber  footprint  at  the  end  of 
each  measurement  period  (1  sample  per  measurement  period),  and  analyzed  for  soil  pH, 
soil  moisture  and  Total  Kjeldahl  Nitrogen  (TKN=organic  N+NH3-N+NH4'*’-N)  by  the 
North  Carolina  State  University  Department  of  Biological  and  Agricultural  Engineering. 
Percent  water  filled  pore  space  (%  WFPS)  is  a  measure  of  soil  water  content  and  can  be 
expressed  as  the  percentage  of  pore  spaces  in  the  soil  that  are  filled  with  water.  The 
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%WFPS  is  a  convenient  expression  to  describe  soil  moisture  because  it  accounts  for  the 
differing  bulk  and  particle  densities  of  soils  and  therefore  allows  for  the  comparison  of 
soil  moisture  from  different  soil  types.  Soil  temperature  was  measured  with  a  Campbell 
Scientific  temperature  probe  (accuracy  ±  3%)  inserted  into  the  soil  to  a  depth  of 
approximately  5  cm.  Air  temperatures  (Campbell  Scientific;  accuracy  ±  3%)  were 
measured  inside  of  a  radiation  shield  at  a  height  of  1.5  meters.  Soil  data  was  stored  in  15 
minute  binned  averages  utilizing  a  Campbell  Scientifie  2 IX  Mierologger. 

Results  and  Discussion 
Environmental  Controls  on  NH3  Flax 

The  ammonia-water  and  ammonia-soil  system  have  been  studied  in  the  past 
because  of  their  industrial  importance  and  as  a  means  for  studying  the 
absorption/desorption  mechanism  (Whitman  and  Davis,  1924;  Godfrey,  1973;  Levenspiel 
and  Godfrey,  1974;  Ibusuki  and  Aneja,  1984;  Leuning  et  al.,  1984;  Warneck,  2000). 
These  previous  studies  indicate  that  the  three  most  important  parameters  in  the  biological 
and  chemical  processes  which  determine  the  NH3  equilibrium  and  production  rate  are 
temperature,  pH  and  nitrogen  content  of  the  soil.  In  the  soil  environment  these  variables 
translate  to  soil  temperature,  soil  pH,  and  the  soil  moisture  content. 

Soil  Temperature 

The  temperature  dependence  is  linked  to  the  NH3  production/emission  through 
chemical  and  biological  processes  occurring  simultaneously  in  the  soil  environment.  In 
the  absenee  of  recent  N  fertilization,  ammonia  appears  in  the  soil  through  a  process  called 
mineralization  or  ammonification  whereby  microorganisms  satisfy  their  energy  needs  in 
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the  soil  by  converting  amino  acids  in  the  dead  organic  matter  to  NH3  as  in  the  following 
reaction  (Delwiche,  1970;  Kinzig  and  Socolow,  1994): 

CH2NH2COOH  +  1  1/202->2C02  +  H2O  +  NH3  (1) 

Given  that  biochemical  reactions  have  been  shown  to  rise  exponentially  with  temperature 
in  the  range  between  288  -  308  K,  one  would  expect  there  to  be  a  corresponding  increase 
in  the  soil  NH3  concentration  as  soil  temperature  increases  (assuming  that  the  system  is 
not  limited  by  soil  organic  matter  content)  (Wameck,  2000).  In  fact,  this  exponential 
dependence  of  both  reduced  and  oxidized  nitrogen  trace  species  on  soil  temperature  has 
been  repeatedly  demonstrated  in  both  field  and  laboratory  studies  (Sherlock  and  Goh, 
1985;  Van  der  Molen  et  al.,  1990;  Roelle  et  ak,  1999). 

The  NH3  and  NH4'^  in  the  soil  solution  are  in  aqueous  equilibrium, 

NH/(aq)^NH3(aq)  (2) 

and  if  the  NH3  vapor  pressure  in  solution  is  greater  than  the  vapor  pressure  of  NH3  in  the 
surrounding  air  than  the  NH3  will  be  volatilized  (Sherlock  and  Goh,  1985).  This 
relationship  can  be  expressed  in  the  form  of  Henry’s  Law,  where  the  dimensionless 
Henry  coefficient  (H)  can  be  written  as: 

H=(NH3(aq))/(NH3(g))and  (3) 

Log  H- -1.69 +  1477.7/T  (4) 

(Hales  and  Drewes,  1979;  Sherlock  and  Goh,  1985). 

From  equations  3  and  4  above,  it  can  be  shown  that  an  increase  in  soil 
temperature  will  produce  a  corresponding  increase  of  the  NH3(g)  concentrations  in  the 
soil.  This  same  Henry’s  law  equilibrium  applies  in  water  bodies  (such  as  hog  waste 
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lagoons)  and  Aneja  et  al.  (2001a)  have  developed  a  fundamental  mechanistic  ammonia 
model  to  predict  lagoon  NH3  emissions  based  on  this  equilibrium.  The  Aneja  et  al. 
(2000)  study  confirmed  that  lagoon  temperature  was  the  dominant  parameter  affecting 
NH3  emissions  from  lagoons  and  the  temperature  dependence  model  that  they  developed 
was  then  compared  to  a  temperature  dependence  model  developed  in  this  study  for  NH3 
emissions  from  soil  surfaces.  Other  physiochemical  parameters  that  influence  NH3 
emissions  are  pH,  nitrogen  content,  and  moisture  content  of  the  soil. 

The  calculated  NH3  fluxes  (including  rainfall  totals  which  are  discussed  later) 
during  the  spring  and  winter  2000  measurement  period  can  be  seen  in  Figure  5.3.  The 
values  calculated  during  this  study  fall  within  the  range  of  NH3  soil  fluxes  reported  by 
other  researchers  (Table  5.2).  The  daily  averaged  NH3  flux  values  plotted  versus  the 
daily  averaged  soil  temperature  (Figure  5.4)  reveal  a  statistically  significant  relationship 
(p<0.01)  with  NH3  flux  increasing  exponentially  as  soil  temperature  increases.  The  next 
best  fit  to  this  data  was  a  linear  model  which  resulted  in  an  R^=0.62  [ISfH3-N  Flux  (ng  N 
m'^  s’')  =  6.9*  Tsoii  (°C)-35.6].  The  soil  temperature  model  in  this  study  [LogioNH3-N 
Flux  (ng  N  m’^  s’')  =  0.054*Tsoii  +  0.66;  R^=0.71]  and  the  lagoon  temperature  model  in 
the  Aneja  et  al.  (2000)  study  [Log,oNH3-N  Flux  (pg  N  m’^  min')  =  0.048*Tiagoon  +  2.1; 
R^=0.78]  indicate  approximately  the  same  amount  of  variability  in  the  NH3  flux  from  the 
soil  and  lagoon  surfaces.  This  strong  exponential  dependence  of  NH3  emissions  on 
temperature  has  been  cited  for  other  nitrogen  (N)  trace  gases  (NO,  NO2  and  N2O)  with 
approximately  the  same  results  (Kim  et  al.,  1994;  Sullivan  et  al.,  1996;  Roelle  et  al., 
1999).  In  fact,  the  US  Environmental  Protection  Agency  (EPA),  currently  utilizes  this 
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exponential  temperature  dependence  in  the  Biogenic  Emissions  Inventory  System  (BEIS) 
model  to  estimate  the  biogenic  nitric  oxide  emissions  (NO)  which  are  then  used  as  input 
data  for  ozone  and  air  quality  models  (Birth  and  Geron,  1995;  Thornton  et  ah,  1997). 

The  strong  temperature  dependence  reported  in  this  study  suggests  that  a  similar 
approach  may  also  be  effective  in  estimating  the  biogenic  NH3  emissions. 

Soil  pH,  Soil  Moisture  and  Soil  Nitrogen  Content 

As  previously  discussed,  other  parameters  such  as  soil  pH  and  soil  moisture  have 
been  identified  as  controlling  NH3  production.  The  [OH']  produced  as  a  result  of  the 
following  dissociation  in  the  soil  solution 

NH3(aq)ONH4\aq)  +  0H-  (5) 

can  be  represented  by: 

[OH']  =  K\v/[H'^],  (Kw^water  dissociation  constant)  (6) 

As  the  soil  pH  increases  ([OH']  increases),  the  equilibrium  is  shifted  towards  more  NH3 
being  released  (Warneck,  2000;  Li,  2000),  and  studies  have  demonstrated  pH  to  effect 
NH3  emissions  (Singh  and  Nye,  1988;  Aneja  et  al.,  2001a).  However,  in  intensively 
managed  agricultural  soils,  the  pH  value  of  a  soil  column  (~20  cm)  tends  to  remain  fairly 
uniform  (see  Table  5.1),  and  therefore  no  significant  relationships  between  soil  pH  and 
NH3  flux  can  be  discerned.  The  studies  which  identified  soil  pH  as  a  significant  variable 
typically  looked  at  pH  in  the  top  1-3  cm  of  the  soil  and  were  during  fertilization  events, 
both  of  which  found  much  larger  changes  in  pH  values  (1-3  pH  units).  Therefore  it  is 
recommended  that  future  field  studies  should  also  include  an  analysis  of  soil  pH  over  a 
few  different  sampling  depths. 
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The  soil  moisture  conditions  at  this  field  site  are  best  described  by  Figure  5.3, 
which  shows  the  rain  events  in  relation  to  the  sampling  days  and  measured  flux  values. 
The  percent  water  filled  pore  space  (%WFPS)  (Table  5.1)  also  describes  moisture 
conditions,  however  the  effects  of  our  soil  sampling  technique  (20  cm  depth),  has  the 
potential  to  dilute  the  actual  moisture  content  in  the  top  few  centimeters  of  the  soil 
column,  where  the  largest  concentration  of  ammoniacal  nitrogen  has  been  shown  to 
reside  (Singh  and  Nye,  1986).  As  more  water  is  introduced,  pores  in  the  soil  matrix 
begin  to  fill  and  hinder  the  diffusion  of  NH3  gas  from  the  soil  to  the  air  (Kirk  and  Nye, 
1991).  Additionally,  from  an  equilibrium  standpoint,  as  the  water  content  in  the  soil 
increases,  the  NH3  decreases  as  the  equilibrium  moves  towards  the  right  hand  side  of 
equation  (5). 

Increases  in  emissions  have  been  observed  when  soils  with  high  moisture  content 
are  subjected  to  drying  via  high  winds  or  temperature.  This  increase  or  “pulse”  is 
believed  to  be  caused  by  the  combination  of  an  increase  in  the  ammoniacal  nitrogen 
concentration  of  the  soil  and  greater  diffusion  through  the  relatively  drier  soil  (Burch  and 
Fox,  1989;  Battye  et  al,  1994).  A  similar  “pulse”  in  emissions  (day  3  in  Figure  5.3)  may 
be  causing  the  large  increase  in  emissions  which  occurred  after  the  sharp  decrease  in  soil 
moisture  (decreased  from  40.5  %WFPS  to  24.9  %WFPS).  While  this  observation  is 
based  on  only  one  data  point,  prior  to  which  there  was  a  lapse  of  several  days  of  data,  it  is 
supported  by  similar  observations  from  other  researchers  (Burch  and  Fox,  1989;  Battye  et 
al.,  1994). 
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The  relationship  of  NH3  volatilization  and  N  content  of  the  soil  can  be  seen  in 
Figure  5.5(a,b).  Several  investigators  have  attempted  to  develop  models  describing  the 
physical  and  chemical  processes  affecting  the  release  of  NH3  from  soil  surfaces  (Sherlock 
and  Goh,  1985;  Singh  and  Nye,  1986).  Based  on  their  mechanistic  models,  a  linear 
dependence  of  NH3  flux  on  soil  nitrogen  content  was  expected  and  found  to  explain  more 
variability  in  the  data  than  other  attempted  relationships.  In  both  plots  (a  and  b),  there 
was  a  relatively  weak  dependence  of  NH3  flux  on  NH3-N  and  TKN  content  of  the  soil 
(R^=0. 12  and  0.02  for  NH3-N  and  TKN  respectively)  when  all  data  points  are  considered 
in  the  regression.  Given  the  strong  influence  that  soil  moisture  has  on  NH3,  both  in  terms 
of  its  equilibrium  and  in  its  control  of  diffusion  through  the  soil,  both  a  and  b  were 
reanalyzed  taking  the  moisture  conditions  into  account.  When  2  sampling  days  (rain 
events  >  0.3  cm  and  standing  water  evident  in  some  parts  of  the  field)  were  removed, 
there  was  an  appreciable  increase  in  the  significance  of  N  content  on  NH3  release.  The 
rain  events  which  occurred  on  December  13*'^  and  16*'’  were  both  light  mists  with  the 
majority  of  the  total  rainfall  occuring  after  the  sampling  period.  Interestingly,  the  largest 
rainfall  event  (December  17*'’)  did  not  suppress  NH3  emissions  on  December  18*'’,  relative 
to  the  other  wintertime  measurements.  However  strong  winds  associated  with  the  frontal 
passage  following  the  rainfall  event  did  act  to  dry  out  the  surface  layer  even  though  the 
soil  core  had  the  highest  %  WFPS.  The  greater  dependence  of  NH3  volatilization  rates  on 
NH3  content  (R^=0.86)  as  compared  to  TKN  (R^=0.27)  is  expected  as  TKN  is  the  sum  of 
both  ammoniacal  N  and  organic  N. 
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Budget  for  Site 

Aneja  et  al.  (2000)  described  a  temperature  based  model  similar  to  the  model 
described  in  this  study,  although  their  model  related  NH3  emissions  from  hog  waste 
lagoons  to  lagoon  temperature.  While  the  processes  regulating  NH3  emissions  from 
lagoon  surfaces  differ  from  those  regulating  soil  NH3  emissions,  the  temperature 
relationship  is  found  to  be  applicable  to  both.  The  physical  and  chemical  processes 
regulating  the  emissions  from  the  lagoon  and  soil  surfaces  and  their  dependence  on  soil 
temperature  are  described  in  detail  in  Aneja  et  al.  (2001a)  and  Wameck  (2000), 
respectively.  Utilizing  the  temperature  model  developed  in  this  study  to  estimate  NH3 
emissions  from  soils  and  the  temperature  based  algorithm  developed  by  Aneja  et  al. 
(2000)  to  estimate  NH3  emissions  from  lagoons  (See  Table  5.2  for  measured  lagoon  NH3 
emission  averages),  it  is  possible  to  estimate  the  relative  seasonal  source  strengths  of  the 
soil  and  the  lagoon. 

Although  this  yearly  estimate  is  based  on  a  model  developed  from  two  seasons  of 
data,  the  temperature  dependence  is  assumed  to  remain  consistent,  during  periods  when 
the  soil  has  not  been  recently  fertilized,  throughout  diurnal  cycles  and  throughout  the  year 
(Van  der  Molen  et  al.,  1990;  Aneja  et  al.,  2000;  Warneck,  2000).  Using  the  daily 
averaged  temperatures  (Information  obtained  from  North  Carolina  State  Climate  Office), 
the  seasonally  averaged  emissions  from  the  101  hectares  of  soil  surface  at  this  site  during 
spring,  summer,  fall  and  winter  were  determined  to  be  324  kg  NH3-N,  933  kg  NH3-N, 

383  kg  NH3-N  and  83  kg  NH3-N  respectively  (Figure  5.6).  Similarly,  the  seasonally 
averaged  lagoon  emissions  from  the  1  hectare  of  lagoon  surface  during  spring,  summer. 
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fall  and  winter  were  determined  to  be  1,140  kg  NH3-N  2,953  kg  NH3-N,  1,315  kg  NH3-N 
and  340  kg  NH3-N  respectively.  Therefore,  the  NH3  emissions  from  soil  surfaces 
represent  approximately  28%,  32%,  29%  and  24%  of  the  spring,  summer,  fall  and  winter 
lagoon  emissions.  The  typical  practice  in  most  hog  operations  is  to  spray  the  fields  with 
the  hog  waste  effluent  instead  of  fertilizing  them  with  commercially  derived  fertilizers. 
This  field,  however,  was  not  sprayed  and  therefore  this  budget  may  be  biased  low,  as 
emission  factors  for  land  spreading  of  slurry  are  often  cited  as  being  larger  (-15-76%) 
than  the  factors  for  the  commercially  derived  fertilizers  (Misselbrook  et  ah,  2000). 
Conclusions  and  Recommendations 

Utilizing  a  dynamic  flow-through  chamber  technique,  NH3  flux  values  were 
calculated  for  the  spring  and  winter  (2000)  at  an  upper  coastal  plain  site  in  North 
Carolina.  Soil  pH  remained  relatively  constant  throughout  the  measurement  period  and 
therefore  was  not  useful  as  a  NH3  flux  predictor.  The  NH3  flux  values  were  most 
strongly  correlated  with  soil  temperature  [LogioNH3-N  Flux  (ng  N  m'^  s"')  ==  0.054*Tsoii  + 
0.66;  R^=0.71],  which  may  help  to  steer  the  way  towards  developing  a  temperature  and 
land  use  type  model  (similar  in  design  to  the  EPA’s  Biogenic  Emissions  Inventory 
System  (BEIS)  model  for  estimating  biogenic  NO  emissions)  to  estimate  biogenic  NH3 
emissions.  However,  when  the  major  rain  events  were  eliminated,  the  role  of  soil 
nitrogen  (both  NH3-N  and  TKN)  in  explaining  the  variability  in  NH3  flux  improved 
significantly. 

The  average  NH3  flux  values  from  this  study  corresponded  well  with  other 
reported  values  and  confirmed  that  soils  have  smaller  fluxes  of  NH3  than  lagoons  (Table 
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5.2).  However,  as  in  most  commercial  animal  operations  the  land  area  used  for  crops  is 
significantly  larger  than  the  surface  area  for  lagoons.  A  preliminary  analysis  revealed 
that  given  the  relative  sizes  of  the  agricultural  soils  in  comparison  to  the  lagoons,  the  soils 
(soils  represent  ~28%  of  the  lagoon  NH3  emissions)  cannot  be  neglected  when 
developing  and  apportioning  NH3  budgets.  Further,  the  soil  algorithm  developed  in  this 
study  was  based  on  a  crop  which  had  not  been  fertilized  for  several  months  prior  to  the 
experimental  period.  Given  that  land  spreading  of  slurry  is  estimated  to  release 
approximately  15-76%  of  the  applied  nitrogen  may  indicate  that  estimates  presented  in 
this  study  should  be  considered  as  a  lower  limit.  It  should  be  noted  that  the  data  used  to 
produce  this  inventory  was  only  from  this  one  field  site  during  the  time  periods 
discussed.  The  exponential  temperature  relationship  reported  here,  however  has  been 
consistently  reported  throughout  the  literature  and  therefore  provides  some  basis  for  this 
study  to  be  extended  throughout  the  year  and  at  different  soil  and  crop  types.  It  should 
not  be  assumed,  however  that  this  temperature  dependence  can  be  extended  to  all 
temperatures,  as  temperatures  outside  the  range  of  15-35  °C  are  often  found  to  alter  the 
often-cited  exponential  relationship. 

While  temperature  has  often  been  shown  to  be  a  controlling  variable  in  nitrogen 
trace  gas  emissions,  other  variables  (pH,  N-content,  %  WFPS)  have  displayed  more 
mixed  results  (Williams  et  al.,  1992;  Matson  et  al,  1997).  Therefore,  future  work  should 
consist  of  field  and  laboratory  studies  to  further  investigate  these  relationships  and  data 
following  slurry  application  should  also  be  obtained  to  help  refine  the  NH3  budget  for 
intensively  managed  agricultural  soils.  While  this  research  does  help  to  quantify  the  soil 
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source  strength,  it  does  not  shed  much  light  on  transport  through  canopies  or  on  the  effect 
ambient  NH3  concentrations  have  on  deposition  versus  emissions.  Even  though  the 
values  reported  in  this  study  are  within  the  range  of  other  reported  values,  a  side-by-side 
comparison  of  the  differing  NH3  flux  methodologies  would  be  extremely  beneficial  in 
furthering  our  knowledge  of  this  species.  A  simple  model  based  on  first  principles  will 
be  proposed  and  analyzed  in  the  following  chapter. 
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Researcher 

Site  Description 

Reported  NH3  Flux 
Values 

Harrison  et  al.,  1989 

Grass  and  Crop  Surfaces 

-20  to  1 00  ng  N  m'^  s'' 

Meixner  et  al.,  1991 

Wheat 

-12  to  25  ngN  m'^  s'* 

Weber  et  al.,  2001 

Post-Fertilization  (80  kg  N  ha"') 
Winter  Wheat 

1999 

2000 

64  ng  N  m'^  s'' 
178  ngN  m'^  s'' 

This  Study 

■  Spring 

Corn  Crop,  pre-fertilization 

38  to  271  ngN  m'^  s'' 

■  Winter 

No  Crop  Planted 

3  to  26  ng  N  m'^  s'' 

Aneja  et  al.,  2000 

Anaerobic  Lagoon  Surface 

■  Spring 

1706  +  552  pgNm'^nin'' 

■  Summer 

4017  +  987  pgNm'^min'' 

■  Fall 

844  ±401  pg  N  m'^  min'' 

■  Winter 

305  ±  154  pgN  m'^  min'' 

Table  5.2.  List  of  researchers  and  reported  NH3  flux  values  measured  under  various 
crop  and  fertilization  scenarios. 
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Forage/Feed  N 
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Figure  5.1.  Major  routes  for  NH3  emissions  from  intensively  managed  animal  operations  in  North  Carolina,  USA.  Source:  Aneja 
et  al.,  2001b. 


upper  Coastal  Plain  Research  Station 
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Figure  5.2.  Site  of  Upper  Coastal  Plain  Research  Station. 


•  NH3  Flux 


one  standard  deviation  of  the  ammonia  flux. 
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Soil  Temperature  (°C) 

Figure  5.4.  Plot  of  NH3-N  Flux  versus  soil  temperature.  Vertical  lines  represent  one  standard  deviation  of  the  average  NH3 
Flux,  and  horizontal  bars  represent  one  standard  deviation  of  the  temperature  during  that  flux  measurement  period. 


NHs-N  (fig/g)  in  the  soil 

figure  5.5a.  Plot  of  NH3-N  Flux  versus  NH3-N  content  of  the  soil.  Note  that  (■)  represents  the  di 
events  >0.3  cm  and  (•)  represents  those  data  points  that  were  measured  on  days  having  rain  >0.3  ci 


■  Flux  for  all  data  points  (excluding  rain 
events>0.3  cm) 


600  650  700  750  800  850  900  950  1000  1050 

TKN  (pg/g)  in  the  soil 

Figure  5.5b.  Plot  of  NH3-N  Flux  versus  Total  Kjeldahl  Nitrogen  content  of  the  soil.  Note  that  (■)  represents  the  data 
points  without  rain  events  >0.3  cm  and  (•)  represents  those  data  points  that  were  measured  on  days  having  rain  >0.3  cm 


■  Emissions  from  Soil 
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Figure  5.6.  Budget  of  NH3-N  estimates  from  soil  and  lagoon  surfaces  at  the  Upper  Coastal  Plain  Research  Site.  Soil 
estimates  were  determined  by  applying  the  flux  algorithm  developed  in  this  study  to  the  total  land  area  (-101  ha).  Lagoon 
estimates  were  determined  by  applying  the  flux  algorithm  developed  by  Aneja  et  al.  (2000)  to  the  total  lagoon  surface  area 
(~  1  ha). 


Chapter  VI.  Modeling  of  Ammonia  Emissions  from  Soils 


Abstract 

Using  a  dynamic  flow-through  chamber  system  in  conjunction  with  a  Thermo 
Environmental  17C  Chemiluminescence  ammonia  (NH3)  analyzer,  emissions  from 
slurry-amended  (~  33  kg  N  ha‘‘)  and  non-amended  soils  were  calculated  at  a  swine  farm 
in  eastern  North  Carolina.  The  average  NH3-N  flux  values  during  the  period  when  the 
soils  were  not  amended  with  any  slurry  were  ~  54  ng  N  m'^  s'*  while  the  average  NH3-N 
flux  values  measured  immediately  following  the  application  of  slurry  to  the  soil  were 
1723.9  ng  N  m'^  s'*.  An  empirical  model  relating  soil  temperature  to  NH3  flux  for  non- 
amended  soils  explained  over  70%  of  the  variability  in  NH3  emissions,  however  a  similar 
empirical  model  relating  soil  temperature  to  NH3  flux  for  slurry  amended  soils  was  able 
to  explain  only  39%  of  the  variability  in  NH3  emissions.  A  mass  transport  model,  based 
on  physical  and  chemical  processes  to  estimate  NH3  emissions  from  recently  amended 
soils  is  also  presented  and  compared  and  contrasted  to  the  empirical  model.  The 
variables  used  in  the  mechanistic  model  are  pH,  soil  temperature  and  total  ammoniacal 
nitrogen  content.  When  using  the  mass  transport  model,  the  percentage  difference 
between  predicted  and  measured  values  for  the  non-amended  and  slurry-amended  soils 
were  164%  and  16  %  respectively,  indicating  that  the  mechanistic  model  is  only 
applicable  for  periods  when  nitrification/denitrification,  plant  uptake  and  immobilization 
are  small  enough  in  comparison  to  the  chemical  and  physical  processes  following  slurry 
application  that  they  can  be  ignored.  The  percentage  of  the  nitrogen  (N)  applied  which 
was  emitted  as  NH3  increased  at  its  greatest  rate  immediately  following  slurry  application 
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(1-2  days)  and  then  began  to  level  out  at  a  value  of  approximately  20%  by  day  4. 

Previous  laboratory  studies  found  these  volatilization  events  to  be  short  lived  (few  days-2 
weeks)  and  this  study  corroborates  those  findings. 
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Introduction 


Emissions  of  ammonia  (NH3)  from  soils  has  been  found  to  be  a  significant  source 
of  NH3  into  the  atmosphere  (Table  1 .2).  The  vast  majority  (approx.  90%)  of  the  NH3  in 
the  atmosphere  is  converted  into  NH/  aerosols  by  the  irreversible  reactions  of  ammonia 
with  sulfuric  acid  (H2SO4),  nitric  acid  (HNO3),  hydrochloric  acid  (HCl),  water  and  to  a 
lesser  extent  (approximately  10%)  the  hydroxyl  radical  (OH)  in  the  atmosphere 
(Wameck,  2000).  These  NH/  aerosol  producing  reactions  can  be  summarized  as 
follows: 


NH3(g)  +  H2S04(1)  ^  NH4HS04(1) 

(1) 

NH3(g)  +  HN03(g)  0  NH4N03(s) 

(2) 

NH3(g)  +  HCl(g)  0  NH4C1(s) 

(3) 

NH3(g)  +  H20(1)  NH4^+  OH' 

(4) 

NH3  +  OH  ^  NH2  +  H20(1) 

(5) 

(Finlayson  Pitts  and  Pitts,  1986). 

The  conversion  of  NH3  to  NH4’^  dictates  the  spatial  scale  of  the  contribution  of 
NH3  to  the  total  atmospheric  nitrogen  input  (Aneja  et  al.,  2001a).  Due  to  ammonia’s 
relatively  shorter  lifetime  in  the  atmosphere  (less  than  5  days),  low  source  height  and 
high  deposition  velocities,  its  distribution  is  usually  limited  to  its  nearby  surroundings 
(Warneck,  2000;  Aneja  et  al.,  2001a).  The  NH4’^  aerosols,  however  have  lifetimes  on  the 
order  of  1 5  days  and  therefore  can  travel  and  deposit  at  larger  distances  from  the 
ammonia  source. 
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Control  strategies  to  minimize  NH3  emissions  require  a  more  thorough 
understanding  of  the  relative  source  strengths  of  the  various  emission  pathways. 
Furthermore,  computers  used  to  model  the  fate  of  this  increased  ammonia  can  only  be  as 
accurate  as  the  input  data,  and  failing  to  capture  any  of  the  spatial  or  temporal  variability 
in  the  NH3  emissions  will  result  in  poor  model  output.  This  study  describes  a  mass 
transport  model,  relating  NH3  release  from  the  soil  to  the  soil  nitrogen  content,  pH  and 
soil  temperature.  The  model  will  then  be  assessed  with  field  measurements  with  the  aim 
of  gaining  a  better  understanding  of  the  physical  and  chemical  processes  controlling  the 
NH3  emissions  from  soils  and  the  applicability  of  this  proposed  mass  transfer  model. 
Methods  and  Materials 
Sampling  Site  and  Sampling  Scheme 

The  NH3  flux  measurements  were  made  at  the  Upper  Coastal  Plain  Research 
Station,  located  in  Edgecombe  County,  North  Carolina  (See  Figure  5.2  for  measurement 
location;  See  Table  6.1  for  site/soil  characteristics).  This  facility  (contains  approximately 
178  hectares,  101  of  which  are  cropland  soils)  is  operated  with  typical  agronomic  and 
husbandry  practices  for  the  respective  crops  and  animals.  The  facility  also  maintains  a 
farrow-to-finish  swine  operation  with  approximately  1250  swine  on  site.  The  waste  from 
the  animals  (urine  and  feces)  is  flushed  from  the  swine  production  houses  into  two 
uncovered  anaerobic  waste  treatment  and  storage  lagoons  (a  primary  and  secondary 
lagoon,  total  acreage  ~  1  hectare).  The  effluent  from  these  lagoons  is  periodically 
sprayed  to  the  crops  as  a  nutrient  source.  A  corn  crop  was  harvested  on  August  2U*  and 
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the  stalks  were  shredded  and  left  on  the  soil  surface.  No  cover  crop  was  planted  nor  was 
the  field  fertilized  throughout  the  winter  and  Spring  of  2001. 

NH3  concentration  measurements  were  made  on  13  randomly  selected  plots 
located  within  a  10  m  radius  of  a  mobile  laboratory  (Modified  Ford  Aerostar  Van, 
temperature  controlled  to  within  the  operating  range  of  the  instruments).  Portions  of  the 
research  facility  were  being  sprayed  with  the  slurry  during  this  measurement  period.  For 
the  purposes  of  this  study  and  to  avoid  contamination  of  the  research  equipment,  the  plots 
used  in  this  study  were  amended  individually.  Twelve  of  the  plots  were  surface  applied 
with  750  ml  of  slurry  collected  from  the  hog  waste  lagoon  and  the  last  plot,  which  was 
used  as  the  control,  was  unamended.  Based  on  chemical  analysis  conducted  on  April 
25“^,  2001  by  the  North  Carolina  Agronomy  Division,  this  equated  to  approximately  33 
kg  N  ha’*,  an  amount  representative  of  typical  agronomic  practices  (Troeh  and 
Thompson,  1993).  The  daily  sampling  scheme  was  identical  to  the  procedures  described 
in  Chapter  V. 

Instrumentation  and  Flux  Calculation 

The  chamber  design,  associated  mass  balance  equation  and  calibration  procedures 
are  described  in  full  in  Chapter  I. 

Soil  Analysis 

A  soil  sample  (top  2  cm)  was  taken  from  the  center  of  the  chamber  footprint  at  the 
end  of  each  measurement  period  (approximately  1  sample  per  day),  and  analyzed  for  soil 
pH,  soil  moisture,  NHx  (NHx=NH/+NH3),  NO3-N  and  Total  Kjeldahl  Nitrogen 
(TKN=organic  N+NH3-N+NH4’*^-N)  by  the  North  Carolina  State  University  Department 
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of  Biological  and  Agricultural  Engineering.  Based  on  previous  studies,  the  unamended 
plots  are  found  to  remain  fairly  consistent  in  terms  of  pH  and  N-content  and  therefore  the 
unamended  plot  was  only  sampled  once.  Soil  temperature  was  measured  with  a 
Campbell  Scientific  temperature  probe  (accuracy  ±  3%)  inserted  into  the  soil  to  a  depth 
of  approximately  5  cm.  Air  temperatures  (Campbell  Scientific;  accuracy  ±  3%)  were 
measured  inside  of  a  radiation  shield  at  a  height  of  1 .5  meters.  Data  was  stored  in  15 
minute  binned  averages  utilizing  a  Campbell  Scientific  21X  Micrologger,  in  conjunction 
with  a  Toshiba  laptop  computer. 

Model 

The  exchange  of  NH3  gas  from  the  soil  into  the  atmosphere  is  assumed  to  be 
related  to  the  resistances  in  both  the  liquid  and  gas  phases  and  the  gas  concentration 
gradient  between  the  soil-air  interface.  The  mechanistic  model  describing  this  exchange 
can  be  written  as: 

F1ux(NH3),„  =  (6) 

where  K  is  the  transfer  coefficient  from  the  NH3  gas  in  the  soil  to  the  NH3  gas  in  the  air 
(Singh  and  Nye,  1986a).  Total  ammonia  (NHx=NH3  +  NH4^)  in  the  soil  is  dependent  on 
the  plant  uptake  rate,  immobilization  rate,  nitrification/denitrification  rates  and 
volatilization  of  the  gas  from  the  soil  (Sherlock  and  Goh,  1985).  Previous  researchers 
have  shown  that  during  the  time  period  immediately  following  fertilization  (4  days-2 
weeks),  the  volatilization  process  is  sufficiently  strong  so  that  other  processes  (plant 
uptake,  immobilization,  nitrification/denitrification)  can  be  neglected  (Sherlock  and  Goh, 
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1985).  Furthermore,  following  slurry  applieation  the  son)  will  generally  be 

significantly  larger  than  the  the  NH3  flux  equation  can  be  rewritten  as: 

Flux(NH3\„=K[NH3t^„„„  (7) 


Therefore  the  determination  of  the  NH3  flux  from  the  soil  is  dependent  on  knowing  the 
NH3  gas  concentration  in  the  soil/slurry  environment  and  the  exchange  coefficient  K. 

The  NH3  gas  concentration  in  the  soil  solution  can  be  calculated  by  examining  the 
chemical  equilibrium  between  NH3  and  NH4’^.  As  slurry  is  initially  applied  to  the  soil  it 
is  rapidly  hydrolyzed  to  produce  ammonium  ions  in  the  soil: 

(Urea) - >NH4^aq)  (8) 

The  subsequent  dissociation  of  in  the  solution  can  be  described  by: 

NH^^+H^O  S  NH3(3q)+H30^aq)  (9) 

where  K.,„  ,  is  the  equilibrium  constant  and  is  given  by: 

INii  t 


K 


NH3* 


[NH; 

5(aq)l 

H3O 

'^(aq) , 

NH, 

+ 

4  (aq)_ 

(10) 


The  ammonium  equilibrium  constant  is  found  to  be  dependent  on  temperature  through 
the  following  equations  described  by  Bates  and  Pinching,  1950;  Hales  and  Drewes,  1979; 
Aneja  et  ah,  2001b. 


=  5.67*10”'®  ♦expl 


-6286* 
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298.15 


or  similarly 
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I  \  2729.92 

Log(KHHr)= -0-09018 - — 


(12) 


where  T  is  the  soil  temperature  in  Kelvin. 

This  first  order  equilibrium  between  NH/(aq)  and  NH3(aq)  is  extremely  fast  and  in  terms  of 
the  volatilization  it  will  not  be  rate  limiting  (Sherlock  and  Goh,  1985). 

The  total  ammoniacal  nitrogen  content  of  the  soil  (NHx(aq))  can  be  rewritten  as: 


and  substituting  equation  (13)  into  equation  (10)  yields: 


(13) 


(14) 


Solving  for  NH3(aq),  equation  (14)  can  be  rewritten  as: 

[nh„„,] 


[NH3,„)]  = 


3(aq)J  f  \ 

i.iv: 

^  ^NH/  J 


(15) 


Recognizing  that  H30^  will  also  hydrolyze  to  produce  ions  and  that  H’^=10’'’”, 
equation  (15)  can  be  rewritten  in  the  form: 


[nH3,.„]  = 


1  + 


^x(aq)J 

10~P” 


(16) 


Substituting  for  the  equilibrium  constant  ^  from  equation  (12),  we  are  left  with: 


NH 


[nh.i.,,1 


3(aq) 


(0.09018+^™— -pH) 

1  +  10  T 


(17) 
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At  the  interface  of  the  air  and  liquid  film  in  the  soil  solution,  the  aqueous  NH3 
concentration  is  assumed  to  be  related  to  the  gaseous  NH3  concentration  through  Henry’s 
Law  equilibrium  where; 

NH3„=HNH,<,„  (18) 

and  the  Henry’s  Law  equilibrium  constant  H  can  be  expressed  as: 

1477  7 

Log(H”' )  =  -1 .69  +  — (Hales  and  Drewes,  1979). 

Substituting  the  Henry’s  Law  relationship  (equation  18)  into  the  flux  equation  (7)  yields 
the  following  expression  for  NH3  Flux: 

FIux(nH,„,)=k[nH„.,,]h  (19) 

Substituting  the  from  equation  (17)  into  equation  (18)  yields  the  following  NH3 

flux  equation: 

F1ux(NH )= - ^20) 

^  ''  (0.0901 8+=^=^-pH) 

(I  +  IO  T 

The  mass  transport  model  given  in  equation  (20)  is  the  same  model  used  by  Aneja 
et  al.  (2001b)  for  estimating  NH3  flux  from  lagoon  surfaces.  While  all  parameters  in 
equation  20  can  be  measured  or  calculated  in  both  the  lagoon  and  soil  environments,  the 
flux  is  critically  dependent  upon  accurate  representation  of  the  mass  transfer  coefficient 
(K),  and  the  value  of  the  Henry’s  law  coefficient  (H).  Aneja  et  al.  (2001b)  used  two-film 
theory  which  accounts  for  the  role  of  air  velocity  and  temperature  in  the  gas  (kg)  and 


178 


liquid  phase  (Icl)  resistance  to  NH3  transfer  across  an  interface  and  is  given  by  the 
following  equation: 


In  a  review  of  over  30  publications,  it  was  discovered  that  the  reported  values  of 
this  overall  mass  transfer  coefficient  ranges  from  1.3  x  10'  ms'  to  1.2  x  10'  ms'. 
Several  researchers  have  attempted  to  relate  this  overall  mass  transfer  coefficient  to 
various  properties  such  as  roughness  length,  friction  wind  velocity  and  height  of  the 
internal  boundary  layer  (van  der  Molen  et  al.,  1990;  Olesen  and  Sommer,  1993).  Using  a 
micrometeorological  technique,  Svensson  and  Perm  (1993)  calculated  mass  transfer 

•  <-3 

coefficients  from  soils  amended  with  manure  and  reported  values  ranging  from  4.3  x  10' 
m  s'’  to  1.2  X  10'^  m  s'’ .  Researchers  in  a  laboratory  using  swine  manure,  however 
reported  the  overall  mass  transfer  coefficient  to  range  between  1.3  x  10'®  m  s'’  and  5.2  x 
10'®  m  s'’  (Zhang,  1992).  Great  variability  does  exist  in  the  range  of  reported  values. 
However  conditions  during  which  these  measurements  were  made  do  differ.  In  general,  a 
review  of  these  published  reports  indicated  that  for  the  slurry  mixture  alone,  values  of  the 
overall  mass  transfer  coefficient  were  on  the  order  of  10'®  m  s'’  -  10'®  m  s’’  while  field 
and  laboratory  experiments  on  slurry/soil  mixtures  were  on  the  order  of  10'"’  m  s'’  -  10'^ 
m  s'’. 

Singh  and  Nye  (1986b)  described  a  controlled  laboratory  experiment  in  which  the 
exchange  coefficient  was  measured  under  a  variety  of  flow  and  pH  conditions.  These 
researchers  reported  that  the  value  of  K  increased  linearly  as  the  flow  rate  over  the  soil 
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surface  increased  from  0  to  1  1  min'^  However,  they  reported  that  as  flow  rates  increased 
beyond  1  1  min"',  the  rate  of  increase  becomes  smaller  and  begins  to  level  out  at  a  value 
of  3.7  X  10’^  m  s*'.  Sherlock  and  Goh  (1985)  theororized  that  in  soils,  water  movement 
and  the  diffusion  of  ions  is  probably  more  restricted  than  in  flooded  soils  or  water  bodies. 
Therefore  they  concluded  that  volatilization  of  NH3  from  non-flooded  soils  is  more  likely 
to  be  controlled  by  the  rate  of  diffusion  through  the  soil  than  on  wind  speed. 

Numerous  studies  conducted  by  Denmead  et  al.  (1974)  and  Beauchamp  et  al. 
(1978,  1982)  reported  no  positive  relationships  between  wind  speed  and  NH3  emissions 
from  non  flooded  soils.  Given  the  good  reproducibility  of  the  measurements  in  the 
experiments  conducted  by  Singh  and  Nye  (1986b)  (+6%),  their  experimentally  measured 
value  of  3.7  x  10'^  m  s'*  will  be  used  in  this  study.  For  comparison  purposes,  the  NH3 
flux  using  the  empirically  determined  K  value  based  on  temperature  and  wind  speed 
(equation  21)  and  as  described  by  Aneja  et  al.  (2001b)  will  also  be  presented. 

Results  and  Discussion 
Temperature 

Positive  relationships  between  the  release  of  nitrogen  trace  gases  from  soils  and 
soil  temperature  are  well  established  in  the  literature.  For  example,  the  solid  line  (plotted 
on  a  logarithmic  scale)  in  Figure  6.1  shows  this  temperature  dependence  in  an  empirical 
model  [Logio(NH3-N  Flux)  =  0.054*Tsoii  +  0.66;  R^=0.71]  developed  by  Roelle  and 
Aneja  (2001).  This  empirical  model  was  developed  from  measurements  made  at  the 
same  measurement  location  as  this  slurry-amended  study  and  was  conducted  during  2 
different  seasons  (winter  and  spring,  2000)  when  no  fertilizer  (slurry  or  chemically 
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derived)  was  applied  to  the  soil  (within  3  months  of  data  collection).  The  R  of  0.71  in 
this  study  is  consistent  with  many  nitrogen  trace  gas  experiments  from  various  soil  and 
crop  types  which  have  reported  values  ranging  between  0.42  and  0.87  (Sullivan  et  al., 
1996;  Thornton  et  al.,  1997;  Roelle  et  al.,  1999;  Roelle  and  Aneja,  2001). 

The  four  data  points  (solid  squares)  surrounding  the  solid  line  in  Figure  6.1 
represent  the  control  plots  (no  slurry  applied  to  the  plot)  from  this  (Spring  2001) 
measurement  campaign.  These  data  points  fit  the  general  form  of  the  model,  and  in  fact 
including  these  four  new  data  points  into  the  model  changes  the  R^  only  slightly  from 
0.71  to  0.70.  The  12  data  points  (solid  triangles)  in  this  graph  represent  the  data  from  the 
spring  2001  slurry  amended  plots.  Although  these  12  data  points  did  fall  within  the  range 
of  soil  temperatures  used  to  develop  the  empirical  model,  attempting  to  estimate  the 
fluxes  from  the  slurry  amended  plots  using  the  empirical  model  in  Figure  6.1  would 
result  in  an  underestimation  of  at  least  an  order  of  magnitude.  An  empirical  model 

'y 

developed  from  only  the  slurry  amended  soil  data  results  in  an  R  =  0.26  and  a  new 
empirical  relationship  taking  all  data  points  (amended  and  non-amended)  into  account 
results  in  a  decreased  R^  value  from  0.71  to  0.39. 

These  results  highlight  that  temperature  will  explain  over  70%  of  the  variability  in 
NH3  emissions  when  it  is  developed  from,  and  applied  to,  soils  that  have  not  been 
recently  amended  (fertilized).  However,  similar  empirical  models  relating  temperature  to 
NH3  flux  developed  from  recently  amended  soils  or  a  combination  of  both  amended  and 
non-amended  soils  reduces  the  explanatory  capability  of  the  model  to  26%  and  39% 
respectively.  This  suggests  that  when  nitrification/denitrification,  plant  uptake  and 


181 


immobilization  are  sufficiently  small  to  be  neglected  in  comparison  to  the  chemical  and 
physical  processes  during  time  periods  immediately  following  slurry  application,  than 
temperature  alone  can  not  adequately  estimate  the  NH3  flux.  During  non-fertilized 
episodes,  this  relationship  can  be  attributed  to  the  fact  that  the  biochemical  reaction  rates 
of  microorganisms  responsible  for  the  production/consumption  of  NH3  respond  to 
changes  in  temperature  (Sullivan  et  al.,  1996;  Roelle  et  al.,  1999;  Roelle  et  al.,  2000; 
Wameck,  2000).  However,  during  recently  fertilized  events,  as  shown  through  the 
development  of  the  mass  transfer  model,  parameters  other  than  temperature  such  as 
ammoniacal  nitrogen  content  and  pH  must  be  considered. 

Similarly  to  the  empirical  model,  the  mass  transport  model  relates  the  NH3  flux  to 
soil  temperature  (exponential  dependence).  Maintaining  all  other  parameters  constant 
(typically  observed  agronomic  values  before(i)  and  after(2)  slurry  application:  NHx(i)=6 
pg/g;  pH(i)=5.5;  K(i)=3. 69x10'^  m  s'’;  NHx(2)=90  pg/g;  pH(2)=6.5;  K(2)=3. 69x10'^  m  s'') 
and  varying  the  temperature  between  16.3  and  20.6  °C  (typical  range  during  a 
measurement  period)  in  both  the  mechanistic  and  empirical  model  produces  similar 
profiles,  yet  at  different  magnitudes  (Figure  6.2).  The  average  NH3  flux  measured  before 
and  after  slurry  application  was  54  and  1724  ng  N  m'^  s'',  respectively.  Although  soil 
temperatures  from  the  different  seasons  did  overlap,  it  is  apparent  that  a  temperature 
based  empirical  relationship  will  fail  to  capture  the  magnitude  of  the  emissions.  While 
Figure  6.2  reemphasizes  the  exponential  dependence  of  NH3  flux  on  soil  temperature,  it 
also  highlights  the  role  that  other  parameters,  namely  pH  and  ammoniacal  nitrogen  must 
play  in  the  mass  transfer  equation. 
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Applying  the  mechanistic  model  to  the  non-amended  field  site  data  collected 
during  the  spring  and  winter  2000  measurement  campaign  (Figure  6.3-first  10  data 
points)  results  in  much  weaker  model  performance  than  when  the  mechanistic  model  is 
used  for  time  periods  immediately  following  slurry  application  (Figure  3 -last  4  data 
points).  The  average  percent  difference  between  modeled  and  measured  values  for  the 
slurry  amended  plots  was  16%  while  the  average  percent  difference  for  non-amended 
plots  was  164%.  This  weaker  relationship  for  the  non-amended  plots  is  explained  by  the 
fact  that  the  assumptions  in  deriving  the  mechanistic  model  were  that 
nitrification/denitrification,  plant  uptake  and  immobilization  processes  were  much  slower 
than  the  volatilization  processes,  and  that  ambient  [NH3]  were  negligible,  all  of  which  are 
invalid  assumptions  during  the  unamended  measurement  period.  Therefore,  during  the 
slurry  amended  measurement  periods,  it  is  important  to  consider  the  other  controlling 
parameters  from  the  mass  transfer  equation. 

Based  on  measured  wind  speeds  and  temperature  during  this  measurement  period, 
the  mass  transfer  coefficient  ranged  from  1 .5  x  lO’^m  s’*  to  2.4  x  10'^  m  s’*  as  opposed  to 
the  experimentally  measured  value  of  3.69  x  10’^  m  s’*.  Results  using  the  same  soil 
conditions  and  the  calculated  mass  transfer  coefficient  from  equation  21  are  also  plotted 
in  Figure  6.3.  From  the  results  of  this  study,  the  experimentally  measured  mass  transfer 
coefficient  is  more  accurate  in  estimates  of  the  NH3  flux.  This  may  confirm  that  in  soil 
environments,  NH3  emissions  are  limited  more  by  rate  of  diffusion  of  the  NH3  to  the 
soil/air  interface  than  by  windspeed  or  temperature,  as  suggested  by  Sherlock  and  Goh 
(1985). 
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Ammoniacal  Nitrogen  Content  and  pH 

Figure  6.4  shows  both  the  daily  trend  of  the  NH3-N  content  of  the  soil  (primary 
axis)  and  the  averaged  NH3  flux  from  the  3  sample  plots  on  each  of  the  successive 
measurement  days  (secondary  axis).  The  solid  squares  in  Figure  6.4  represent  the 
average  flux  from  the  control  (unamended)  plots.  The  decreasing  trends  in  both  NH3-N 
flux  and  NHx  content  of  the  soil  are  expected  based  on  the  fact  that  as  volatilization  of 
NH3  from  the  soil  continues,  in  the  absence  of  other  NH3  production  mechanisms,  there 
will  be  progressively  smaller  concentrations  of  NH3  in  the  soil.  To  examine  how  changes 
in  [NHx]  effect  the  NH3  flux  predicted  by  the  mechanistic  model,  all  parameters  in  the 
mechanistic  model  are  kept  constant  (T=20  °C;  pH=6.5;  K=3.69xl0‘^  m  s’*)  and  only  the 
[NHx]  is  varied  (Figure  6.5).  This  modeled  linear  relationship  is  supported  by  an  earlier 
study  conducted  by  Roelle  and  Aneja,  (2001)  during  a  non-fertilized  period  which 
showed  the  NH3  flux  from  the  soil  being  highly  correlated  to  the  NH3-N  content  of  the 
soil  [NH3  Flux=55.5*(NHx)-160;  R^=0.86].  Using  this  empirical  relationship,  which  is 
based  on  soil  conditions  with  much  smaller  [NHx-N]  (<9  pg/g  compared  to  ~  105  pg/g) 
results  in  significantly  higher  NH3  flux  estimates  (solid  line)  than  the  mechanistic  model 
for  values  of  [NHx]  greater  than  approximately  5  pg/g. 

Producing  a  new  empirical  model  [NH3  Flux=22.7*(NHx)-  91.4;  R^  =  0.75]  taking 
both  the  slurry-amended  and  non-slurry  amended  soils  into  account  shows  that  the  soils 
ammoniacal  nitrogen  content  is  significant  (p<0.01)  in  estimating  NH3  emissions  both 
prior  to  and  after  slurry  application.  The  shaded  and  dashed  lines  in  Figure  6.5  differ  by 
approximately  40%  at  lower  values  of  NHx  (~10  pg/g)  and  less  than  7%  at  the  higher 
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values  of  NHx  (-120  |ag/g).  Although  it  appears  that  the  mechanistic  and  empirical 
models  yield  approximately  equivalent  flux  estimates,  it  should  be  noted  that  in  this 
example  the  pH  was  kept  constant  in  the  mechanistic  model.  While  many  field  studies 
have  found  pH  to  remain  fairly  uniform  throughout  non-fertilized  periods,  this 
assumption  is  invalid  during  periods  immediately  following  fertilization  (Sullivan  et  al., 
1996;  Roelle  et  al,  1999;  Li  et  al,  1999). 

As  slurry  is  initially  applied  to  the  soil  it  is  rapidly  hydrolyzed  to  produce 
ammonium  and  carbonate  ions  which  causes  an  increase  in  the  soil  pH.  As  the  ammonia 
volatilizes,  H"^  ions  accumulate  in  the  soil  and  throughout  the  volatilization  event  there  is 
a  gradual  decrease  of  the  soil  pH.  Therefore,  during  recently  fertilized  episodes,  pH  and 
[NHx]  can  not  be  considered  in  isolation  because  a  change  in  [NHx]  causes  a  change  in 
the  pH.  This  effect  can  be  seen  clearly  in  Table  6. 1 ,  which  shows  the  pH  begirming  at 
approximately  6.6  for  the  amended  plots  and  gradually  decreasing  to  6.2,  while  the 
control  plot  is  assumed,  based  on  previous  studies,  to  remain  fairly  constant  at  5.4. 

While  the  overall  change  during  these  4  days  is  less  than  ‘A  of  a  pH  unit,  the 
effects  on  the  estimated  ammonia  flux  can  be  large.  Based  on  the  mass  transport  model, 
changing  only  the  soil  pH  by  1  unit  and  maintaining  all  other  variables  constant  at 
(T=20°C,  NHx=90  |ag/g,  K=3.69  x  10'^  m  s'*)  results  in  changes  of  NH3  flux  of 
approximately  an  order  of  magnitude  (Figure  6.6-solid  line).  However,  in  a  more 
realistic  situation,  varying  both  the  pH  and  [NHx]  (as  shovm  on  the  secondary  axis) 
results  in  the  NH3  flux  as  seen  in  Figure  6.6  (solid  squares). 
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Figure  6.1  shows  the  measured  NH3-N  volatilized  as  a  %  of  the  N  applied.  At  a 
rate  of  N-application  of  33  Kg  N  ha'*,  nearly  20%  of  the  applied  N  is  lost  as  NH3  within 
the  first  4  days  after  application.  The  predicted  NH3-N  volatilized  as  a  %  of  the  N 
applied  is  shown  by  the  solid  line  in  Figure  6.7,  which  is  based  on  the  daily  averaged  soil 
parameters  measured  and  the  mass  transfer  equation  with  the  constant  mass  transfer 
coefficient  (3.69  x  10'^  m  s'*).  The  solid  line  appears  to  accurately  predict  the  percent  of 
N  lost  as  NH3  during  the  first  three  days  following  slurry  application.  However,  while 
the  measured  rate  of  loss  begins  to  level  out  at  approximately  20%,  the  solid  line 
(predicted  value)  tends  to  be  increasing  although  at  a  slower  rate.  Other  modeling  and 
experimental  studies  have  found  that  the  percent  of  N  applied,  which  is  lost  as  NH3, 
typically  levels  out  at  approximately  30%  and  that  this  usually  occurs  within  the  first  1-2 
weeks  after  fertilizer  application  (Singh  and  Nye,  1986b  and  1988).  While  the  NH3 
volatilized  (as  a  percent  of  the  N  applied)  in  this  study  appears  to  level  out  slightly  faster 
than  other  reported  studies,  this  may  be  a  factor  of  the  amount  of  nitrogen  initially 
applied  to  the  soil.  Whereas,  the  Singh  and  Nye  (1986a)  study  applied  nearly  210  kg  N 
ha'*,  this  study  applied  33  kg  N  ha'*.  Furthermore,  the  [NHx]  was  still  elevated  (in 
comparison  to  background  levels)  at  the  conclusion  of  the  measurement  period  indicating 
that  this  volatilization  event  may  have  persisted  for  a  few  more  days.  Regardless,  the 
consistent  trend  among  the  model,  measurements  and  previous  lab  and  field  studies 
indicates  that  these  volatilization  events  are  short-lived  (few  days-2  weeks),  and  generally 
result  in  emissions  of  20-30%  of  the  applied  nitrogen. 
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Conclusions  and  Recommendations 


Using  a  mechanistic/mass  transport  model  and  measured  soil  parameters, 
ammonia  emissions  were  estimated  and  compared  to  calculated  flux  values  at  an  eastern 
North  Carolina  swine  farm.  A  mechanistic  model  developed  for  volatilization  events 
performed  well  (average  of  16%  difference  between  modeled  and  measured  flux  values) 
immediately  following  a  slurry  application,  however  performed  poorly  on  the  non- 
amended  soils  (average  of  1 64%  difference  between  modeled  and  measured  flux  values). 
This  relationship  is  not  surprising  given  that  the  assumptions  in  the  mass  transport  model, 
namely  negligible  plant  uptake  and  ammonification  processes  relative  to  the  ammonia  in 
the  soil/slurry  mixture,  both  of  which  (plant  uptake  and  ammonification  processes)  are 
not  valid  assumptions  during  periods  prior  to  or  long  after  slurry/fertilizer  application. 
Further,  during  time  periods  prior  to  slurry  application  the  assumption  of  ammonia 
concentration  in  the  air  being  negligible  in  comparison  to  the  ammonia  concentration  in 
the  soil  cannot  be  assumed  to  be  valid  and  therefore  ammonia  concentrations  in  the  air 
must  be  included  in  the  mass  transfer  model. 

In  this  study,  the  value  of  the  mass  transfer  coefficient  was  kept  constant  at  3.7  x 
10'^  m  s'^  based  on  laboratory  results  reported  by  Singh  and  Nye  (1986b)  and  research 
conducted  by  other  investigators.  Using  an  equation  to  describe  the  exchange  coefficient 
in  terms  of  wind  speed  and  temperature  resulted  in  calculated  exchange  coefficients 
ranging  from  1.5  x  lO’^m  s"'  to  2.4  x  10'^  m  s'\  The  literature  currently  describes  mass 
transfer  coefficients  ranging  from  1.3  x  10'^  m  s"*  to  1.2  x  10'^  m  s"'  and  therefore  a  more 
thorough  understanding  of  this  parameter  must  be  obtained.  The  fact  that  the 
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experimentally  measured  exchange  coefficient  performed  better  than  the  calculated 
exchange  coefficient  may  indicate  that  for  soil  systems  the  NH3  flux  is  rate  limited  by 
diffusion  to  the  soil/air  interface  as  opposed  to  wind  speed  or  temperature.  Soil 
ammoniacal  nitrogen  content  was  found  to  be  linearly  related  to  NH3  while  pH  and 
temperature  were  both  found  to  have  an  exponential  dependence.  The  measurements 
revealed  that  the  applied  N  is  lost  at  the  greatest  rate  in  the  first  2  days  following 
application  and  begin  to  level  out  by  day  4,  although  additional  field  data  is  required  to 
confirm  this  relationship. 

While  the  currently  used  emission  factor  approach  may  adequately  capture  the 
total  NH3  emitted  to  the  atmosphere  on  a  yearly  basis,  these  results  show  that  they  would 
perform  poorly  in  resolving  any  temporal  or  spatial  trends.  In  an  effort  to  further  refine 
global  Nitric  Oxide  (NO)  emission  estimates,  Yienger  and  Levy  (1995)  proposed  a 
“pulsing”  mechanism  to  account  for  the  large  bursts  of  NO  following  the  wetting  of  dry 
soils.  Similarly,  it  would  appear  that  a  mechanism  should  be  incorporated  into  the 
emission  estimate  process  to  further  refine  the  budget  of  ammonia  emissions  from  soils, 
especially  intensively  managed  soils  which  are  consistently  amended  with  both 
commercially  derived  fertilizers  and  animal  waste.  As  a  first  approach  during  non- 
amended  periods,  the  temperature-based  model  appears  to  capture  the  majority  of  the 
variation  in  NH3  emissions.  A  possible  approach  for  the  amended  periods  would  be  to 
adopt  the  same  procedures  commonly  used  in  the  estimates  of  biogenic  NO  emissions, 
which  is  to  apply  a  temperature  algorithm  and  a  factor  to  adjust  for  the  amount  of 
fertilizer  the  crop  receives.  Unfortunately,  this  approach  requires  a  large  set  of  data 
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conducted  over  many  different  soil  and  crop  types  to  empirically  determine  these  factors, 
which  to  date  is  still  unavailable. 
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Site:  Upper  Coastal  Plain  Research  Station 
35,9°  N  Latitude 
11, T  W  Longitude 
Elevation:  35.5  Meters 
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Model  for  Unfertilized  Soils 
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Figure  6.1.  Measured  NH3-N  flux  versus  soil  temperature.  The  solid  line  represents  a  temperature  based  aigonthm 
developed  as  part  of  an  earlier  study  (Roelle  and  Aneja,  2001).  The  solid  squares  represent  the  control  plots  from  spring 
2001  measurements  while  the  solid  triangles  represent  the  NH3-N  flux  from  the  slurry  applied  soils.  Vertical  bars  represent 
one  standard  deviation  of  the  NH3  flux. 
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Mechanistic  Model  on  slurry  amended  soil 
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Figure  6.2.  Modeled  NH3-N  flux  versus  soil  temperature.  The  solid  line  represents  the  empirical  model  and  the  shaded(i) 
and  dashed(2)  lines  represent  a  mechanistic  model  with  all  parameters  constant  at:  NHx(i)=6  pg/g;  pH(i)=5.5;  K^i)=3.69  x  10 
m  s‘*;  NHx(2)=90  pg/g;  pH(2)=6.5;  K(2)“3.69  x  10’^  m  s'*.  Soil  temperature  was  varied  between  16.3  and  20.6  °C. 


♦  Measured  (non-amended) 
■  Predicted(l) 
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Figure  6.3.  Plots  of  measured  and  predicted  (mechanistic  model)  NH3-N  flux  on  different  measurement  days.  The  first  10 
data  points  represent  the  non-amended  plots,  while  the  last  4  data  points  represent  the  slurry  amended  plots.  Predicted(l) 
are  estimated  using  K=3.69  x  10'^  m  s’'  while  Predicted(2)  calculates  K  from  equation  (21). 
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averaged  NH3  flux  from  the  control  (unfertilized)  plots. 


Mechanistic  Model 


s  J.IU  N  Bu)  xnid  ^hN 


198 


Figure  6.5.  NH3  flux  versus  ammoniacal  nitrogen  content  of  the  soil.  The  solid  line  represents  the  mechanistic  model  and 
the  dashed  and  shaded  lines  represents  empirical  models  from  the  non-amended  and  combined  data  sets  repectively.  All 
values  kept  constant  at  T=20°C;  pH=6.5;  K=3.69  x  10'^  m  s‘\ 
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Figure  6.6.  NH3  flux  versus  soil  pH.  Solid  line  represents  the  mechanistic  model,  and  only  varying  the  pH.  All  values  in 

the  model  are  kept  constant  at  T=20°C;  NHx=90  pg/g;  K=3. 69x10’^  m  s’^  The  solid  squares  (■)  represent  the  mechanistic 
model,  varing  both  the  pH  and  the  [NHx]  in  the  soil  (•),  as  indicated  on  the  secondary  axis. 
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Figure  6.7.  NH3-N  volatilized  as  a  percent  of  the  N  applied.  Solid  squares  represent  the  measured  values  from  each  plot  on 
the  4  measurement  days.  The  solid  line  represents  the  modeled  NH3-N  volatilized  using  the  soil  parameters  measured. 


Chapter  VII.  Summary  and  Suggestions  for  Future  Research 


A  dynamic  flow-through  chamber  technique  was  used  to  measure  oxidized  and 
reduced  biogenic  nitrogen  compound  emissions  into  the  rural  troposphere.  This  research 
initially  focused  on  soils  amended  with  chemically  derived  fertilizers  and  investigated  the 
chemical  and  physical  parameters  which  are  often  cited  as  being  significant  in  controlling 
the  emissions  from  soils.  NO  emissions  and  soil  properties  were  studied  at  several 
croplands  throughout  North  Carolina  and  were  generally  found  to  follow  a  diurnal  profile 
with  maximum  emissions  coinciding  with  maximum  soil  temperatures.  The  exponential 
dependence  of  NO  flux  on  soil  temperature  existed  at  all  sites,  but  to  different  levels  of 
significance.  It  was  also  observed  that  NO  flux  did  respond  to  varying  amounts  of  both 
total  extractable  nitrogen  (TEN)  in  the  soil  and  soil  moisture  content. 

Throughout  the  1990’s  two  processes  were  identified  as  potentially  contributing 
significant  nitrogen  to  the  ecosystem  throughout  North  Carolina.  Considering  the 
experience  that  our  research  group  had  in  measuring  nitrogen  trace  gas  emissions  it  was 
decided  that  these  two  relatively  new  nitrogen  sources  should  be  characterized  to 
determine  their  source  strength  and  modeled  to  determine  possible  impacts.  The  first  of 
these  processes  was  the  land  application  of  biosolids,  which  are  the  byproduct  of  the 
wastewater  treatment  process.  These  biosolids  are  land  applied  for  their  nutrient  content 
and  as  a  cost  effective  way  for  the  wastewater  treatment  plant  to  dispose  of  this 
byproduct.  Measurements  were  made  at  a  field  site  amended  with  biosolids  throughout 
1999-2000  and  NO  emissions  were  found  to  be  dependent  primarily  on  soil  temperature. 
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The  exponential  dependence  of  NO  emissions  on  soil  temperature  is  the 
relationship  most  often  cited  in  the  literature  and  has  been  found  to  be  consistent  across 
soil  and  crop  types.  Using  this  temperature  relationship  and  also  taking  into 
consideration  the  nitrogen  fertilizer  being  applied  to  the  crop  results  in  an  NO  emissions 
model  which  is  described  by  Williams  et  al.  (1992)  (Chapter  III).  The  Williams  model 
combines  soil  temperatures  and  a  land  use  database  to  produce  NO  emission  inventories 
which  can  then  be  used  in  air  quality  models.  Currently,  land  use  databases  have 
multiple  categories  for  different  crop  and  vegetation  types  but  there  is  no  category  for 
soils  which  are  amended  with  biosolids.  Considering  that  this  study  found  biosolid 
amended  soils  to  have  high  emissions  of  NO,  relative  to  soils  amended  with  chemically 
derived  fertilizers,  indicates  that  biogenic  emission  inventory’s  can  be  underestimated 
and  can  therefore  cause  biases  in  air  quality  models. 

In  this  study  we  have  determined  how  many  acres  of  biosolid  amended  soils  exist 
in  the  various  counties  of  North  Carolina  and  have  applied  our  observationally  based 
temperature  algorithm  (Chapter  III)  to  determine  a  modified  NO  emissions  inventory.  It 
was  not  possible  to  determine  where  within  the  counties  these  biosolids  were  being 
applied,  so  therefore  we  had  to  use  our  discretion  in  apportioning  this  land-use 
throughout  the  grid  cells  of  the  various  counties.  We  used  two  approaches,  the  first  of 
which  was  to  evenly  distribute  the  bio  solids  throughout  all  the  grid  cells  of  the  county 
and  the  other  was  to  concentrate  them  within  the  county.  The  Multiscale  Air  Quality 
Simulation  Platform  (MAQSIP)  was  then  run  in  two  scenarios,  the  first  of  which  we  call 
the  base  case  and  the  second  being  the  modified  case.  In  the  base  case,  the  current 
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biogenic  NO  emissions  inventory  is  used  and  in  the  modified  case  MAQSIP  is  run  with  a 
modified  NO  emissions  inventory  which  takes  the  biosolid  amended  acreages  into 
consideration.  It  was  found  that,  in  general,  the  greatest  change  in  the  model  output 
occurred  in  the  evening  and  consisted  of  ozone  being  depleted.  Concentrating  the 
biosolids  within  a  few  grid  cells  of  the  respective  counties  resulted  in  changes  where 
ozone  was  depleted  by  as  much  as  1 1%  in  the  evening  and  increased  by  approximately 
2%  in  the  afternoon.  While  neither  scenario  is  likely  to  be  representative  of  actual 
conditions,  this  approach  does  serve  as  boundaries  with  which  the  significance  of  this 
land  use  class  can  be  assessed.  These  results  indicate  that  the  biosolids  can  have 
significant  impacts  on  model  results  and  therefore  data  needs  to  be  made  available 
describing  where  these  biosolids  are  being  applied  and  then  they  need  to  be  included  in 
future  land  use  databases. 

Future  research  concerning  the  emissions  of  nitric  oxide  from  biosolid  amended 
soils  should  focus  on  emissions  from  various  crop  and  soil  types  amended  with  bioslids 
to  verify  the  relationships  reported  in  this  study.  A  comparison  of  these  field  results  and 
a  similarly  designed  laboratory  study  found  the  differences  between  the  measured 
emissions  to  be  statistically  significant.  While  this  initial  comparison  looked  at  soil 
temperature  and  soil  moisture,  it  did  not  consider  nitrogen  content  and  microbial 
population.  Therefore  a  more  detailed  comparison  should  be  conducted  to  assess  how 
these  factors  (nitrogen  content  and  microbial  population)  in  the  field  and  laboratory 
environments  may  act  to  control  the  emissions  of  NO. 
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The  second  process  which  is  contributing  increased  nitrogen  to  the  environment 
of  North  Carolina  is  the  land  application  of  slurry  from  swine  facilities.  During  the 
1990’s  the  swine  population  in  North  Carolina  increased  from  approximately  3  million 
animals  to  approximately  1 0  million  today,  making  North  Carolina  the  second  largest 
swine  producing  state  in  the  U.S.  Therefore,  in  an  attempt  to  quantify  the  strength  of  the 
ammonia  being  emitted  via  the  swine  operations,  ammonia  emissions  were  calculated 
during  various  seasons  and  for  time  periods  prior  to  and  after  the  slurry  was  applied  to  the 
soil.  An  observationally  based  model  to  estimate  ammonia  emissions  prior  to  slurry 
application  was  developed  and  then  compared  to  a  similar  model  for  estimating  emissions 
from  lagoon  surfaces.  A  comparison  of  these  two  models  revealed  that  on  a  yearly  basis, 
ammonia  emissions  from  the  soil  contribute  approximately  28%  of  the  lagoon  emissions. 
Before  the  results  of  this  one  study  are  extrapolated  to  farms  throughout  North  Carolina, 
similar  studies  should  be  conducted  at  other  swine  facilities.  However,  given  that  the 
temperature  dependence  of  nitrogen  trace  gases  has  consistently  been  shown  to  follow 
this  exponential  relationship  across  soil  and  crop  types  and  during  different  times  of  the 
year,  it  is  expected  that  these  results  will  be  similar  at  other  locations. 

Ammonia  emissions  were  also  measured  for  time  periods  immediately  following 
slurry  application  and,  as  expected,  the  emissions  were  larger  by  approximately  an  order 
of  magnitude.  For  time  periods  prior  to  slurry  application,  soil  temperature  alone  was 
found  to  be  the  best  predictor  of  ammonia  flux,  explaining  over  70%  of  the  variability  in 
the  ammonia  emissions.  However  for  time  periods  after  slurry  application,  soil 
temperature  explained  approximately  40%  of  the  ammonia  emissions.  These  results 
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indicate  that  microbial  processes  are  responsible  for  the  majority  of  the  ammonia  in  the 
soil  prior  to  slurry  application.  In  an  attempt  to  estimate  ammonia  emissions  for  time 
periods  after  slurry  application,  a  fundamental  mechanistic  mass  transfer  model  was 
investigated  and  assessed  with  the  field  data  from  this  site.  This  mass  transfer  model  was 
found  to  predict  the  ammonia  emissions,  within  approximately  15%,  but  only  for  time 
periods  immediately  following  slurry  application. 

Future  research  concerning  the  emissions  of  ammonia  from  soils  amended  with 
lagoon  effluent  should  consist  of  additional  measurements  at  various  swine  farms  to 
verify  the  relationships  reported  here.  Further,  an  extended  study  for  time  periods  after 
slurry  application  should  be  conducted  to  determine  the  length  of  time  that  emissions 
remain  elevated  following  application  of  the  slurry.  Once  this  relationship  is  determined, 
a  better  ammonia  inventory  can  be  developed  for  this  region  which  takes  into 
consideration  not  only  microbial  production  of  ammonia  in  the  soil  but  also  the 
volatilization  events  associated  with  lagoon  effluent  application. 
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Appendix  1.  Variations  of  Nitric  Oxide  Fluxes  from  Diverse 
Physiographic  Agricultural  Soils  in  North  Carolina 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Emissions  of  nitric  oxide  (NO)  were  determined  during  late  spring  and  summer  1995  and 
the  spring  of  1996  from  four  crop  types,  located  at  four  different  physiographic  regions  in 
North  Carolina.  Emission  rates  were  calculated  using  a  dynamic  flow-through  chamber 
system  coupled  to  a  state-of-the-art  mobile  laboratory  for  in-situ  analysis.  Average  NO 
fluxes  during  late  spring  1995  were;  50.9  ±  47.7  ng  N  m'^  s'*  for  com  in  the  lower  coastal 
plain.  Average  NO  fluxes  during  summer  1995  were:  6.4  ±  4.6  and  20.2  ±  19.0  ng  N 
s'*  respectively  for  com  and  soybean  in  the  coastal  region;  4.2  ±  1.7  ng  N  m'^  s'*  for 
tobacco  in  the  piedmont  region;  and  8.5  +  4.9  ng  N  m'^  s'*  for  corn  in  the  upper  piedmont 
region.  Average  NO  fluxes  for  spring  1996  were:  66.7  ±  60.7  ng  N  m'  s'  for  wheat  in 
the  lower  coastal  plain;  9.5  ±  2.9  ng  N  m'^  s'*  for  wheat  in  the  coastal  plain;  2.7  ±  3.4  ng 
N  m'^  s'*  for  wheat  in  the  piedmont  region;  and  56.1  ±  53.7  ng  N  m's'*  for  corn  in  the 
upper  piedmont  region.  An  exponential  dependence  of  NO  flux  on  soil  temperature  was 
present  at  all  of  the  locations.  Further,  all  locations  displayed  a  diurnal  trend  of  NO 
emissions  which  revealed  a  peak  in  NO  emissions  that  coincided  with  the  maximum  soil 
temperature  for  the  day.  The  composite  data  of  all  the  research  sites  revealed  a  general 
positive  trend  of  increasing  NO  flux  with  soil  water  content  and  extractable  nitrogen. 

Presented  at: 

Workshop  on  Atmospheric  Nitrogen  Compounds  II:  Emissions,  Transport, 
Transformation,  Deposition  and  Assessment 
Chapel  Hill,  North  Carolina 
June  7 -9, 1999 
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Appendix  2.  Biogenic  Nitric  Oxide  Source  Strengths 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospherie  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Emissions  of  nitric  oxide  (NO)  were  determined  during  late  spring  and  summer  1995  and 
the  spring  of  1 996  from  four  crop  types,  located  at  four  different  physiographie  regions  in 
North  Carolina.  Emission  rates  were  calculated  using  a  dynamic  flow-through  chamber 
system  coupled  to  a  state-of-the-art  mobile  laboratory  for  in-situ  analysis.  Average  NO 
fluxes  during  late  spring  1995  were:  50.9  ±  47.7  ng  N  m'^  s''  for  com  in  the  lower  coastal 
plain.  Average  NO  fluxes  during  summer  1995  were:  6.4  ±  4.6  and  20.2  ±  19.0  ng  N  m’^ 
s  respectively  for  com  and  soybean  in  the  coastal  region;  4.2  ±  1.7  ng  N  m‘  s'  for 
tobacco  in  the  piedmont  region;  and  8.5  ±  4.9  ng  N  m'^  s''  for  com  in  the  upper  piedmont 
region.  Average  NO  fluxes  for  spring  1996  were:  66.7  ±  60.7  ng  N  m'  s'  for  wheat  in 
the  lower  coastal  plain;  9.5  ±  2.9  ng  N  m'^  s''  for  wheat  in  the  coastal  plain;  2.7  ±  3.4  ng 
N  m'^  s''  for  wheat  in  the  piedmont  region;  and  56.1  ±  53.7  ng  N  m's''  for  corn  in  the 
upper  piedmont  region.  An  exponential  dependence  of  NO  flux  on  soil  temperature  was 
present  at  all  of  the  loeations.  Further,  all  loeations  displayed  a  diurnal  trend  of  NO 
emissions,  which  revealed  a  peak  in  NO  emissions  that  coincided  with  the  maximum  soil 
temperature  for  the  day.  The  composite  data  of  all  the  research  sites  revealed  a  general 
positive  trend  of  increasing  NO  flux  with  soil  water  content  and  extractable  nitrogen 

Presented  at: 

Air  &  Waste  Management  Association’s  92nd  Annual  Meeting  and  Exhibition 
St.  Louis,  Missouri 
June  20 -24,  1999 
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Appendix  3.  Biogenic  Air  Pollution _ 

Viney  P.  Aneja  and  Paul  A.  Roelle 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

and  Jeffrey  Peirce 

Department  of  Civil  and  Environmental  Engineering 
Duke  University 
Durham,  NC  27708-0287 

Microorganisms  in  soil  and  in  aquatic  environments  play  an  extremely  important 
role  in  the  cycling  of  atmospheric  pollutants,  but  one  which  is  difficult  to  assess 
quantitatively  and  not  been  examined  under  laboratory  controlled  environments. 

Chemical  transformations  performed  by  certain  bacteria  in  a  media  may  be  important  to 
indoor  air  pollution  and  conservation  efforts.  These  organisms  break  down  the  complex 
organic  molecules.  In  general,  moisture  and  temperature  promote  large  populations  of 
decomposers  and  high  rates  of  metabolism. 

Laboratory  experiments  are  conducted  to  measure  pollutant  flux  from  natural  and 
engineered  multimedia  systems  under  controlled  conditions.  New  test  chambers 
designed  and  fabricated  in  triplicate  mimic  larger  scale  test  chambers  utilized  in  the  field. 
In  the  most  recent  design,  construction,  and  use  of  this  experimental  apparatus;  pollutant 
flux  is  monitored  in  the  laboratory  at  the  multimedia  (i.e.  soil-air)  interface  in  a  controlled 
dynamic  system.  The  laboratory  test  chambers  are  constructed  with  glass  walls  (10  cm 
diameter  X  20  cm  height),  Teflon  top  and  bottom  plates,  viton  o-rings,  and  Teflon  tubing 
and  valves  throughout  the  test  system  to  minimize  wall  interactions  with  the  materials 
and  the  gases  under  consideration.  Pollutant  concentrations  in  the  head  space  above  the 
media  samples  are  monitored  continually  at  one  second  intervals  using  a  Thermo 
Environmental  Instrument  Model  42S  Chemiluminescence  NO-NOi-NOx  analyzer  with 
computer  data  recording.  Comprehensive  physical,  chemical,  and  microbiological  testing 
of  the  samples  is  performed  in  the  laboratory  to  fully  characterize  the  materials  in  the  test 
chamber.  Selected  conditions  during  controlled  experiments  include  temperature,  pH, 
water  filled  pore  space  emd  nitrogen  loadings. 

Soil  temperatures  were  examined  at  5°  C  intervals  between  3°  C  and  48°  C,  for 
various  levels  of  soil  moisture.  The  soil  moisture  is  represented  as  %  water  filled  pore 
space  and  ranged  in  values  from  1 .82  and  96.76.  In  general  NO  flux  was  found  to 
increase  linearly  with  temperature,  and  also  responded  to  varying  amounts  of  soil 
moisture.  For  a  fixed  temperature,  NO  flux  reached  maximum  values  for  %  water  filled 
pore  space  (WFPS)  between  the  range  of  20-45%,  indicating  that  moisture  stress  (low 
levels  of  %WFPS)  or  saturation  (high  levels  of  %WFPS)  may  inhibit  NO  emissions.  NO 
flux  was  also  examined  for  various  levels  of  pH  between  the  ranges  of  4.3  and  8.3.  No 
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significant  trends  were  observed  for  pH  levels  between  5  and  7,  however  NO  was  found 
to  increase  linearly  with  increasing  pH  above  7.3,  and  maximum  NO  emissions  were 
observed  for  acidic  soil  (pH=4.3).  Further  implications  of  these  controlled  studies  to 
indoor  pollution  and  conservation  will  be  discussed. 


Presented  at: 

An  International  Conference  on  Microbiology  and  Conservation  (ICMC  ’99) 
Florence,  Italy 
June  16-19, 1 999 
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Appendix  4.  Seasonal  Variations  of  Nitric  Oxide  Fluxes  from  Diverse  Physiographic 
Agricultural  Soils  in  North  Carolina 

Paul  A.  Roelle,  Viney  P.  Aneja  and  Wayne  Robarge* 

Department  of  Mcirine,  Earth  and  Atmospheric  Sciences 
^Department  of  Soil  Science,  Box  7619 
North  Carolina  State  University,  Box  8208 
Raleigh,  NC  27695 

and  Bruce  Gay,  Thomas  Pierce  and  Chris  Geron 

US  EPA 
MD-80 

Research  Triangle  Park,  NC 

Emissions  of  nitric  oxide  (NO)  were  determined  during  late  spring  and  summer  1995  and 
the  spring  of  1996  from  four  crop  types,  located  at  four  different  physiographic  regions  in 
North  Carolina.  Emission  rates  were  calculated  using  a  dynamic  flow-through  chamber 
system  coupled  to  a  state-of-the-art  mobile  laboratory  for  in-situ  analysis.  Average  NO 
fluxes  during  late  spring  1995  were:  50.9  ±  47.7  ng  N  m’^  s"'  for  com  in  the  lower  coastal 
plain.  Average  NO  fluxes  during  summer  1995  were:  6.4  ±  4.6  and  20.2  ±  19.0  ng  N  m'^ 
s’*  respectively  for  com  and  soybean  in  the  coastal  region;  4.2  +  1.7  ng  N  m’^  s’*  for 
tobacco  in  the  piedmont  region;  and  8.5  ±  4.9  ng  N  m’^  s'*  for  corn  in  the  upper  piedmont 
region.  Average  NO  fluxes  for  spring  1996  were:  66.7  ±  60.7  ng  N  m’^  s'*  for  wheat  in 
the  lower  coastal  plain;  9.5  ±  2.9  ng  N  m’^  s’*  for  wheat  in  the  coastal  plain;  2.7  ±  3.4  ng 
N  m’^  s’*  for  wheat  in  the  piedmont  region;  and  56.1  ±  53.7  ng  N  m’s’*  for  corn  in  the 
upper  piedmont  region.  An  exponential  dependence  of  NO  flux  on  soil  temperature  was 
present  at  all  of  the  locations.  Further,  all  locations  displayed  a  diurnal  trend  of  NO 
emissions,  which  revealed  a  peak  in  NO  emissions  that  coincided  with  the  maximum  soil 
temperature  for  the  day.  The  composite  data  of  all  the  research  sites  revealed  a  general 
positive  trend  of  increasing  NO  flux  with  soil  water  content  and  extractable  nitrogen. 

Presented  at: 

International  Global  Atmospheric  Chemistry  Conference 
Bologna,  Italy 
September  12-17, 1999 
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Appendix  5.  Measurements  and  Analysis  of  Criteria  Pollutants  in  New  Delhi,  India 

Viney  P.  Aneja,  Paul  A.  Roelle,  Sharon  B.  Phillips,  Quansong  Tong,  Nealson  Watkins 
and  Richard  Yablonsky 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

and  A.  Agarwal 

Centre  for  Science  and  Environment 
New  Delhi  110062,  India 

Ambient  concentrations  of  Carbon  Monoxide  (CO),  Nitrogen  Oxides  (NOx),  Sulfur 
Dioxide  (SO2)  and  Total  Suspended  Particulates  (TSP)  were  measured  from  January 
1997  to  November  1998  in  the  center  of  downtown  (the  Income  Tax  Office  (ITO)  located 
on  B.S.G.  Marg)  New  Delhi,  India.  The  data  consists  of  24-hour  averages  of  SO2,  NOx, 
TSP;  and  8  and  24-hour  averages  of  CO.  The  measurements  were  made  in  an  effort  to 
characterize  air  pollution  in  the  urban  environment  of  New  Delhi  and  assist  in  the 
development  of  an  Air  Quality  Index.  The  yearly  average  CO,  NOx,  SO2  and  TSP 
concentrations  for  1997  and  1998  were  found  to  be  4810  +  2287  pg  m'^  and  5772  ±2116 
pg  m'^;  83  ±  35  pg  m'^  and  64  ±  22  pg  m'^;  20  ±  8  pg  m‘^  and  23  ±  7  pg  m'^;  409  ±110 
and  365  ±  100  pg  m'^  respectively.  In  general,  the  maximum  CO,  SO2,  NOx  and  TSP 
values  occurred  during  the  winter  with  minimum  values  occurring  during  the  summer, 
which  can  be  attributed  to  a  combination  of  meteorological  conditions  and  photochemical 
activity  in  the  region.  The  ratio  of  CO  to  NOx  (~50)  indicates  that  mobile  sources  are  the 
predominant  contributors  for  these  two  compounds  in  the  urban  air  pollution  problem  in 
New  Delhi.  The  ratio  of  SO2  to  NOx  (-0.6)  indicates  that  point  sources  are  contributing  to 
SO2  pollution  in  the  city.  The  averaged  background  CO  concentrations  in  New  Delhi 
were  also  calculated  (-1939  pg  m'^)  which  exceed  those  for  Eastern  USA  (~500  pg  m'^). 
Further,  all  measured  concentrations  exceeded  the  US  National  Ambient  Air  Quality 
Standards  (NAAQS)  except  for  SO2.  TSP  was  identified  as  exceeding  the  standard  on  the 
most  frequent  basis. 

Presented  at: 

Workshop  on  Air  Quality  and  Pollution  Inventory  for  the  City  of  Delhi 
New  Delhi,  India 
June  6-8, 2000 
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Appendix  6.  Characterization  of  Nitric  Oxide  Emissions 

from  Sludge  Amended  Soils _ 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

and  Jeffrey  Peirce 

Department  of  Civil  and  Environmental  Engineering 
Duke  University 
Durham,  NC  27708-0287 

Land  spreading  nitrogen-rich  municipal  sludge  to  human  food  and  non-food  chain  land  is 
a  practice  followed  throughout  the  US.  This  practice  may  lead  to  the  recovery  and 
utilization  of  the  nitrogen  by  vegetation,  but  it  may  also  lead  to  emissions  of  biogenie 
nitric  oxide  (NO),  which  may  enhance  ozone  pollution  in  the  lower  levels  of  the 
troposphere.  Recent  global  estimates  of  biogenic  NO  emissions  from  soils  are  cited  in  the 
literature,  which  are  based  on  field  measurements  of  NO  emission  from  various 
agricultural  and  non-agricultural  fields.  However,  biogenie  emissions  of  NO  from  sludge- 
amended  soils  are  lacking.  Using  a  state-of-the-art  mobile  laboratory  and  a  dynamic 
flow-through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were  measured 
during  the  summer  of  1999  from  a  sludge  amended  agricultural  soil.  The  field  site  was 
sampled  both  prior  to  and  immediately  following  an  application  of  the  municipal  waste 
sludge.  The  average  NO  flux  for  the  entire  time  period  prior  to  the  sludge  application 
(June  9  -  September  3)  was  ~  50  ng  N  m'^  s''.  The  results  of  the  research  segregated  into 
time  periods  prior  to  and  immediately  following  the  sludge  application  will  be  presented. 
NO  emissions  responded  to  changes  in  soil  temperature  where  daily  maximums  in  NO 
emissions  coincided  with  daily  maximums  in  soil  temperatures.  Relationships  between 
NO  flux  and  pH  and  percent  water  filled  pore  space  (%WFPS)  will  also  be  presented. 

Presented  at: 

Air  &  Waste  Management  Association’s  93rd  Annual  Meeting  and  Exhibition 
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Appendix  7.  Nitric  Oxide  Emissions  from  Soils  Amended  with 
Municipal-waste  Biosolids 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

and  Jeffrey  Peirce 

Department  of  Civil  and  Environmental  Engineering 
Duke  University 
Durham,  NC  27708-0287 

Land  spreading  nitrogen-rich  municipal  waste  biosolids  (N-NOs’  ~  443  mg/Kg  Dry 
biosolid,  N-NO2'  ~  443  mg/Kg  dry  biosolid,  Phosphorus  ~  11,970  mg/Kg  Dry  biosolid. 
Total  N  ~  42,586  mg/Kg  Dry  biosolid)  to  human  food  and  non-food  chain  land  is  a 
practice  followed  throughout  the  US.  This  practice  may  lead  to  the  recovery  and 
utilization  of  the  nitrogen  by  vegetation,  but  it  may  also  lead  to  emissions  of  biogenic 
nitric  oxide  (NO),  which  may  enhance  ozone  pollution  in  the  lower  levels  of  the 
troposphere.  Recent  global  estimates  of  biogenic  NO  emissions  from  soils  are  cited  in  the 
literature,  which  are  based  on  field  measurements  of  NO  emission  from  various 
agricultural  and  non-agricultural  fields.  However,  biogenic  emissions  of  NO  from  soils 
amended  with  biosolids  are  lacking.  Using  a  state-of-the-art  mobile  laboratory  and  a 
dynamic  flow-through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were 
measured  during  the  Summer/Fall  of  1999  and  Winter/Spring  of  2000  from  an 
agricultural  soil  which  is  routinely  amended  with  municipal  waste  biosolids.  The  average 
NO  flux  for  the  entire  time  period  prior  to  a  biosolids  application  (June  9  -  September  3) 
was  57.8  ±  34.6  ng  N-NO  m'^  s''.  Field  experiments  were  conducted  which  indicated  that 
the  application  of  biosolids  increase  the  emissions  of  NO.  The  results  of  the  research 
segregated  into  time  periods  prior  to  and  immediately  following  the  biosolids  application 
will  be  presented.  NO  emissions  responded  to  changes  in  soil  temperature  where  daily 
maximums  in  NO  emissions  coincided  with  daily  maximums  in  soil  temperatures. 
Relationships  between  NO  flux  and  pH  and  percent  water  filled  pore  space  (%WFPS) 
will  also  be  presented. 

Presented  at: 
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Appendix  8.  Measurement  and  Modeling  of  Biogeochemical 

Cycling  of  Nitrogen  Compounds _ 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Land  spreading  nitrogen-rich  municipal  sludge  to  human  food  and  non-food  chain  land  is 
a  practice  followed  throughout  the  US.  This  practice  may  lead  to  the  recovery  and 
utilization  of  the  nitrogen  by  vegetation,  but  it  may  also  lead  to  emissions  of  biogenic 
nitric  oxide  (NO),  which  may  enhance  ozone  pollution  in  the  lower  levels  of  the 
troposphere.  Recent  global  estimates  of  biogenic  NO  emissions  from  soils  are  cited  in  the 
literature,  which  are  based  on  field  measurements  of  NO  emission  from  various 
agricultural  and  non-agricultural  fields.  However,  biogenic  emissions  of  NO  from  sludge- 
amended  soils  are  lacking.  Using  a  state-of-the-art  mobile  laboratory  and  a  dynamic 
flow-through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were  measured 
during  the  summer  of  1999  from  a  sludge  amended  agricultural  soil.  The  field  site  was 
sampled  both  prior  to  and  immediately  following  an  application  of  the  municipal  waste 
sludge.  The  average  NO  flux  for  the  entire  time  period  prior  to  the  sludge  application 
(June  9  -  September  3)  was  ~  50  ng  N  m'^  s'^  The  results  of  the  research  segregated  into 
time  periods  prior  to  and  immediately  following  the  sludge  application  will  be  presented. 
NO  emissions  responded  to  changes  in  soil  temperature  where  daily  maximums  in  NO 
emissions  coincided  with  daily  maximums  in  soil  temperatures.  Relationships  between 
NO  flux  and  pH  and  percent  water  filled  pore  space  (%WFPS)  will  also  be  presented. 

Presented  at: 

Meteorology  at  the  Millennium,  150th  Anniversary  Conference 
Cambridge,  United  Kingdom 
July  10-14,  2000 
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Appendix  9.  Measurement  and  Modeling  of  Nitric  Oxide  Emissions  from  Biosolid 
Amended  Soils 


Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

and  Sonia  Aneja 
Summer  Intern 
Broughton  High  School 
Raleigh,  NC 

Land  spreading  nitrogen-rich  municipal  biosolids  to  hiunan  food  and  non-food  chain  land 
is  a  practice  followed  throughout  the  US.  This  practice  may  lead  to  the  recovery  and 
utilization  of  the  nitrogen  by  vegetation,  but  it  may  also  lead  to  emissions  of  biogenic 
nitric  oxide  (NO),  which  may  enhance  ozone  pollution  in  the  lower  levels  of  the 
troposphere.  Recent  global  estimates  of  biogenic  NO  emissions  from  soils  are  cited  in  the 
literature,  which  are  based  on  field  measurements  of  NO  emission  from  various 
agricultural  and  non-agricultural  fields.  However,  biogenic  emissions  of  NO  fi’om 
biosolid-amended  soils  are  lacking.  Using  a  state-of-the-art  mobile  laboratory  and  a 
dynamic  flow-through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were 
measured  during  the  summer,  winter  and  spring  of  1999/2000  from  agricultural  soils 
amended  with  biosolids.  The  field  site  was  sampled  prior  to  and  immediately  following 
applications  of  the  biosolids.  The  average  NO  flux  for  summer,  winter  and  spring  were 
57.8  ±  34.6  ng  NO-N  s\  3.62  ±  3.93  ng  NO-N  m'^  s'‘  and  59.05  ±  54.21  ng  NO-N  m' 
^  s'*  respectively.  NO  emissions  responded  to  changes  in  soil  temperature  where  daily 
maximums  in  NO  emissions  coincided  with  daily  maximums  in  soil  temperatures. 

Presented  at: 

International  Symposium  on  the  Measurement  of  Toxic  and  Related  Air  Pollutants 
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Appendix  10.  Effect  of  Environmental  Variables  on  NO  Emissions 
from  Agricultural  Soils 


Viney  P.  Aneja,  Paul  A.  Roelle  and  Yongxian  Li 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Nitric  oxide  (NO)  is  an  important  air  pollutant  which  leads  to  the  production  of  ozone 
and  acidic  precipitation.  Biogenic  emissions  of  NO  are  gaining  in  importance  for  ozone 
formation  in  semi-urban  and  rural  regions  of  the  Southeast  United  States.  Using  a 
dynamic  chamber  system  interfaced  to  a  mobile  laboratory  for  continuous  NO  analysis, 
soil  emissions  of  NO  were  measured  over  typical  row  crops  in  North  Carolina  during 
1994  -  1996.  This  paper  investigates  the  effect  of  soil  temperature,  soil  water  content  and 
soil  extractable  nitrogen  on  soil  NO  emissions  from  com  and  soybean  canopies.  Three  of 
five  sets  of  measurements  from  com  and  two  of  four  sets  of  measurements  from  soybean 
displayed  an  exponential  relationship  between  soil  NO  emissions  and  soil  temperature. 
No  significant  correlation  between  soil  NO  emissions  and  water  filled  pore  space  (WFPS) 
was  observed.  The  best  correlation  observed  between  soil  NO  emissions  and  total 
extractable  nitrogen  (TEN)  was  found  to  be  linear  for  soybean  (NO  emission  =  0.67  + 
1.43*TEN,  R^  =  0.34),  and  logarithmic  for  com  (NO  emission  =  21.20*Ln  (TEN)  - 
27.27,  R^  =  0.17).  Multiple  regression  models  (Ln  (NO  emission)  =  1.5017  + 
0.0786*TEN  -  0.0006*TEN^;  R^  =  0.58  for  com;  Ln  (NO  emission)  =  77.917  -  6.19*T 
+  0.1243*T^  +  0.0068*TEN^;  R^  =  0.65  for  soybean)  appear  to  be  the  best  overall  model 
for  predicting  soil  NO  emissions. 

Presented  at: 

Internet  Conference:  Nitrogen  Emissions  from  Soil 

“Programme  based  on  an  Austrian  research  co-operation  between  the  Magistrat  der  Stadt 
Wien,  MA  49  (City  of  Vienna)  and  the  Bundesministerium  fur  Bildung,  Wissenschaft 
und  Kultur  (Austrian  Federal  Ministry  of  Education,  Science  and  Culture) 
<www.nitro-soil.at> 

September  1 8  -  December  3 1 ,  2000 
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Appendix  11.  Characterization  of  Sources  of  Biogenic  Atmospheric  Ammonia 

Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Atmospheric  ammonia  (NH3)  has  recently  reached  a  heightened  awareness  as  a  result  of 
its  environmental  impact  on  particulate  matter  formation,  soil  acidification,  aquatic 
eutrophication  and  odor  issues.  Animal  waste  (22  Tg  NH3  yf’)  contributes  the  largest 
emissions  to  the  total  global  NH3  budget  (45  Tg  NH3  year'')  (Dentener  and  Crutzen, 
1994),  hence  areas  of  intensive  livestock  operations  are  potential  sources  of  significant 
ammonia  emissions.  North  Carolina  is  currently  the  second  largest  hog  producing  state  in 
the  USA  and  the  current  technology  to  manage  the  waste  is  known  as  the  lagoon  and 
spray  system,  which  consists  of  an  exposed  waste  lagoon  and  mechanisms  through  which 
the  waste  is  periodically  sprayed  onto  the  soil  as  a  nutrient  source.  It  has  been  reported 
that  a  substantial  fraction  of  the  total  NH3  emissions  in  NC  are  from  these  intensive 
livestock  operations  (9.21  kg  NH3  animal''  year'')  (Aneja  et  ah,  1998;  Battye  et  ah, 
1994).  The  objective  of  this  study  is  to  characterize  the  sources  of  biogenic  atmospheric 
ammonia.  Fluxes  of  NH3  were  determined  in-situ  using  a  dynamic  flow-through  chamber 
system,  in  conjunction  with  a  high  sensitivity  chemiluminescence  ammonia  analyzer.  Soil 
and  water  samples  from  the  center  of  the  chamber  footprint  were  collected  and  analyzed 
for  pH,  Total  Kjeldahl  Nitrogen  (TKN)  content,  and  percent  water  filled  pore  space. 
Results  indicate  that  NH3  fluxes  from  the  agricultural  field  sites  ranged  from  45-  100  ng 
N  m'^  s''  for  soils  and  305-4017  pg  N  m'^  min''  for  lagoons.  NH3  emissions  responded  to 
changes  in  temperature  where  daily  maximums  in  NH3  emissions  coincided  with  daily 
maximums  in  temperatures.  Relationships  between  NH3  flux,  temperature,  pH  and  TKN 
will  also  be  presented. 

Presented  at; 
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Appendix  12.  Nitric  Oxide  Emissions  from  Biosolid  Amended  Soils 


Paul  A.  Roelle  and  Viney  P.  Aneja 

Department  of  Marine,  Earth  and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27695-8208 

Land  spreading  nitrogen-rich  municipal  biosolids  to  human  food  and  non-food  chain  land 
is  a  practice  follov^^ed  throughout  the  US.  This  practice  may  lead  to  the  recovery  and 
utilization  of  the  nitrogen  by  vegetation,  but  it  may  also  lead  to  emissions  of  biogenic 
nitric  oxide  (NO),  which  may  enhance  ozone  pollution  in  the  lower  levels  of  the 
troposphere.  Recent  global  estimates  of  biogenic  NO  emissions  from  soils  are  cited  in  the 
literature,  which  are  based  on  field  measurements  of  NO  emission  from  various 
agricultural  and  non-agricultural  fields.  However,  biogenic  emissions  of  NO  from 
biosolid-amended  soils  are  lacking.  Using  a  state-of-the-art  mobile  laboratory  and  a 
dynamic  flow-through  chamber  system,  in-situ  concentrations  of  nitric  oxide  (NO)  were 
measured  during  the  summer,  winter  and  spring  of  1999/2000  from  agricultural  soils 
amended  with  biosolids.  The  field  site  was  sampled  both  prior  to  and  immediately 
following  an  application  of  the  biosolids.  The  average  NO  flux  for  the  entire  time  period 
prior  to  the  sludge  application  (June  9  -  September  3)  was  ~  50  ng  N  m‘^  s’*.  The  results 
of  the  research  segregated  into  time  periods  prior  to  and  immediately  following  the 
sludge  application  will  be  presented.  NO  emissions  responded  to  changes  in  soil 
temperature  where  daily  maximums  in  NO  emissions  coincided  with  daily  maximums  in 
soil  temperatures.  Relationships  between  NO  fiux  and  pH  and  percent  water  filled  pore 
space  (%WFPS)  will  also  be  presented. 

Presented  at: 

2000  Annual  Meeting  and  International  Conference  of  the  American  Institute  of 
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